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PREFACE 


The Institute for Defense Analyses (IDA) prepared this paper for the Defense 
Threat Reduction Agency (DTRA), in partial fulfillment of the task “Support for DTRA 
in the Validation Analysis of Hazardous Material Transport and Dispersion Prediction 
Models.” The objective of this effort was to conduct analyses and special studies 
associated with the verification, validation, and accreditation of hazardous transport and 
dispersion prediction models. This task involves the comparison of two models: DTRA’s 
Hazard Prediction and Assessment Capability (HPAC) and Lawrence Livermore National 
Laboratory’s (LLNL) National Atmospheric Release Advisory Center (NARAC) 
modeling system. 

The IDA Technical Review Committee was chaired by Thomas P. Christie and 
consisted of Arthur Fries, Jeffrey H. Grotte, Ira Kohlberg, David A. McWhorter, and 
Jozsef A. Toth. The authors would like to thank Allan Reiter (DTRA), Leon Wittwer 
(DTRA), Ian Sykes (ARAP Group, Titan Corporation), and Don Ermak (LLNL) for their 
comments, critiques, and support throughout this effort. 

The LLNL work was performed under the auspices of the U.S. Department of 
Energy by University of California Lawrence Livermore National Laboratory under 
contract No. W-7405-Eng-48. 
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SUMMARY 


A. BACKGROUND 

During fiscal year 2000, a series of studies in support of the Defense Threat 
Reduction Agency (DTRA) was begun. The goal of these studies is to improve the 
verification, validation, and accreditation (VV&A) of hazard prediction and assessment 
models and capabilities. These studies are part of a larger joint VV&A collaborative 
effort that DTRA and the Department of Energy (DOE), via the Lawrence Livermore 
National Laboratory (LLNL), are conducting. This joint effort includes comparisons of 
the LLNL and DTRA transport and dispersion (T&D) modeling systems, NARAC and 
HP AC, respectively. 1 

B. PURPOSE 

The purpose of this work is to compare, in a systematic way, HPAC and NARAC 
model predictions for a set of controlled hypothetical release scenarios. Only “model- 
versus-model” comparisons are addressed in this work. Model-to-field trial comparisons 
for HPAC and NARAC have been addressed in a recent companion study, 2 in support of 
the same joint VV&A effort. 

C. METHODOLOGY 
1. Compared Releases 

In January 2000, members of the joint collaboration met at LLNL to establish a 
set of release scenarios for this model-to-model study. Seventeen scenarios were agreed 
upon based on several selection criteria. 


1 NARAC = National Atmospheric Release Advisory Center and HPAC = Hazard Prediction and 
Assessment Capability. 

2 Warner S., Platt, N., Heagy, J. F., Bradley, S., Bieberbach, G., Sugiyama, G., Nasstrom, J. S., Foster, 
K. T, and Larson, D., User-Oriented Measures of Effectiveness for the Evaluation of Transport and 
Dispersion Models, IDA Paper P-3554, 8 January 2001. 
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First and foremost was the desire for simplicity so that the scenarios would test 
only the core transport and dispersion algorithms within each model. Consequently, all 
releases were chosen to have idealized weather - a defined single unidirectional vertical 
wind speed profile. 3 Terrain was similarly elementary - flat terrain with a constant 
surface roughness of 0.008 m. Downwind transport distances of 40 km were used for all 
but one “high altitude” release. Initial source geometries were also kept simple - 
spherical sources (Gaussian with standard deviation = 1 m in all directions) for 16 of the 
17 releases and a 500 m line source in the remaining release. 

A mix of atmospheric conditions - stable, near-neutral, and unstable - were 
chosen through suitable choices of the atmospheric boundary layer height, z l7 and Monin- 
Obukhov length, L. A range of release heights (relative to the boundary layer height) 
was also chosen. 4 

Two agent types were considered: “neutrally-buoyant gas” (9 releases) and 
“particles” (8 releases). Particle runs were carried out for three particle size distributions, 
two using fixed diameter particles of 5 microns (p) and 50 p, and the third with a 
log-normal distribution with a mass median diameter (MMD) of 500 p and a geometric 
standard deviation (a) of 2. Releases were chosen to be either instantaneous or 
continuous, and in all cases the mass of released agent was 1 kg. 5 

The chosen comparison scenarios do not represent an exhaustive set of potential 
cases and should not be considered as the basis for a complete statistical analysis of 
model comparisons. Rather, the trials were chosen to examine the impact of critical 
transport and dispersion parameters on relative model behavior. Specifically, the 17 trials 
allowed for the comparison of model predictions as a function of atmospheric stability 
condition, source type (gas or particle), release duration, release height, and particle size. 6 


3 Vertical wind speed variations were inferred from a 5 m/s wind speed observation at 10 meters. See 
Chapter 2, Section B.l for additional details. 

4 The Monin-Obukhov length, L, is a length scale determined from the surface heat and momentum 
fluxes, and is a fundamental scaling parameter for boundary layer turbulence. The boundary layer 
height defines the depth of the turbulent layer generated by the surface, and is another fundamental 
turbulence scaling parameter. See, Pasquill, F., Atmospheric Diffusion, The Dispersion of Windborne 
Material from Industrial and other Sources, Second Edition, Wiley, 1974. 

5 Additional details associated with the examined releases are provided in Chapter 2. Note that these 
releases were not intended to represent actual or hypothetical accidents or attacks. 

6 In order to manage the size and scope of this study, only a limited range of values were chosen for 
several of these variables. 
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Table 1 presents the key input conditions for the seventeen runs that were 
examined. The first column denotes the run number and the second column identifies the 
source type (neutrally-buoyant gas or particles). The next two columns present the values 
used for two critical dimensionless input parameters - the ratio of release height to 
boundary layer height (z r /zi) and the ratio of boundary layer height to Monin-Obukhov 
length (zj/L). 


Table 1. Key Input Conditions for the Model-to Model Comparisons 


Run 

Source 

Release Height / 
Boundary Layer 
Height (z r /Zj) a 

Boundary Layer Height 
/ Monin-Obukhov 
Length (z s /L) b 

1 

neutrally-buoyant gas 

0.02 

4 (stable) 

2 

neutrally-buoyant gas 

0.50 

4 (stable) 

3 

neutrally-buoyant gas 

0.80 

4 (stable) 

4 

neutrally-buoyant gas 

0.004 

-1 (near-neutral) 

5 

neutrally-buoyant gas 

0.50 

-1 (near-neutral) 

6 

neutrally-buoyant gas 

0.80 

-1 (near-neutral) 

7 

neutrally-buoyant gas 

0.004 

-10 (Unstable) 

8 

neutrally-buoyant gas 

0.50 

-10 (Unstable) 

9 

neutrally-buoyant gas 

0.80 

-10 (Unstable) 

10 

5 [i 

0.004 

-1 (near-neutral) 

11 

5 [i 

0.50 

-1 (near-neutral) 

12 

5 [i 

0.80 

-1 (near-neutral) 

13 

50p 

0.004 

-1 (near-neutral) 

14 

50p 

0.004 

-1 (near-neutral) 

15 

50p 

1.50 

-1 (near-neutral) 

16 

log normal 0 

0.004 

-1 (near-neutral) 

17 

log normal 

0.004 

-1 (near-neutral) 


a z, = release height, z ; = boundary layer height 

b L = Monin-Obukhov length. For this study, atmospheric stability categories were defined as follows: 
stable runs Zi/L = 4, near-neutral runs Zj/L = -1, and unstable runs Zi/L = -10. 

c For the log-normal particle distribution, the MMD = 500 |i and the geometric standard deviation = a = 
2 . 

2. Measures Used for These Comparisons 

For each run, we compared contoured regions at three to seven dosage levels. 
Dosage contours provide immediate graphical information about the relative performance 
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of model-to-model results; model agreement or lack thereof can be quickly assessed. The 
following quantitative measures were also used for comparisons of model outputs. 7 

. Downwind Contour Distances: For a given near-surface (10m) dosage value, 
the downwind distance (along the symmetry axis) to that dosage was 
computed for both models. This metric, in particular, addresses the relative 
transport characteristics of the model predictions. 

. Crosswind Dosage Plume Widths: For these comparisons, we used a dose- 
weighted width (described in Chapter 2). This metric allows for direct 
comparisons of model crosswind dispersion features. 

. Area-Based Measure of Effectiveness (MOE): This measure considers three 
regions of interest: 1) HPAC and NARAC predictions agree (overlap), 2) 
HPAC predicts a larger region for a given dosage than NARAC, and 3) 
NARAC predicts a larger region for a given dosage than HPAC. 8 A two- 
dimensional vector composed of the fraction of the NARAC prediction that 
corresponds to overlap and the fraction of the HPAC prediction that 
corresponds to overlap represents the area-based MOE. 9 The MOE provides 
information associated with both transport (direction and downwind distance) 
and dispersion. 

D. FINDINGS 

1. Model Intercomparisons 

One of the chief findings of this study was that, with the selection of consistent 
input parameters (see Section D.2), the predictions of HPAC and NARAC for many of 
these simple scenarios agreed quite well. Figure 1 displays three such cases. Shown in 
the figure are the predicted HPAC (blue) and NARAC (red) near-surface (10m) dosage 
contours for three neutrally-buoyant gas release runs: Model versus Model run 1 (MvM 
1), MvM 4, and MvM 7. These three runs, each with release height 2 meters (m), 
correspond to three meteorological stability categories, stable, near-neutral, and unstable, 
respectively. 


7 In addition, to the metrics described here, we examined the crosswind maximum dosages (i.e., the 
maximum predicted dosages of each model, at several downwind distances) and the dosage along the 
downwind symmetry axis. 

8 For simplicity when discussing the area-based MOE, we will use the terminology 
“overprediction/underprediction” when one model predicts a larger/smaller area for a specified contour 
level. 

9 Chapter 2, Section C.5 provides additional discussion. 
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MvMOl: Dosage (kg s/nr), Height = 10 m, Time = 60 min. 


Stable 


Near- 

Neutral 





Downwind (km) 

Figure 1. Results for Three Nominal Gas Release Scenarios as a Function of 
Meteorological Stability Category: at 60 Minutes with 
NARAC (Red -) and HPAC (Blue - -) Dosage Contours 


Close agreement could also be found in many of the particle releases. Figure 2 
shows near-surface (10m) dosage contours for two particle releases, one with 5 |i 
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particles and the other with 50 p particles. Each release occurred at a height of 2 meters. 
Good agreement is seen in both cases, with better agreement for the 5 p particle release 
(relative to the 50 p particle release) at least at the higher dosages 


MvMI 0: Dosage (kg s/m' 3 ), Height = 10m,Time = 60 min. 


5 p 


50 




Figure 2. Results for Two Nominal Particle Release Scenarios as a Function of Particle 
Size (5 and 50 p): at 60 Minutes with NARAC (Red -) and HPAC (Blue - -) Dosage Contours 


a. Influence of Release Height on Comparisons 

A significant finding of this study was that release height (relative to boundary 
layer height) had dramatic effects on model agreement, with a clear trend for better 
model agreement with lower release heights. The highest height releases displayed 
marked differences between model predictions. Figure 3 compares HPAC and NARAC 
surface dosage contours for MvM 3 at 30 and 60 minutes after the release (80 meter gas 
release height with a boundary layer height of 100 meters) and for MvM 15 at 120 and 
240 minutes (750 meter release height with a boundary layer height of 500 meters). 
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For MvM 3, the HP AC 10' kg-s/m 30-minute dosage contour completely 
encloses the NARAC 10' 13 kg-s/m 3 contour, indicating relative overprediction by 
HP AC. 10 This trend continues at 60 minutes. Similar conclusions are seen for MvM 15 
at 120 and 240 minutes. 


MvM03: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. MvM03: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 



MvM15: Dosage (kg s/m 3 ), Height = 10 m, Time = 120 min. MvM15: Dosage (kg s/m 3 ), Height = 10 m, Time = 180 min. 



Figure 3. Results for Two “Higher” Altitude Release Scenarios: MvM 3 at 30 and 60 
Minutes and MvM 15 at 120 and 240 Minutes with NARAC (Red -) and HPAC (Blue - -) 

Dosage Contours 

Differences in the vertical diffusivity models used, as well as in the blending of 
the diffusivity across the boundary layer, are likely contributors to the marked differences 
in the predictions for these higher altitude releases. 

Figure 4 displays the area-based MOE for all 17 compared releases for various 
contour levels. In the figure, the red circles denote releases with 2-meter release heights 
(the lowest release height considered), and blue triangles denote all other releases (with 
release heights ranging from 50 to 750 meters). Each symbol is associated with a single 


Equivalently, NARAC underpredicts HPAC. Throughout this study when drawing comparisons we 
adopt the nomenclature that one model overpredicts the other. This terminology is not to be construed 
as a statement about the essential correctness or incorrectness of either model. 
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near-surface (10 m) contour level region (therefore, there are more symbols than runs). 
For each run, between three and seven contour levels were examined. 



Figure 4. Area-Based MOE for All Gas and Particle Releases. Red Circles Denote 
Releases at 2-Meter Height (Lowest Height Considered) and Blue Triangles Denote 
Releases At All Other Heights. MOE Values Separate as a Function of Height. A OL = the 
Area of Overlap, A narac = the Area of the NARAC Prediction, and A HPA c = the Area of the 
HPAC Prediction All at a Specific Contour Level. 

The symbols in this plot show a very definite separation, with the 2-meter releases 
falling mainly above the diagonal (NARAC overprediction relative to HPAC) and the 
higher altitude releases falling mainly below the diagonal (HPAC overprediction relative 
to NARAC). Extreme regions where NARAC contours enclose HPAC contours (y = 1) 
and where HPAC contours enclose NARAC contours (x = 1) are marked. There are a 
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number of points that lie close to the diagonal and near the point x = y = 1 (the point of 
“perfect” model agreement, with respect to this measure). These MOE comparisons are 
consistent with the conclusions drawn from the graphical comparisons discussed 
previously. 

b. Influence of Atmospheric Stability on Gas Release Comparisons 

Closer inspection of the three contour plots for neutrally-buoyant gas releases 
shown in Figure 1 reveals that model agreement improves slightly with decreasing 
meteorological stability. Contour differences are also most pronounced at the longer 
ranges, where the trend is for longer (downwind) and wider (crosswind) HP AC contours. 
Figure 5 shows a scatter plot of the near-surface downwind contour distances for the gas 
releases at 60 minutes, with color indicating stability condition (blue = stable, red = 
neutral, and green = unstable). 

Overall, agreement is good, with the bulk of the points falling within the 15 
percent overprediction limits. Taken by groups, the unstable (green) points are closest to 
the diagonal line, followed by the neutral (red) points, then the stable (blue) points, 
reflecting the graphical differences displayed in Figure 1. Other comparison metrics 
(e.g., crosswind plume widths and area-based MOEs) show a similar trend. 

For the three gas cases examined with a source release height of 2 meters, the 
relative 1 km downwind crosswind plume widths for the HP AC and NARAC predictions 
are shown in Figure 6. NARAC relatively overpredicts HP AC crosswind plume width 
for stable and near-neutral conditions. For the unstable case, MvM 7, the 1 km crosswind 
dosage plume widths are similar. This is consistent with the conclusions drawn from the 
Prairie Grass comparisons that involved similar downwind distances. 11 

Differences between models for the runs done with stable meteorological 
conditions may not be too surprising, since it is for these conditions, relative to neutral 
and unstable, that the understanding of transport and dispersion phenomena is least 
complete. 


11 Warner S., Platt, N., Heagy, J. F., Bradley, S., Bieberbach, G., Sugiyama, G., Nasstrom, J. S., Foster, 
K. T, and Larson, D., User-Oriented Measures of Effectiveness for the Evaluation of Transport and 
Dispersion Models, IDA Paper P-3554, 8 January 2001, pp. 3-31 and pp. 1-18, Figure 1-18, e.g., see the 
“800 meter arc, 95%” chart. 
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Gas Runs: Downwind Contour Distances, 60 Minutes 



NARAC Dosage Contour Distance (km) 


Figure 5. Scatter Plot Comparison of Downwind Contour Distances at 60 Minutes for All 
Gas Runs: Colors Indicate Meteorological Stability Categories and Dashed Lines are 15% 

Overprediction Limits 
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x 10 5 MvMOI : Crosswind Dosage Profiles, Time = 30 min. 



Figure 6. Dosage Contours at 30 Minutes and 1 km Downwind for Three 2-Meter Height 
Neutrally-Buoyant Gas Releases: MvM 1 (Stable), MvM 4 (Near-Neutral), and MvM 7 
(Unstable). Of the Seventeen Releases Studied, These Are Considered the Most Similar to 
the Previously Examined Prairie Grass Field Trials. 
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c. Influence of Particle Size on Comparisons 

Figure 7 presents a scatter plot of HP AC and NARAC near-surface downwind 
contour distances (60 minutes after the release) with particle size indicated by color (blue 
= 5 |t, red = 50 p, and green = 500 |i). Generally, HPAC and NARAC downwind 
contour distances agree to within 20 percent. HPAC slightly overpredicts NARAC for 
this measure at the longer ranges/lower dosages. If all distance/dosages are considered, 
the 5 |i particle results show the fewest large differences. For the shorter ranges/higher 
dosages (inside 20 km) the 5 |i points show the best agreement. 

For the 50 |i (red) and 500 |i (green) particles, the disagreement between the two 
models is greatest for the shorter distances (higher dosages) with NARAC relatively 
overpredicting HPAC (e.g., see the 10~ 7 dosage contours for MvM 13 in Figure 2). For 
longer distances (lower dosages), there is a slight opposite trend, with the 500 |i points 
(green) lying closest to the diagonal and the 5 |i points (blue) lying farthest from the 
diagonal. 

The particle settling models employed by HPAC and NARAC are somewhat 
different (there are several approaches to the settling of particles in the transport and 
dispersion literature). It is therefore reasonable to expect some differences in model 
performance as a function of particle size, such as those observed in this study. 

2. HPAC Prediction Excursions: Resolution Issues 

In the initial stages of this study, a set of “baseline” parameters that control spatial 
and temporal resolution was agreed to after a series of test comparisons. During these 
explorations, some significant resolution-dependent effects were noted. These effects 
were used to guide the final choice for the baseline resolution settings. Choice of proper 
resolution is always critical to model simulations. The discussion in this section is 
intended to underscore the need to thoroughly understand and account for the intrinsic 
computational framework of the models in a model-to-model comparison. Three main 
effects were observed, two concerning HPAC spatial resolution and one concerning 
HPAC temporal resolution. 12 


12 Related considerations apply to setting up NARAC meteorological and concentration grids but are not 
discussed here. (See Chapter 2 for additional discussion.) 
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Particle Runs: Downwind Contour Distances, 60 Minutes 



Figure 7. Scatter Plot Comparison of Downwind Contour Distances at 60 Minutes for 
Particle Runs MvM 10 through 14,16, and 17: Colors Indicate Particle Size and Dashed 

Lines are 20% Overprediction Limits 

a. Vertical Resolution 

Substantial changes in HP AC downwind plume transport were observed when the 
HP AC vertical resolution was changed. This effect was most prominent in the 2-meter 
stable gas release, MvM 1. Figure 8 shows dosage contour comparisons with two values 
of the HP AC vertical resolution: a “coarse” value of 40 meters (the default HP AC value 
for a grid top height of 2,000 meters) and a “fine” value of 10 meters (the adopted 
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baseline setting). The coarse resolution HP AC plume (blue) propagates much farther 
downwind than the fine resolution HP AC plume. 

This effect is caused by the inadequate resolution of the vertical wind profile for 
the coarse case, resulting in an artificially high transport speed at elevations near the 
release height of 2 meters. This effect motivated the choice to adopt the finer vertical 
resolution value of 10 meters for all of the releases. Retention of the default HP AC 
resolution value of 40 meters would have clearly (and inappropriately) skewed the 
conclusions of this study 


Vertical Resolution 

^ A 


Coarse (40 meters) Fine (10 meters) 

MvMOl: Dosage (kg s/m 3 ), Height = 10m, Time = 60 min. MvMOl: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure 8. MvM 1 at 60 Minutes with Upper: Dosage Contours for HPAC (Blue - -) and 
NARAC (Red — ); Lower: Dosage on Downwind Symmetry Axis. HPAC Runs Are Done at 
Coarse and Fine Vertical Resolution Values of 40 and 10 Meters as Noted. 
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b. Horizontal Resolution 


A second resolution-induced effect was observed in HPAC dosage output, when 
the horizontal resolution was changed from the default value of approximately 1 km to 
the finer value of 0.4 km. 13 Two differences between the contour sets done with different 
horizontal resolution were apparent. First, the downwind plume extent for the fine (i.e., 
higher resolution) contours was somewhat reduced. Next, the finer resolution setting also 
had the undesirable effect of introducing high-frequency spatial modulation in the 
contours that distorted the plume shape, particularly at larger dosages. We hypothesize 
that, computational artifacts introduced by “puff-splitting” routines in SCIPUFF 14 that 
control the creation of new Gaussian concentration puffs when plumes cross resolution 
cells, cause these modulations. 

We adopted the default HPAC horizontal resolution value of 1 km as the baseline 
value in all model-to-model comparisons. This choice, combined with the 10-meter 
vertical resolution choice, produced HPAC contours that had comparable transport 
speeds and acceptable (i.e., sensible) smoothness properties. Contrary to the vertical 
resolution case, retention of the default HPAC resolution value in the horizontal case was 
required to ensure a reasonable comparison. 

c. Time-Step Resolution 

Significant spatial oscillations in the HPAC dosage were observed near the 
release for all of the 2-meter instantaneous releases (6 out of the 17 releases). These 
oscillations eventually decay with downwind distance (~ 1 to 3 km downwind). While 
the precise cause for the oscillations is unclear (as is their absence in the 11 other 
releases), we have demonstrated that the magnitude of the oscillations is related to the 
output time-step. Figure 9 shows the symmetry axis dosage for MvM 4 at 30 minutes, 
computed with two time-steps: 60 seconds (the baseline value, in blue) and 30 seconds 
(in red). Decreasing the time-step by a factor of two diminishes the magnitude of the 
oscillations by as much as a factor of 100. 15 


* 3 It was initially thought that finer horizontal resolution would result in more favorable comparisons. 

14 SCIPUFF = Second-Order Closure Integrated Puff. SCIPUFF is the set of transport and dispersion 
algorithms used in the version of HPAC that we examined. 

15 An explanation for this effect is given in Chapter 3, Section C.3. We note that the above oscillations 
are a consequence of our particular comparative methodology and stress that the internally computed 
HPAC dosages, that are typically used to create a dosage plot, do not show these oscillations near the 
source. 
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Because we did not focus on short-range dosage comparisons, the presence of the 
dosage oscillations did not adversely affect the comparisons. Future model-to-model 
studies that involve comparisons of short-range HP AC dosages would likely benefit from 
further reductions in the time-step. 


MvM04, 30 Min: Dosage on Downwind Symmetry Axis 



Figure 9. HPAC Dosage Along Downwind Symmetry Axis for MvM 4 at 30 Minutes. Blue 
Trace Is for 60 Second Time-Step (Baseline Value) and Red Trace Is for 30 Second Time- 
Step. The Magnitude of the Dosage Oscillations Decreases Substantially by Decreasing 

the Time-Step. 


3. Conclusions 

Controlled comparisons of HPAC and NARAC predictions of 17 notional releases 
were completed. In general, good agreement was obtained between the models, with the 
majority of comparisons resulting in predicted downwind contour distances and 
crosswind plume widths within 15 percent and 50 percent of each other, respectively. 

To achieve close agreement between many of the HPAC and NARAC runs, we 
implemented a stringent protocol for this exercise (e.g., identical fixed wind profiles. 
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careful control of all physical input parameters and model-specific parameters controlling 
various resolution aspects). Close agreement also required several iterations of model 
predictions in order to arrive at a set of baseline input parameter settings for all of the 
runs. These baseline settings were adopted and agreed to by all members of the 
collaboration. 

HP AC and NARAC take substantially different mathematical and numerical 
approaches to transport and dispersion of material: HP AC uses a second-order closure 
equation set and a Gaussian puff method, while NARAC uses the diffusion equation and 
a Lagrangian, Monte Carlo particle method. Thus, the fact that the core transport and 
dispersion subsystems of HP AC and NARAC agree closely on a range of releases should 
inspire users with a certain degree of confidence in physical validity of both of these 
models. 

The scenario releases that we examined allowed us to probe the effects on 
differences between model predictions of release height, particle size, meteorological 
stability category, and to a lesser degree, source term geometry and duration of the 
release (“instantaneous” versus “continuous”). We found that major differences between 
model predictions occurred for greater release heights with increasing atmospheric 
stability and larger particle size contributing as secondary factors. In particular, runs 3 
and 15 (Table 1), with z r /zi values of 0.80 and 1.50, demonstrated dramatic differences in 
model predictions (due to vertical diffusivity differences, which were associated with 
release height and boundary layer changes). Runs 1 through 9 were used to explore 
differences in model predictions as a function of atmospheric stability. Runs 10 through 
14, 16, and 17 allowed for the examination of model prediction differences related to 
particle size. 

E. OUTLINE OF THIS PAPER 

This paper is organized as follows. Chapter 1 describes the goals of the overall 
VV&A comparison effort. Brief descriptions of the HP AC and NARAC modeling 
systems are given. Chapter 2 gives detailed descriptions of 17 baseline release scenarios 
that make up the core of the model-to-model comparison. Comparison metrics are 
introduced and illustrated in some detail. Chapter 3 reviews the results of the 
comparisons, including graphical comparisons. Excursions from the baseline model 
parameters and settings (e.g., spatial and temporal resolution) for HPAC are also 
considered. 
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In addition to the main body of this paper, there are five appendices, A through E. 
Appendix A provides an acronym list, and Appendix B presents details associated with 
the model predictions that were used for this study’s comparisons. Appendices C and D 
provide plots and graphical displays that describe the comparisons of the specific gas 
release and particle release predictions, respectively. Finally, Appendix E includes an 
extract from the pertinent task order. 
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CHAPTER 1 
INTRODUCTION 



1. INTRODUCTION 


For fiscal year 2000, The Institute for Defense Analyses (IDA) began a series of 
studies in support of the Defense Threat Reduction Agency (DTRA). The goal of these 
studies is to improve the verification, validation, and accreditation (VV&A) of hazard 
prediction and assessment models and capabilities (e.g., HPAC and NARAC). 1 These 
studies are part of a larger joint VV&A effort that DTRA and the Department of Energy 
(DOE), via Lawrence Livermore National Laboratory (LLNL), are conducting. This joint 
effort includes comparisons of LLNL and DTRA transport and dispersion (T&D) 
modeling systems, NARAC and HPAC, respectively, and their predictions. 

IDA’s role is to conduct independent analysis and special studies associated with 
this VV&A effort. This role includes conducting comparisons between the models, 
providing analysis and discussions associated with these examinations, and exploring and 
developing measures of effectiveness (MOE) that can aid hazard prediction model 
validation and accreditation. 2 

The focus of this paper is on model-to-model comparisons for a collection of 
relatively simple release scenarios. 3 Comparison metrics that have been developed in the 
course of this work are also discussed in some detail. 

A. BACKGROUND 

This VV&A effort represents a cooperative comparison in that DTRA and LLNL 
have set up a team that works closely together to provide comparable model runs and 
analysis. The two models were developed independently and against different 
requirements. The purpose of this cooperative program is to leverage resources and 
expertise for the goal of developing better modeling and simulation tools. This 
cooperative spirit is consistent with recent newsletter articles suggesting that the 


1 HPAC = Hazard Prediction and Assessment Capability and NARAC = National Atmospheric Release 
Advisory Center. 

2 Appendix E of this document contains an extract from the pertinent fiscal year 2000 task order. 

2 Seventeen release scenarios are considered. These scenarios were defined in a joint effort at the outset 
of this study. 
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challenges of hazard prediction model YV&A [Ref. 1-1] are better addressed through 
joint efforts. These joint efforts are expected to bring a certain synergy to the Chemical 
and Biological (CB) warfare modeling and simulation community [Ref. 1-2]. 

Model-to-model parametric studies can lead to the elucidation of model 
performance and quantification of relative model differences. Previous comparisons of 
the predictions of transport and dispersion model outputs demonstrated the potential for 
parametric sensitivity studies to identify, clarify, and communicate performance 
differences, including differences arising from operational assumptions used in each of 
the models [Ref. 1-3]. 

The paper focuses on model-to-model comparisons, and represents a portion of a 
larger three-year program. Recent related analyses include a careful comparison of 
HP AC and NARAC predictions of short-range, well-characterized field trial data (i.e., the 
Prairie Grass field trials). These results were completed and documented elsewhere 
[Ref. 1-4]. These initial first-year studies are expected to serve as a basis for future, more 
complex comparisons (e.g., complex terrain and weather, longer transport ranges). 

For this study, we adopt Department of Defense definitions for VV&A [Ref. 1-5]: 

. Verification - The process of determining the degree to which a model or 
simulation implementation accurately represents the developer’s conceptual 
description and specification. Verification also evaluates the extent to which 
the model or simulation has been developed using sound and established 
software engineering techniques. 4 [Department of Defense (DOD) Directive 
5000.59] 

. Validation - The process of determining the degree to which a model or 
simulation is an accurate representation of the real world from the perspective 
of the intended uses of the model or simulation. [DOD directive 5000.59] 

. Accreditation - The official certification that a model or simulation is 
acceptable for use for a specific purpose. [DOD Directives 5000.59 and 
5000.59-P], 

Model-to-model comparison studies can only address VV&A in a limited sense. 
For instance, it is quite possible that two models agree perfectly with one another, yet; 
perhaps due to deficiencies in understanding of the relevant phenomenology, disagree 
with real-world observations. However, agreement between models, when one of the 
models has already been validated against field trial data, can improve the case for 


4 This effort does not examine software engineering technique issues. 
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validation. In addition, even without specific field trial data, agreement between two 
models can boost confidence in both models’ ability to adequately predict observations. 
This is particularly so in the case of HP AC and NARAC, whose transport and dispersion 
methodologies (described below) differ distinctly. Further, identification of model 
discrepancies serves a valuable purpose. Such identifications can clarify limitations of a 
model, direct more appropriate parameter settings for a particular class of releases, or 
point to the need for new modeling techniques and/or new phenomenological 
investigations. 

B. BRIEF HPAC DESCRIPTION 

HP AC is composed of a suite of software modules that generate source terms for 
hazardous releases, retrieve and prepare meteorological information for use in a 
prediction, model the T&D of the hazardous release over time, and plot and report the 
results of these calculations. HPAC has been applied to various national defense 
problems including military studies and operational planning. 

For hazardous material T&D, HPAC uses the Second-Order Closure Integrated 
Puff (SCIPUFF) model and an associated wind field model. SCIPUFF is a Lagrangian 
model for atmospheric dispersion that uses the Gaussian puff numerical method - an 
arbitrary time-dependent concentration field is represented by a superposition of three- 
dimensional Gaussian distributions. The downwind concentration is calculated from a 
turbulent diffusion parameterization based on second-order closure theory. This 
methodology provides a link between measurable wind-flow field velocity statistics and 
predicted dispersion rates. The “second-order” feature allows concentration variance to 
be estimated (in addition to mean concentration), and this uncertainty estimate can be 
used as the basis for a probabilistic description of the dispersion prediction. 5 

This model-to-model study uses HPAC software version 3.2 [Ref. 1-7]. 

C. BRIEF NARAC DESCRIPTION 

The ADAPT/LODI 6 modeling system, within NARAC, is used for both real-time 
operational applications and detailed assessments of events involving atmospheric 
releases of hazardous material. 


5 See Reference 1-6 for details of HPAC design, functionality, capabilities, and V&V. 

6 ADAPT = Atmospheric Data Assimilation and Parameterization Techniques. LODI = Lagrangian 
Operational Dispersion Integrator. 
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The ADAPT meteorological data assimilation model constructs fields of such 
variables as the mean winds, pressure, precipitation, temperature, and turbulence, using a 
variety of interpolation methods and atmospheric parameterizations [Ref. 1-8]. Non- 
divergent wind fields are produced by an adjustment procedure based on the variational 
principle and a finite-element discretization. 

The LODI dispersion model solves the 3-D advection-diffusion equation by 
integrating a stochastic differential equation (SDE) for many (typically hundreds of 
thousands) particle trajectories in a Lagrangian Monte Carlo method [Ref. 1-9]. LODI 
includes methods for simulating the processes of mean wind advection, turbulent 
diffusion, radioactive decay and production, first-order chemical reactions, wet 
deposition, gravitational settling, dry deposition, and buoyant/momentum plume rise. 

The ADAPT/LODI models are coupled to NARAC databases providing 
topography, geographical data, chemical-biological-nuclear agent properties and health 
effects, real-time meteorological observational data, and global and mesoscale forecast 
model predictions. Graphical output is typically constructed by first smoothing raw 
LODI output with the NARAC system tool, smoothBinDat, and then plotting. 

This model-to-model study used version 2.9 of ADAPT and version 9k of LODI 
within the NARAC software. 

D. RECAP 

This work is part of a larger VV&A effort comparing two widely used transport 
and dispersion model systems, HP AC and NARAC. The focus of this work is model-to- 
model comparisons for a set of relatively simple releases (discussed in detail in the 
following chapter). The modeling architectures of HPAC and NARAC are quite distinct. 
HP AC treats the concentration field as a continuum (represented as a sum of Gaussian 
“puffs”) and propagates these puffs in time, while NARAC treats the concentration field 
as a collection of particles and propagates each particle in time. It is shown below that 
HPAC and NARAC predictions can be brought into general agreement, but only after full 
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consideration of wind profile parameterizations, spatial and temporal resolution issues, 
and other internal model-dependent parameters. 7 


7 Comparing the outputs of two codes is not necessarily a straightforward procedure. Ensuring that both 
codes’ methodologies and performance are fairly applied and assessed across any set of input 
conditions is typically a challenge. The seemingly straightforward task of comparing model outputs is 
fraught with traps that the analyst must guard against. 


1-5 



REFERENCES 


1-1. Merkle, P. B., “Analysis and Validation: A Perspective for Hazard Models,” 
Chemical and Biological Defense Information Analysis Center (CBIAC) 
Newsletter, Spring 1998, Vol. 4, No. 2. 

1-2. Gibbs, R. L., “Improving the Development Process for Chemical and Biological 
Warfare Modeling and Simulation,” CBIAC Newsletter, Fall 1999, Vol. 5, No. 4. 

1-3. Warner, S., Carpenter, J. N., Cook, J. M., Miller, R. S., and Hegemann, B. E., 
NBC Hazard Prediction Model Capability Analysis, IDA Document D-2245, 
September 1999. 

1-4. Warner S., Platt, N., Heagy, J. F., Bradley, S., Bieberbach, G., Sugiyama, G., 
Nasstrom, J. S., Foster, K. T, and Larson, D., User-Oriented Measures of 
Effectiveness for the Evaluation of Transport and Dispersion Models, IDA Paper 
P-3554, 8 January 2001. 

1-5. Director, Operational Test and Evaluation and Director, Test, Systems 
Engineering and Evaluation OUSD(A&T), Simulation, Test and Evaluation 
Process: STEP Guidelines, 4 December 1997. 

1-6. Bradley, S., Mazzola, T., Ross, R., Srinivasa, D., Fry, R., and Bacon, D., 
Verification and Validation of HP AC 3.0, for Defense Special Weapons Agency, 
June 1998, and references therein; Sykes, R. I., “HPAC/SCIPUFF: Kamisiyah 
Modeling Issues,” 3rd Annual GMU/DTRA Transport and Dispersion Modeling 
Workshop, Fairfax, VA, 28-29 July 1999; and Nappo, C. J., Eckman, R. M., 
Shankar Rao, K., Herwehe, J. A., and Gunter, L., Second Order Closure 
Integrated Puff (SCIPUFF) Model Verification and Evaluation Study, Air 
Resources Laboratory, NOAA, May 1998. 

1-7. DTRA, The HP AC User’s Guide: Version 3.2, October 1999. 

1-8. Sugiyama, G. and Chan, S. T., “A New Meteorological Data Assimilation Model 
for Real-Time Emergency Response,” 10th Joint Conference on the Applications 
of Air Pollution Meteorology, Phoenix, AZ (11-16 January 1998), Am. Met. Soc., 
Boston, MA, 285-289. 

1-9. Ermak, D. L. and Nasstrom, J. S., “A Lagrangian Stochastic Diffusion Method for 
Inhomogeneous Turbulence,” Atmospheric Environment, 2000 and Nasstrom, J. 
S., Sugiyama, G., Leone, J. M. Jr., and Ermak, D. L., “A Real-Time Atmospheric 
Dispersion Modeling System,” American Meteorological Society’s 11th Joint 
Conference on the Applications of Air Pollution Meteorology, Long Beach, CA, 
9-14 January 2000. 


1-R-l 



CHAPTER 2 

COMPARED RELEASES 



2. COMPARED RELEASES 


In this chapter, the release scenarios and the underlying rationale for their 
selection are discussed. HPAC and NARAC model specifics, including wind profile 
parameterization and grid selection, are also discussed. Several comparison measures are 
also introduced and illustrated. 

A. OVERVIEW 

In January 2000, members of the joint collaboration met at LLNL to establish a 
set of release scenarios for this model-to-model study. The joint collaboration agreed to 
seventeen releases based upon several selection criteria. 

First and foremost was the desire for simplicity so that the scenarios would test 
only the core transport and dispersion algorithms within each model. Consequently, all 
releases were chosen to have idealized weather - a defined single unidirectional vertical 
wind speed profile. Terrain was similarly elementary - flat terrain with a constant 
surface roughness of 0.008 m. Downwind transport distances of 40 km were used for all 
but one “high altitude” release. Initial source geometries were also kept simple - 
spherical sources (Gaussian with standard deviation = 1 m in all directions) for 16 of the 
17 releases and a 500 m line source in the remaining release. 

A mix of stable, near-neutral, and unstable atmospheric stability conditions was 
chosen through suitable choices of the atmospheric boundary layer height, z l7 and Monin- 
Obukhov length, L. Stable, “neutral,” and unstable conditions were chosen to have Zj/L = 
4, z/L = -l 1 , Zj/L = -10, respectively. A range of release heights (relative to boundary 
layer height) was also chosen. 

Two agent types were considered: “neutrally-buoyant gas” (9 releases) and 
“particles” (8 releases). Particle runs were carried out for three particle size distributions, 
two using fixed diameter particles of 5 microns (jl) and 50 JLL, and the third with a log¬ 
normal distribution with a mass median diameter (MMD) of 500 |l and a geometric 


1 Theoretically, neutral stability corresponds to z/L = 0, but this case was not computationally accessible 
by HPAC. 
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standard deviation (a) of 2. Releases were chosen to be either instantaneous or 
continuous, and in all cases the mass of released agent was 1 kg. Continuous sources 
assumed a constant dissemination rate for 15 minutes for a total 1 kg release. Table 2-1 
summarizes the conditions for the 17 runs that were examined. 


Table 2-1. Input Conditions for the Model-to Model Comparisons 


Run 

Source 

Source 

Geometry 

Release 

Height 

(m) 

Boundary 
Layer 
Height (m) 

L c 

(m) 

Stability d 

Duration 

of 

Release 

1 

NB gas 3 

sphere 6 

2 

100 

25 

stable 

instant e 

2 

NB gas 

sphere 

50 

100 

25 

stable 

instant 

3 

NB gas 

sphere 

80 

100 

25 

stable 

instant 

4 

NB gas 

sphere 

2 

500 

-500 

“neutral” 

instant 

5 

NB gas 

sphere 

250 

500 

-500 

“neutral” 

instant 

6 

NB gas 

sphere 

400 

500 

-500 

“neutral” 

instant 

7 

NB gas 

sphere 

2 

500 

-50 

unstable 

instant 

8 

NB gas 

sphere 

250 

500 

-50 

unstable 

instant 

9 

NB gas 

sphere 

400 

500 

-50 

unstable 

instant 

10 

5 p particles 

sphere 

2 

500 

-500 

“neutral” 

continuous 
(15 min) 

11 

5 p particles 

sphere 

250 

500 

-500 

“neutral” 

continuous 
(15 min) 

12 

5 p particles 

sphere 

400 

500 

-500 

“neutral” 

continuous 
(15 min) 

13 

50 p particles 

sphere 

2 

500 

-500 

“neutral” 

instant 

14 

50 p particles 

sphere 

2 

500 

-500 

“neutral” 

continuous 
(15 min) 

15 

50 p particles 

sphere 

750 

500 

-500 

“neutral” 

instant 

16 

Log normal f 

sphere 

2 

500 

-500 

“neutral” 

instant 

17 

Log normal 

Line 

sources 

2 

500 

-500 

“neutral” 

instant 


a NB = neutrally-buoyant. The gas considered was sulfur hexafluoride (SF 6 ), albeit with some of the 
actual material characteristics suppressed. 

b The initial sphere contains 1 kg of material and has a Gaussian density distribution with o x = G y = o z = 
1 m. 

c L = Monin-Obukhov length. 

d Zi = boundary layer height. For stable runs Zi/L = 4, for near-neutral runs Zi/L = -1, and for unstable 
runs Zi/L = -10. 

e Instant = instantaneous. 

f For the log-normal particle distribution, MMD = 500 [i and the geometric standard deviation (a) 2. 

S The initial line source contains 1 kg of material and is 500 m long. 
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B. MODEL SPECIFICS 


In this section, we review the parameter settings and input values used by the 
HP AC and NARAC models for the predictions examined in this study. Specific attention 
is given to spatial and temporal resolution, meteorological settings, and vertical wind 
profiles. 

1. Meteorology: Vertical Wind Profiles 

Preliminary model-versus-model (MvM) runs revealed differences in output 
concentrations that could be attributed to differences in vertical wind profile 
parameterizations between the two models. For this reason, it was agreed to specify the 
user-controlled NARAC vertical wind profile for each run, according to the documented 
HP AC vertical wind profile parameterization, given by [Ref. 2-1]. 



In (z/z 0 )-yr m (z,L) 
In {zJz 0 )-V n (z„L) 


z<z s 


\z>z s 


( 2 - 1 ) 


with the following identifications: 

. u : wind speed (m/s) 

. z: height (m) 

• zo- surface roughness (m); 0.008 m for all runs 
. Zs- surface layer height (m) 

• L: Monin-Obukhov length (m) 

. u s : wind speed (m/s) at the surface layer height Zs 
. yf m : stability correction. 

The functional form of the stability correction and the value of the surface layer height 
depend on the Monin-Obukhov length used; details are provided in Reference 2-1. The 
seventeen MvM runs used three distinct vertical wind profiles: 

. Profile 1: Stable runs (1 through 3) 

• Profile 2: Near-neutral runs (4 through 6, 10 through 17) 

. Profile 3: Unstable runs (7 through 9). 
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For all runs, the surface layer height is z s = 50 m, above which the wind speed is 
constant. This consistent value of 50 m was essentially a coincidence, resulting from the 
particular parameter set chosen. The three profiles are illustrated in Figure 2-1 below. 
Note that all three wind profiles pass through the observation point, u = 5 m/s at z = 10 m, 
as required. 


Vertical Wind Profiles for Model-to-Model Runs 



Figure 2-1. Vertical Wind Profiles for Model-Versus-Model Runs 
2. HPAC Computational Details 
a. HPAC Parameters and Settings 

HPAC has several user-controlled parameters that define specifics of the release, 
meteorological conditions, and spatial and temporal resolutions in internal computations 
and external output. For the model-to-model comparisons conducted in this study, it was 
agreed to define baseline parameter settings. Four key HPAC parameters in the baseline 
set are vertical resolution, horizontal resolution, maximum time step, and conditional 
averaging time, discussed separately below. Other baseline HPAC parameters are 
discussed in Appendix B. 
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i. Vertical Resolution 


Vertical resolution for HPAC runs is set within the Domain Editor in the Project 
Editor. For all MvM runs, with the exception of MvM 15, the baseline vertical resolution 
was set to 10 m. This value was found to be a reasonable compromise between adequate 
sampling of the vertical wind profile and the extent of the vertical computational 
domain. 2 The vertical resolution setting limits the maximum height of the vertical 
domain (also set in the Domain Editor ) to be less than or equal to 50 times the resolution 
value. The 10 m value therefore restricts the domain height to be no greater than 500 m. 
A vertical domain height of 500 m was selected for all runs except MvM 15 (note that 
maximum release height for these runs is 400 m). For MvM 15, vertical resolution was 
set to 20 m and the vertical domain height was set to 1,000 m. 

ii. Horizontal Resolution 

Horizontal resolution for HPAC runs is set within the Domain Editor in the 
Project Editor. For all MvM runs, the horizontal resolution was set to the default value. 
Fine horizontal resolution values (e.g., 200 m) were found to give rise to high frequency 
spatial oscillations in the HPAC dosage plumes that may be connected with SCIPUFF 
puff-splitting algorithms used in HPAC (discussed in more detail in Section C of Chapter 
3). Note that horizontal resolution is an internal parameter, in that it affects only the 
internal representation of the plume and its evolution. The horizontal resolution does not 
affect the spatial locations of HPAC output values; these are set via the sampler locations 
within the sampler files. 

iii. Maximum Time Step 

The Maximum Time step for HPAC runs is set within the Time Editor in the 
Project Editor. For all runs, the value was set to 60 seconds. This parameter determines 
the time step used in the output sampler (.smp) file. Excursions from the 60-second 
baseline value are taken up in Section C of Chapter 3. 

iv. Conditional Averaging Time 

The conditional averaging time for HPAC runs is set within the Options Editor in 
the Project Editor. Conditional averaging is used by HPAC to incorporate large-scale 
uncertainties in meteorological conditions (turbulent velocity fluctuations) into output 


2 Section C in Chapter 3 addresses excursions from the baseline vertical resolution. 
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concentration estimates [Ref. 2-2], The conditional averaging time controls the extent to 
which the large-scale (meandering) components are included in the concentration 
estimates. Larger conditional averaging times typically produce larger plumes. For all 
runs, the conditional averaging time was set to 1 hour. 

b. HPAC Surface Observation Files 

In order to specify the meteorological quantities listed in Table 2-1, it was 
necessary to create HPAC surface observation files [Ref. 2-3] for the three stability 
classes considered. 3 Each file specifies a single wind observation of 5 m/s in the negative 
x direction (90 degrees) at a height of 10 m, geo-location of 45° N latitude and 0° E 
longitude, and local time of 0:00 hours. 4 Also specified are the boundary layer heights 
(Zi) and Monin-Obukhov lengths (L) for the runs. The three surface observation files are 
displayed in Figures 2-2 through 2-4. 


# TYPE: OBSERVATION 

# TIMEREFERENCE: UTC 

# Written by: Weather File Editor Version 3.2.3 

# Written on: February 21,2000 @ 13:53 

# 

SURFACE 

11 

ID YYMMDD HOUR LAT LON ELEV ZI Z WDIR 
WSPD MOL 
# 

HOURS N E M M M DEG M/S M 
-9999 

mvmS 000401 0.00 45. 0. 10 100 10 90 5.0 25 


Figure 2-2. Surface Observation File for Stable Releases (MvM 1 through 3) 


3 For example, the boundary layer height and the Monin-Obukhov length cannot be directly specified 
via the HPAC user interface. 

4 45° N latitude and 0° E longitude is the release point for all HPAC releases. The time for all releases is 
0 hours local time. Both geo-location and time of the release are arbitrary choices and do not influence 
the HPAC results. This would not be the case, however, if real terrain was added (depends on geo¬ 
location) or agent decay effects were included (depends on time of day). 


2-6 



# TYPE: OBSERVATION 

# TIMEREFERENCE: UTC 

# Written by: Weather File Editor Version 3.2.3 

# Written on: February 21,2000 @ 13:53 

# 

SURFACE 

11 

ID YYMMDD HOUR LAT LON ELEV Zl Z WDIR 
WSPD MOL 
# 

HOURS N E M M M DEG M/S M 
-9999 

mvmN 000401 0.00 45. 0. 10 100 10 90 5.0 -500 


Figure 2-3. Surface Observation File for Neutral Releases (MvM 4-6,10-17) 


# TYPE: OBSERVATION 

# TIMEREFERENCE: UTC 

# Written by: Weather File Editor Version 3.2.3 

# Written on: February 21,2000 @ 13:53 

# 

SURFACE 

11 


ID YYMMDD HOUR LAT 

LON 

ELEV 

Zl Z WDIR 

WSPD 

# 

HOURS 

MOL 





N E M M 

M DEG 

M/S 

M 

-9999 






mvmU 000401 0.00 45. 0. 

10 100 

10 

90 

5.0 -50 


Figure 2-4. Surface Observation File for Unstable Releases (MvM 7-9) 


c. HPAC Output and Dosage Calculation 

At the outset of this study, it was agreed that the model-to-model comparison 
would focus on comparisons of dosage, sometimes referred to as (time) integrated 
concentration. HPAC output is held in a sampler output file (.smp file) [Ref. 2-4] that 
contains the concentration c(iv k ) at each sampler location r; and time step t k . For the 
purposes of this study, the dosage at a sampler point r; is defined to be the time integral 
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(sum) of the concentration at that point from the time of the release, t = 0, to the current 
time t, 

N 

d(r i ,t) = At^c(r i ,t k ), (2-2) 

k =0 

where At is the output time step (typically 60 sec), t = N At is the current time, and N is 
the number of time steps (in MKS units, dosage has units of kg-s/m 3 ). A routine to read 
HP AC sampler output files and compute dosages for all sampler locations was written in 
MATLAB. Dosages for all runs were computed at t = 30 min and t = 60 min. Dosages 
for MvM run 15 were computed at the additional times t = 120, 180, and 240 min, 
because of the longer run duration. 

3. NARAC Computational Details 5 

a. NARAC Meteorological Grids, Parameters, and Settings 

NARAC meteorological fields including mean winds and turbulence were 
generated by ADAPT 6 on a user-specified meteorological data grid generated by the 
Gridgen utility with flat topography specified. For MvM runs 1 through 14, 16, and 17 
the meteorological data grid was constructed with 0.5 km horizontal resolution covering 
an area 40 km in the downwind direction (x) and 8 km in the transverse (y) direction. For 
the high altitude run, MvM 15, the meteorological grid was again constructed with 0.5 
km resolution, but covered an area 100 km by 20 km in x and y, respectively. The 
horizontal resolution is unimportant, since there is no horizontal variation of the 
meteorology in these simulations. For all runs, a 26-level graded vertical grid was used 
with levels at 0.0, 10.0, 21.6, 35.1, 50.7, 68.8, 89.9, 114.2, 142.6, 175.3, 213.4, 257.6, 
308.9, 368.4, 437.4, 517.4, 610.3, 718.1, 843.2, 988.2, 1156.5, 1351.8, 1578.4, 1841.2, 
2156.2, and 2500.0 m. 

The HP AC vertical wind profile was input to ADAPT via a NARAC standard 
format observational file. The main input file to ADAPT is a namelist specifying 
information about the grid, the meteorological data, and the choice of assimilation 
methods used to generate the 3D meteorological fields used by LODI. The critical 
parameters for this study are the surface roughness, surface layer depth, boundary layer 


5 See Reference 2-5 for additional details associated with the NARAC modeling system. 

6 See Chapter 1, Section C and Reference 2-6. 
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depth, and the inverse Monin-Obukhov length. ADAPT also generates eddy diffusivities 
based on similarity theory parameterizations. 

A 2-meter cutoff distance was selected for all runs using a spherical Gaussian 
source. This is the distance at which the Gaussian distribution is cut off, yielding a finite 
area over which to initialize source particles. (Several runs done with a cutoff of 5 m 
displayed only minor differences in the results.) 

All sources were initialized with 500,000 particles. 7 This allowed convergence in 
the concentration field to roughly three significant figures along the downwind axis. 

Name lists (input files) for all of the NARAC predictions can be found in 
Appendix B. 

b. NARAC Dosage Calculation 

Integrated air concentration (i.e., dosage) is defined by: 

t 

IC(x,y,z,t)= jc(x,y,z,t')dt' (2-3) 

t-At 

where t is the output time and At is the sampling period. NARAC integrated air 
concentrations in the xy plane at z = 10 m 8 were computed every 30 and 60 minutes, 
using a 30- or 60-minute sampling period, respectively. (Note that MvM 15 also 
produced integrated concentrations using 120-, 180-, and 240-minute sampling periods.) 
Air concentrations were also calculated in the xz plane at y = 0 m at particular times 
every 30 minutes. Runs with particle sources also computed instantaneous surface 
deposition every 30 minutes. 

4. NARAC and HPAC Output Comparison Grids 

For each MvM run, HPAC and NARAC output data (typically concentration data) 
were generated on spatial grids. Grids were chosen to respect the spatial extent of the 
plume for the duration of the release. Two grids were adopted: 

. Grid 1: used for all MvM runs except MvM run 15 (high altitude release) 

. Grid 2: used for MvM run 15. 


7 The 500,000 particles used in the NARAC calculation do not correspond to individual molecules. 
Rather, they represent computational constructs that have the appropriate properties for the material 
being modeled. 

8 The concentration, C , is determined using a sampling depth of 20 meters. See Section 4 below. 
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Specifics on the construction of the grids and their conversion to sampler files follow. 


a. NARAC Grids 

The NARAC suite of modeling tools includes the utility Gridgen, which is used to 
generate variable-resolution concentration grids. The grids are constructed to resolve the 
spatial distribution of plume concentration as a function of downwind distance from the 
release point. This is achieved by satisfying the relationship Ax = <J x ln for ri> 2, where Ax 
is the grid spacing (or resolution) and cr x (t) is the standard deviation of the plume 
concentration distribution and n determines the desired resolution. Using this 
relationship and the approximations that <7, grows linearly with time, 

^M = o x (0) + aj, (2-4) 

we can define a grading factor, f x = Ax/x such that 

/,=?*-. (2-5) 

un 

where u is the mean wind speed, x is the downwind distance, t-xlu, and <7„ is the 

standard deviation of the x component of the wind velocity. This assumes that <7 X ( 0) is 
much less than a u t (which is valid for small source sizes or for sufficiently large t). 

Typical one-hour-average values of <7 U are 0.5 to 1.0 m/s. For the runs we have 
defined, u = 5 m/s, so using n = 2 and Equation 2-5, f x = (0.5 m/s)/[(5 m/s)(2)] = 0.05. 

For example, at a downwind distance corresponding to one hour after the release, this 
grading factor leads to a grid spacing of Ax =f x ut = (0.05)(5 m/s)(3600 s) = 900 m. The 

same grading factor was used for both the downwind (x) and the crosswind (y) directions. 

Grid 1 has a downwind (x) extent of ~ 40 km and transverse (y) extent of ~ 8 km. 
The grid was constructed with 50 m resolution out to a distance of 1 km in both the 
downwind and transverse dimensions. This grid was then graded outward with a grading 
factor of 0.05, yielding a mesh with 130 (96) points in the x (y) direction, for a total of 
12,480 grid points. The maximum cell size in the x direction is 1,435 m and in the y 
direction it is 189 m. This resolution is sufficient to resolve the plume at all distances 
greater than 1 km in both the downwind and transverse dimensions. 

MvM run 15 requires a larger concentration grid, due to the height of the release 
and the longer duration. The grid for this run, grid 2, has a downwind (x) extent of ~ 90 
km and a transverse (y) extent of ~ 18 km. The grid uses the same base grid and grading 
factor as defined above, but the resulting mesh is 116 points by 136 points in the x and y 
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directions, respectively, giving a total of 15,776 grid points. The maximum cell size in 
the x direction is 4,607 m and in the y direction it is 455 m. 

Because LODI is a Lagrangian particle code, a sampling volume must be defined 
in order to calculate concentrations. This volume is defined by specifying a sampling 
depth perpendicular to the two-dimensional bin already defined. We chose a sampling 
depth of 20 meters implying that the volume used to calculate the air concentrations at 
10-meter height extends from z = 0 to z = 20 m. 

The NARAC smoothBinData utility was used to smooth the concentration field 
prior to plotting. This utility uses an area-weighted filter to smooth the data on the 
graded concentration grid. For concentration grids with constant spacing the filter 
reduces to a 1-2-1 filter along each dimension [Ref. 2-7]. Use of smoothing is desirable 
due to the rapid fall-off in particle number as the edge of the plume is approached. 

b. HPAC Grids 

All HPAC results in this comparison were generated using HPAC version 3.2 
operated in the Extended Mode. In order to compute HPAC concentration values on the 
NARAC-defined grid points, it was first necessary to convert grids 1 and 2 to HPAC 
sampler (.sam ) files [Ref. 2-8]. HPAC sampler files are limited to 10,000 grid points. 
This limitation required the NARAC grids to be suitably reduced before conversion to 
sampler files. Point reduction was done by eliminating every other row of downwind 
points, starting near the edges in the crosswind direction and working inward, 
symmetrically about the centerline. This process was continued until the remaining 
number of grid points dropped below 10,000. After reduction, the grid sizes were: 

. Grid 1: (130 downwind) x (76 crosswind) for a total of 9,880 points 

. Grid 2: (116 downwind) x (86 crosswind) for a total of 9,976 points. 

Figures 2-5 and 2-6 show the two grids, before and after reducing the number of points. 
Throughout this study, output was compared only on those points common to the HPAC 
and NARAC grids, that is, only on the reduced grids. Unless otherwise mentioned, the 
sampler height was chosen to be 10 m for all samplers for all runs. This value coincides 
with the height used in the NARAC runs for integrated concentrations. 
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Grid 1 for MvM Runs 1-14,16,17 



Figure 2-5. Grid Points for Grid 1: 

Blue Points Denote Those Retained for HPAC Sampler File 


Grid 2 for MvM Run 15 
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Figure 2-6. Grid Points for Grid 2: 

Blue Points Denote Those Retained for HPAC Sampler File 
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C. COMPARISON MEASURES 


In this section, we discuss the measures used to compare the HPAC and NARAC 
output dosages. In selecting comparison measures, we decided to adopt measures that 
answer basic questions about the dosage plumes: 

. How far downwind is the plume? 

. How wide is the plume? 

. What is the maximum dosage at a given downwind distance? 

. To what extent do the two model plumes overlap? 

For this purpose, four quantitative measures were adopted: downwind contour 
distances for selected dosages, crosswind dosage plume widths, crosswind maximum 
dosages, and an area-based measure of effectiveness (MOE). These measures are easily 
computed and communicate in a quantitative fashion much of the information one learns 
from simple graphical comparisons of output dosages. For the high altitude release 
(MvM 15), vertical plume structure was also investigated. Details are provided in 
Chapter 3. 

1. Dosage Contours 

Dosage contours provide immediate graphical information about the relative 
performance of model-to-model results; model agreement or lack thereof can be quickly 
assessed. Within a single run, the dosage can span a large range of values (for the 
scenarios that we examined, typically 10' 12 -10' 6 kg-s/m 3 ). Because of the wide dynamic 
range, contours are usually spaced at powers of 10 (or the logio of the dosage is viewed). 
To avoid processing extremely small or zero dosage values, all dosage data for both 
models were truncated at a minimum threshold value of 10' kg-s/m . A MATLAB 
routine was written to process the dosage data for both models and construct the contour 
plots. Contours were computed using the built-in MATLAB contouring algorithm, 
contour. Contour values were selected by first viewing the histogram of the log of the 
dosage values for the run, then choosing upper and lower contours at powers of 10 near 
the tails of the histogram, and filling in with intermediate powers of 10. This technique 
takes much of the guesswork out of contour value selection. 
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2. Downwind Contour Distances 

Dosages along the centerline (x axis) for both models are computed by 
interpolating 9 two-dimensional dosage data onto the centerline. For a given dosage value 
the downwind distance to that dosage is computed for both models. 

3. Crosswind Dosage Plume Widths 

There are many ways to define the plume width. The width used in this 
comparison is a dose-weighted width, defined as follows. Let the dosage values at a 
given crosswind slice be denoted by di, i = 1,2, ..., n, with corresponding crosswind 
positions y{. Then the dose-weighted mean location is given by 

( 2 - 6 ) 
i =1 

where the weights a, are given by 

a i =-^~- (2-7) 

7=1 

Similarly, the dose-weighted mean square position is 

(>’ 2 ) = S«, yf • (2-8) 

i =1 

The dose-weighted width is then taken to be 

W = 2^{y 2 )-{y) 2 . (2-9) 

Dosage plumes are frequently approximated by a Gaussian in the crosswind 
direction. 10 If the points (yi, di) describe a Gaussian, then the mean position (2-6) and 
width (2-9) correspond to the mean, ju, and twice the standard deviation, 2<j, respectively 
for the Gaussian. Crosswind dosage plume widths are compared for several downwind 
locations for each run. 


9 Interpolation was carried out with the MATLAB routine interpl. 

10 This is no surprise in the case of HP AC, since the HP AC plume is represented as a superposition of 
Gaussian distributions. 
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4. Crosswind Maximum Dosages 

Using the same crosswind dosage values as in the plume width calculation above, 
the maximum dosage, c/ max = in ax( <:(.), is found for both models at each downwind 

i 

location. 

5. Area-Based MOE 

The area-based MOE was introduced in Reference 2-9, where it was used to 
compare model predictions against field trial data. With a minor change in interpretation, 
this MOE can be applied to model-to-model only (i.e., no field trial data) comparisons. A 
description of the MOE as it applies to model comparisons follows. 

In any comparison of two model predictions, there are three regions of interest: 

. Model 1 and Model 2 predictions agree (e.g., overlap for some dosage 
contour) 

. Model 1 overpredicts Model 2 (Model 1 predicts a larger area at the specified 
contour level than Model 2) 

• Model 1 underpredicts Model 2 (Model 1 predicts a smaller area at the 
specified contour level than Model 2). 

Figure 2-7 illustrates dosage contours, Mi and M 2 , for two models at the same 
contour level. The green (solid) overlap area, A OL , represents the region where the model 
predictions agree. The blue (cross-hatched) region, A 0 p, represents the region where 
model 1 overpredicts model 2. The red (diagonal lined) region, Aup, represents the 
region where model 1 underpredicts model 2. 



Figure 2-7. Illustration of Three Regions for Two Model Predictions M! and M 2 
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The area-based MOE is a two-dimensional quantity, defined in terms of the three 
areas, Aql, Aup, A 0 p, as follows: 


MOE = 


1- 


A, 


VP 


A up A ( ,l 


1 --—- 1 . 

Aqp Aql ) 


( 2 - 10 ) 


Noting that the areas inside contours Mi and M 2 can be expressed as A M i = A 0 p + A 0 l 
and A M 2 = Aup + A 0 l, respectively, the MOE can be equivalently expressed as 

( 


MOE = 


i_jV 5 i_. 


A 


■M2 


A*.] 


' A ol Aol ' 

A m) 


4 ’ 4 

V /W 2 ^M\ / 


( 2 - 11 ) 


Geometrically, the MOE occupies a point in the two-dimensional region shown in 
Figure 2-8. If the model predictions agree perfectly, then Ami = A M 2 = A OL , and the MOE 
takes on the value (1,1). Alternatively, if the model predictions disagree completely, then 
there is no overlap region; in this case A 0 l = 0 and the MOE is (0,0). Models that are in 
close agreement will have MOE values concentrated in the vicinity of the (1,1) point. 
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A model that consistently overpredicts or underpredicts the other can be easily 
diagnosed from MOE behavior. Consider the diagonal line in Figure 2-8. On this line, 
the areas A M i and A M 2 are the same. This in turn implies that the areas Aup and A 0 p are 
the same. In other words, underprediction and overprediction are equally likely. Below 
the diagonal, A M i > A M 2 or A 0 p > Aup, represents the region of consistent overprediction 
of model 1 when compared to model 2. Above the diagonal, A M i < A M 2 or A 0 p < Aup, 
corresponds to consistent underprediction of model 1 relative to model 2. 

In practice, the contours for HP AC and NARAC, generated using the MATLAB 
routine contour, are polygons. Contour areas are computed using the MATLAB routine 
polyarea. Without loss of generality, we choose model 1 to be HP AC and model 2 to be 
NARAC throughout this work; therefore A M i = A H pac and A M 2 = A NA rac- MOEs for all 
runs and several contour levels are computed once the basic areas A OL , A H pac, and 
Anarac are found. 
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RESULTS OF COMPARISONS 



3. 


RESULTS OF COMPARISONS 


In this chapter, results of the model-to-model comparison are presented. Gas 
releases (runs 1 through 9) are discussed first, followed by the particle releases (runs 10 
through 14, 16, and 17). The high altitude particle release, run 15, is considered 
separately. This chapter concludes with a discussion of excursions from the baseline 
HP AC resolution settings. 

A. NEUTRALLY-BUOYANT GAS RELEASES 

Nine neutrally-buoyant gas releases were compared: three stable cases (runs 1 
through 3), three near-neutral cases (runs 4 through 6), and three unstable cases (runs 7 
through 9). Overall agreement between the HP AC and NARAC gas release predictions 
was good, and in some cases excellent, with the exception of run 3, which had the highest 
release height of all the stable releases. In general, agreement for the stable releases was 
less favorable than for the near-neutral and unstable releases. Also, there was generally 
better agreement for low altitude releases than for higher altitude releases. Details are 
presented below. 

1. Representative Contour Plots 

Comparative dosage contour plots form a mainstay of this comparison effort, and 
examples are discussed below. The full collection of contour plots can be found in 
Appendix C. Contour plots were constructed for all gas runs at 30 minutes and 60 
minutes after the release. Contour values were chosen by first viewing histograms of 
HP AC and NARAC dosages (after taking the logarithm), then selecting contours that 
cover the bulk of the overlap region of the two histograms. Figure 3-1 shows the dosage 
histogram for run 1 (stable atmospheric conditions) at 60 minutes. 1 The lower and upper 
contours are selected at 10~ 10 kg-s/m 3 and 10~ 5 kg-s/m 3 , respectively, with intermediate 
contours spaced in powers of 10. Figure 3-2 shows the corresponding contour plot; the 
lower portion of the figure shows the dosage along the downwind symmetry axis. The 


1 This histogram is based on all dosages in the output domain (that is, dosages at each output grid point 
are considered). Throughout this comparison solid red lines denote NARAC quantities, while dashed 
blue lines denote HP AC quantities. 
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release location for all runs is at (x,y) = (0,0) and the downwind direction is the negative 
x direction. Unless otherwise stated, dosages are computed 10 m above the surface. 


MvMOl: Dosage Histogram Comparison, Time = 60 min. 



Figure 3-1. MvM 1 (Stable) at 60 Minutes: Histogram of NARAC (Red -) and HPAC (Blue - -) 
Dosages with Lower and Upper Contours at 10' 10 kg-s /m 3 and 10' 5 kg-s /m 3 , Respectively 


MvMOl: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure 3-2. MvM 1 (Stable) at 60 Minutes with Upper: NARAC (Red -) and HPAC (Blue - -) 
Dosage Contours; Lower: Log Dosages on Downwind Symmetry Axis 
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Figures 3-3 and 3-4 show similar dosage contour plots for runs 4 (near-neutral) 
and 7 (unstable) at 60 minutes, respectively. 2 These three figures point to a general trend 
in the results of this comparison: model agreement improves as atmospheric conditions 
become increasingly unstable . 3 This and other trends are substantiated in the following 
sections. 


MvM04: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure 3-3. MvM 4 (Near-Neutral) at 60 Minutes with Upper: NARAC (Red -) and HPAC 
(Blue - -) Dosage Contours; Lower: Log Dosages on Downwind Symmetry Axis 


2 The HPAC dosage oscillations near the origin are discussed in Section C.3 below. These oscillations 
do not affect the conclusions drawn in this Chapter. 

3 Note that the release height for runs 1, 4, and 7 is 2 m. Other than atmospheric stability, all other 
parameters are the same. 
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MvM07: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure 3-4. MvM 7 (Unstable) at 60 Minutes with Upper: NARAC (Red -) and HPAC 
(BLUE - -) Dosage Contours; Lower: Log Dosages on Downwind Symmetry Axis 

2. Scatter Plot Comparison of Downwind Contour Distances 

Transport (and downwind diffusion) of the plume between HPAC and NARAC 
was quantified by comparing the downwind distances to given dosage contours along the 
symmetry axis. The contour distances were computed by first interpolating dosages onto 
the symmetry axis using the MATLAB interpolation routine interpl. A bisection 
algorithm was then used to find the distances where the dosages crossed the given 
contour values. Figures 3-5 and 3-6 show scatter plots of HPAC versus NARAC contour 
distances at 30 minutes and 60 minutes, respectively. Colors indicate stability category: 
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blue = stable, red = near-neutral, and green = unstable. 4 The dashed lines above and 
below the diagonal correspond to fixed overprediction limits. For example, in the 30- 
minute case (Figure 3-5), the upper dashed line represents the HP AC distance exceeding 
the corresponding NARAC distance by 20 percent. Similarly, the lower dashed line 
corresponds to the NARAC distance exceeding the corresponding HP AC distance by 20 
percent. In general, the HP AC contours tend to overpredict 5 (extend farther than) 
corresponding NARAC contours, but by no more than ~ 20 percent at 30 minutes and by 
no more than ~ 15 percent at 60 minutes. Agreement between the models is seen to 
improve with decreasing stability; that is, the green = unstable points are closest to the 
diagonal, while the blue = stable dots are the farthest. 


Gas Runs: Downwind Contour Distances, 30 Minutes 



NARAC Dosage Contour Distance (km) 

Figure 3-5. Scatter Plot Comparison of Downwind Contour Distances at 30 Minutes for All 
Gas Runs: Colors Indicate Meteorological Stability Categories and Dashed Lines are 20% 

Overprediction Limits 


4 This color convention is used throughout this chapter. 

5 Equivalently, NARAC underpredicts HP AC. Throughout this study when drawing comparisons we 
adopt the nomenclature that one model overpredicts the other. This terminology is not to be construed 
as a statement about the essential correctness or incorrectness of either model. 
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Gas Runs: Downwind Contour Distances, 60 Minutes 



Figure 3-6. Scatter Plot Comparison of Downwind Contour Distances at 60 Minutes for All 
Gas Runs: Colors Indicate Meteorological Stability Categories and Dashed Lines are 15% 

Overprediction Limits 

Based on the above results, it appears that the transport (and downwind diffusion) 
components of the two modeling systems for neutrally-buoyant gases compare favorably. 

3. Scatter Plot Comparison of Crosswind Dosage Plume Widths 

Figures 3-7, 3-8, and 3-9 show crosswind plume profiles at four selected 
downwind distances for runs 1, 4, and 7, respectively; the time after release in all plots is 
60 minutes. These figures are representative of the crosswind plume profiles for all gas 
runs (with the exception of run 3, which is considered separately below). Appendix C 
contains crosswind plume profiles for all gas runs. Dose-weighted plume widths were 
computed from expression (2-9) at four downwind distances for the 30-minute and 
60-minute plumes for all gas runs. For the 30- (60-) minute plume, the downwind 
distances were typically: 1 km, 3 km, 5 km, and 10 km (5 km, 10 km, 15 km, and 20 km). 
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At shorter distances, the NARAC plume is wider than the HP AC plume. Further 
downwind, HP AC begins to, and then increasingly, overpredicts NARAC plume widths. 


1Q s MvMOl: Crosswind Dosage Profiles, Time = 60 min. 



Crosswind Distance (km) 


Figure 3-7. MvM 1 at 60 Minutes with Crosswind Plume Profiles at Selected Downwind 
Distances: NARAC (Red -) and HPAC (Blue - -) 


Figures 3-10 and 3-11 show scatter plot comparisons of the HPAC and NARAC 
plume widths at 30 minutes and 60 minutes, respectively. Agreement for the unstable 
cases is best, with the bulk of the unstable (green) points falling within the 20 percent 
limits. The stable and neutral cases show greater differences, with the bulk of the points 
(blue and red) falling outside the 20 percent limits, and a few points falling outside the 50 
percent limits. 
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Dosage (kg s/m ) 


x 10 


MvM04: Crosswind Dosage Profiles, Time = 60 min. 



Figure 3-8. MvM 4 at 60 Minutes with Crosswind Plume Profiles at Selected Downwind 
Distances: NARAC (Red -) and HPAC (Blue - -) 


10 -<s MvM07: Crosswind Dosage Profiles, Time = 60 min. 



Crosswind Distance (km) 


Figure 3-9. MvM 7 at 60 Minutes with Crosswind Plume Profiles at Selected Downwind 
Distances: NARAC (Red -) and HPAC (Blue - -) 
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The trend for NARAC plumes to be wider than HP AC plumes at the smaller 
downwind distances (1 km), but narrower at the larger downwind distances at 30 min is 
clearly seen in Figure 3-10. This observation is consistent with the results from 
comparisons of HP AC and NARAC predictions of the Prairie Grass field trials, which 
incorporated a maximum downwind range of 800 meters [Ref. 3-1]. 

Gas Runs: Crosswind Plume Widths, 30 Minutes 



NARAC Plume Width (km) 


Figure 3-10. Scatter Plot Comparison of Crosswind Plume Widths at 30 Minutes for All 
Gas Runs: Colors Indicate Meteorological Stability Categories and Dashed Lines Show 

20% and 50% Overprediction Limits 
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For the three gas release cases examined with a source release height of 2 meters 
(which was most consistent with the near-surface Prairie Grass field trials), the 
crosswind plume widths for the HP AC and NARAC predictions 1 km downwind are 
shown in Figures C-2 (MvM 1, stable), C-20 (MvM 4, neutral), and C-38 (MvM 7, 
unstable). Consistent with the comparative Prairie Grass results, NARAC’s largest 
overpredictions of HP AC (crosswind plume width) occur for the stable and neutral 
conditions [Ref. 3-2], For the unstable case, MvM 7 (Figure C-38), the 1 km crosswind 
dosage plume widths are similar (as was true for the Prairie Grass comparisons). 


Gas Runs: Crosswind Plume Widths, 60 Minutes 



Figure 3-11. Scatter Plot Comparison of Crosswind Plume Widths at 60 Minutes for All 
Gas Runs: Colors Indicate Meteorological Stability Categories and Dashed Lines Show 

20% and 50% Overprediction Limits 
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4. Scatter Plot Comparison of Crosswind Maximum Dosages 

Maximum crosswind dosages were computed at the downwind plume locations 
used in the plume width calculations above. HP AC versus NARAC dosage data for all 
gas runs at 30 minutes and 60 minutes are displayed in Figures 3-12 and 3-13, 
respectively. In each plot, the size of the markers is proportional to the downwind 
distance. Colors indicate the meteorological stability categories. Also shown on the plots 
are three HP AC overprediction limits: 10 times, 100 times, and 1,000 times the NARAC 
dosage, and 1 NARAC overprediction limit, 10 times the HPAC dosage. 

Several general trends can be gleaned from the two plots. There is reasonable 
agreement for the unstable (green) and neutral (red) cases at both times, with a slight 
tendency for HPAC to overpredict NARAC at 30 minutes and for NARAC to overpredict 
HPAC at 60 minutes. For all stability categories, there is a tendency for closer agreement 
at shorter downwind distances 

The large HPAC overpredictions for the stable (blue) cases are primarily due to 
run 3; these points are noted in the figures. Run 3 is problematic. Dosage contours and 
crosswind profiles for this run are shown in Figures 3-14 through 3-17. There is 
substantial disagreement between the predictions at both 30 and 60 minutes. The height 
of this release plays an important role in the model differences. The release height for 
run 3 is the highest of the stable runs - 80 m - versus 2 m for run 1 and 50 m for run 2. 
The boundary layer height for all three stable runs is 100 m. Since the 
80 m-release height is near the boundary layer, it is reasonable to ask to what extent the 
vertical diffusivity parameterizations and blending above and below the boundary layer 
height affect the HPAC and NARAC results. The resolution of this discrepancy, which is 
a topic of ongoing investigations, is discussed in more detail in Section B.3 of this 
chapter. 
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HPAC Dosage (kg s/m ) 


Gas Runs: Max Crosswind Dosages, 30 Minutes 



Figure 3-12. Scatter Plot Comparison of Maximum Crosswind Dosage at 30 Minutes for All 
Gas Runs: Colors Indicate Meteorological Stability Categories, 

Marker Size is Proportional to Downwind Distance, 

Dashed Lines Show lOx, lOOx, and lOOOx HPAC Overprediction Limits and lOx NARAC 
Overprediction Limit, and Indicated Outliers Are from Stable Run 3 
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HPAC Dosage (kg s/m 3 ) 


Gas Runs: Max Crosswind Dosages, 60 Minutes 



Figure 3-13. Scatter Plot Comparison of Maximum Crosswind Dosage at 60 Minutes for All 
Gas Runs: Colors Indicate Meteorological Stability Categories, 

Marker Size is Proportional to Downwind Distance, 

Dashed Lines Show lOx, lOOx, and lOOOx HPAC Overprediction Limits and lOx NARAC 
Overprediction Limit, and Indicated Outliers Are from Stable Run 3 
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MvM03: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure 3-14. MvM 3 at 30 Minutes with Upper: NARAC (Red -) and HPAC (Blue - -) Dosage 
Contours; Lower: Log Dosages on Downwind Symmetry Axis 

MvM03: Crosswind Dosage Profiles, Time = 30 min. 
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Figure 3-15. MvM 3 at 30 Minutes with Crosswind Plume Profiles at Selected Downwind 
Distances: NARAC (Red -) and HPAC (Blue - -) (Note Log Dosage on Vertical Axis) 
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MvM03: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure 3-16. MvM 3 at 60 Minutes with Upper: NARAC (Red -) and HPAC (Blue - -) Dosage 
Contours; Lower: Log Dosages on Downwind Symmetry Axis 

MvM03: Crosswind Dosage Profiles, Time = 60 min. 
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Figure 3-17. MvM 3 at 60 Minutes with Crosswind Plume Profiles at Selected Downwind 
Distances: NARAC (Red -) and HPAC (Blue - -) (Note Log Dosage on Vertical Axis) 
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5. Area-Based MOE 


The final comparison measure applied to the neutrally-buoyant gas runs is the 
area-based MOE introduced in Chapter 2. For a given contour level, the calculation of 
the area-based MOE requires three areas: the HPAC contour area, A H pac; the NARAC 
contour area, Anarac; and the area of the overlap region, Aol- In MATLAB, the dosage 
contours and overlap region are represented as polygons, whose vertices are stored in 
user-accessible arrays (see Figure 3-18). The three areas are estimated by computing the 
enclosed areas of the polygons via the MATLAB routine polyarea. 6 Here, and 
throughout, we take the x-axis of the two-dimensional MOE to be Aol/Anarac and the y- 
axis to be Aol/Ahpac- Perfect model agreement then corresponds to the point (x,y) = 
(Aol/Anarac,Aol/Ahpac) = (1,1), while “perfect” model disagreement corresponds to the 
point (x,y) = (0,0). The diagonal line from (0,0) to (1,1) divides the MOE space into two 
regions: a region of HPAC overprediction, occurring below the diagonal, and a region of 
NARAC overprediction, occurring above the diagonal. Points on the line x = 1 represent 
cases in which the HPAC contour envelops the corresponding NARAC contour, while 
points on the line y = 1 represent cases in which the NARAC contour envelops the 
corresponding HPAC contour. 

Figure 3-18 shows an example set of contours and overlap regions for four 
contour levels from run 1 at 30 minutes. The color convention for the contours is as 
before: solid red denotes NARAC contours and dashed blue denotes HPAC contours. 
The circles represent the boundary of the overlap region; red circles denote boundary 
points belonging to the NARAC contour, and blue circles denote boundary points 
belonging to the HPAC contour. Table 3-1 gives the contour areas, overlap areas, and the 
x and y components of the MOE for the four contours. 


6 Proper ordering of the polygon vertices is crucial for this routine to give sensible answers. Ordering is 
checked before submitting the vertices to polyarea. 


3-16 




Figure 3-18. Dosage Contours for MvM 1 at 30 Minutes: Contour Levels Are Indicated on 
the Left, Red (Blue) Lines Denote NARAC (HPAC) Contours, and Red (Blue) Circles Denote 
Boundary of Overlap Region Belonging to NARAC (HPAC) Contours 


Table 3-1. MOE Areas for MvM 1 at 30 Minutes 


Contour Level 
(kg s/m 3 ) 

HPAC 

Area 

(km 2 ) 

NARAC 

Area 

(km 2 ) 

Overlap 

Area 

(km 2 ) 

MOE (x) 

MOE (y) 

1.00E-10 

42.53 

44.05 

36.13 

0.820 

0.850 

1.00E-09 

34.79 

36.10 

29.47 

0.816 

0.847 

1.00E-08 

26.90 

27.55 

22.19 

0.806 

0.825 

1.00E-07 

18.74 

18.77 

14.86 

0.792 

0.793 

1.00E-06 

9.75 

9.68 

7.26 

0.750 

0.745 
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Figures 3-19 and 3-20 are plots of the area-based MOE for all gas runs at 30 and 
60 minutes, respectively. The marker size is proportional to the contour level. Grouped 
by stability class, the unstable points (green) are closest to the (1,1) point, again showing 
the best model agreement. HP AC contours tend to envelop NARAC contours at both 
times, particularly for the smaller dosages. NARAC contours envelop HP AC contours 
for some of the larger dosages. Points belonging to run 3 are indicated, and show 
extreme HP AC overprediction, as evidenced in the plots shown in Figures 3-14 through 
3-17. Omitting the results associated with the MvM 3 case, the remaining stable and 
neutral cases show a similar level of agreement, with somewhat better agreement for the 
neutral runs at 60 minutes. 


Gas Runs: Area MOE, 30 Minutes 



Figure 3-19. Area-Based MOE for All Gas Runs at 30 Minutes. Colors Denote Stability: 

Red = Stable, Blue = Neutral, and Green = Unstable. Circle Sizes are Proportional to 
Dosage Contour Level. By Stability Type, Unstable Cases Show Best Model Agreement. 
MvM 3 Points are Indicated, Showing Extreme HPAC Overprediction. 
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Gas Runs: Area MOE, 60 Minutes 



Figure 3-20. Area-Based MOE for All Gas Runs at 60 Minutes. Colors Denote Stability: 

Red = Stable, Blue = Neutral, and Green = Unstable. Circle Sizes are Proportional to 
Dosage Contour Level. By Stability Type, Unstable Cases Show Best Model Agreement. 
MvM 3 Points are Indicated, Showing Extreme HPAC Overprediction. 

B. PARTICLE RELEASES 

Eight particle releases were compared, all for near-neutral atmospheric stability. 
Particle size, release height, release duration, and release geometry were varied (see 
Table 2-1). Of these four variables, release height had the most dramatic effect on model 
agreement, with better agreement for lower height releases. The “very high” altitude 
release at 750 m (MvM 15) is considered separately from the other particle runs. Particle 
size has a moderate effect on model agreement, with smaller particles yielding somewhat 
better agreement. Release duration and release geometry did not induce major trends in 
model agreement. 
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Log(Dosage (kg s/m )) Crosswind (km) 


1. Influence of Particle Size on Comparisons 

Figures 3-21, 3-22, and 3-23 show dosage contours at 60 minutes for particle runs 
10, 13, and 16 respectively. 7 


MvMIO: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure 3-21. MvM 10 at 60 Minutes with Upper: NARAC (Red -) and HPAC (Blue - -) Dosage 
Contours; Lower: Log Dosages on Downwind Symmetry Axis. Particle Diameter is 5 p, 
Release Height is 2 m, and Release Duration is 15 min. 


7 The full collection of particle run plots can be found in Appendix D. 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM13: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure 3-22. MvM 13 at 60 Minutes with Upper: NARAC (Red -) and HPAC (Blue - -) Dosage 
Contours; Lower: Log Dosages on Downwind Symmetry Axis. Particle Diameter is 50 n, 
Release Height is 2 m, and Release Duration is “Instantaneous.” 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM16: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure 3-23. MvM 16 at 60 Minutes with Upper: NARAC (Red -) and HPAC (Blue - -) Dosage 
Contours; Lower: Log Dosages on Downwind Symmetry Axis. Particle Diameter (Mass 
Median Diameter) is 500 Release Height is 2 m, and Release Duration is 

“Instantaneous.” 
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The particle size for these three releases was 5 |i for MvM 10, 50 |i for MvM 13, 
and 500 (i for MvM 16. 8 The release height for all three runs was 2 meters. MvM 10 
was a continuous release (15 min, constant rate), while the other two releases were 
instantaneous. 9 All other input parameters for the three runs were the same. Overall 
agreement for the thee runs is very good; however, a close examination of the figures, 
particularly the symmetry axis dosage plots, shows that model agreement is best for the 5 
(i release (Figure 3-21). This is further supported by Figure 3-24, which shows the area- 
based MOE values for runs 10, 13, and 16 colored by particle size (with the symbol size 
proportional to dosage contour value). It is apparent from this figure that the 5 |l (blue) 
circles, as a group, are closest to the diagonal line; the 50 |i (red) and 500 |i (green) 
circles are progressively farther away. 10 The bulk of the circles are above the diagonal 
line, i.e., in the NARAC overprediction region; this is consistent with the contour plots 
shown in Figures 3-21 through 3-23. 


Particle Runs 10,13,16: Area MOE 



Figure 3-24. Area-Based MOE for Particle Runs 10,13, and 16 at 30 and 60 Minutes. 
Colors Indicate Particle Size and Circle Size is Proportional to Dosage Contour Level. 


° For MvM 16, a log-normal particle distribution with mass median diameter (MMD) = 500 (I and 
geometric standard deviation = a = 2 was used. 

9 Note the appearance of spatial oscillations in the HP AC dosage along the downwind symmetry axis for 
the instantaneous releases (Figures 3-22 and 3-23) and the absence of oscillations in the continuous 
release (Figure 3-21). This is discussed in more detail in Section C.3 below. 

10 We acknowledge that these trends could be an artifact of the small sample size examined. 
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Figure 3-25 and 3-26 present scatter plots of HP AC and NARAC downwind 
contour distances as a function of particle size for 30 and 60 minutes, respectively. At 
both times, the plots show that HP AC and NARAC downwind contour distances agree to 
within 20 percent. In general, HP AC overpredicts NARAC for this measure, with the 
results matching more closely with increasing downwind distance (decreasing dosage 
levels). 

For the higher dosages (shortest downwind distances) the 5-p points show the 
best agreement. All of the 16 5-p points (blue) shown in Figure 3-25 indicate 
HPAC/NARAC model agreement to within 20 percent (for this measure). At 60 minutes, 
2 of 15 5-p points indicate differences larger than 20 percent. These two points are 
associated with higher dosage contours (i.e., shorter downwind distances). For the 50 p 
particle runs, differences larger than 20 percent are indicated for 3 of 14 and 3 of 13 
dosage contours at 30 and 60 minutes, respectively. For the largest particles, the 
disagreement between the two models is greater, with the biggest differences reflecting 
NARAC overpredictions of HP AC at shorter distances and hence larger dosages (e.g., see 
the 10' 9 and 10' 10 dosage contours for MvM 16 in Figure 3-23). Differences larger than 
20 percent (for the 500 p particles) are shown for 4 of 11 and 5 of 11 dosage contours at 
30 and 60 minutes, respectively. 

At the longer downwind distances (i.e., lower dosage contours), the relative 
model behavior is similar for all particle sizes - HP AC slightly (within 20 percent) 
overpredicts NARAC with respect to downwind contour distance. In fact, for downwind 
distances >10 km, the 500 pdata show the best match (closest to the diagonal line) 
between HP AC and NARAC results. This result is true at both 30 and 60 minutes after 
the release. 

Further investigation is needed to determine the most important causes of the 
differences shown by the particle scenario comparisons. The observation that the greatest 
differences are seen in comparisons involving the largest particle sizes is consistent with 
the hypothesis that the differences are due to the algorithms used by NARAC and HP AC 
to compute settling velocity. The two models use different algorithms for determining 
the Reynolds number of the falling particle. HP AC computes the local air density, while 
NARAC uses a constant value of 1.225 kg/m 3 . 
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HPAC Dosage Contour Distance (km) 


Particle Runs: Downwind Contour Distances, 30 Minutes 



NARAC Dosage Contour Distance (km) 


Figure 3-25. Scatter Plot Comparison of Downwind Contour Distances at 30 Minutes for 
Particle Runs MvM 10 through 14,16, and 17: Colors Indicate Particle Size and Dashed 

Lines are 20% Overprediction Limits 
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HPAC Dosage Contour Distance (km) 


Particle Runs: Downwind Contour Distances, 60 Minutes 



NARAC Dosage Contour Distance (km) 


Figure 3-26. Scatter Plot Comparison of Downwind Contour Distances at 60 Minutes for 
Particle Runs MvM 10 through 14,16, and 17: Colors Indicate Particle Size and Dashed 

Lines are 20% Overprediction Limits 
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2. Influence of Release Height on Comparisons 

Release height plays a major role on model agreement for particle releases, with a 
trend for better agreement with lower release heights. Figures 3-27, 3-28, and 3-29 show 
dosage contours at 60 minutes for the three 5 p runs, with increasing release heights: run 
10 (release height = 2 m), run 11 (release height = 250 m) and run 12 (release height = 
400 m). The release duration for these 3 runs is 15 minutes. 


MvMIO: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Downwind (km) 


Figure 3-27. MvM 10 at 60 Minutes with Upper: NARAC (Red -) and HPAC (Blue - -) Dosage 
Contours; Lower: Log Dosages on Downwind Symmetry Axis. Particle Diameter is 5 p, 
Release Height is 2 m, and Release Duration is 15 min. 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM11: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure 3-28. MvM 11 at 60 Minutes with Upper: NARAC (Red -) and HPAC (Blue - -) Dosage 
Contours; Lower: Log Dosages on Downwind Symmetry Axis. Particle Diameter is 5 jx. 
Release Height is 250 m, and Release Duration is 15 min. 
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Log (Dosage (kg s/m )) Crosswind (km) 


MvM12: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure 3-29. MvM 12 at 60 Minutes with Upper: NARAC (Red -) and HPAC (Blue - -) Dosage 
Contours; Lower: Log Dosages on Downwind Symmetry Axis. Particle Diameter is 5 (i, 
Release Height is 400 m, and Release Duration is 15 min. 
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As the release height increases, the HP AC dosage contours tend to envelop the 
corresponding NARAC contours. This is shown by the examination of the area-based 
MOE for the three runs, as shown in Figure 3-30. The MOE values clearly separate by 
height, with the blue circles (2-meter release) showing the closest agreement, followed by 
the red circles (250-meter release), and finally the green circles (400-meter release). The 
trend for increasing HP AC overprediction (below the diagonal) with increasing release 
height is evident from the figure. NARAC tends to overpredict HP AC for the 2-meter 
release height, which suggests that there is an intermediate release height where the two 
models have maximal agreement. 


Particle Runs 10,11,12: Area MOE 



Figure 3-30. Area-Based MOE for Particle Runs 10,11, and 12 at 30 and 60 Minutes. 
Colors Indicate Release Height and Circle Size is Proportional to Dosage Contour Level. 
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3. High Altitude Release: Run 15 

The 50 |i particle release at a height of 750 meters, MvM 15, shows considerable 
model disagreement. This was the only release that was studied where the release height 
was greater than the boundary layer height (500 m). The run was chosen to have a 4-hour 
duration, in order to yield significant dosages at the 10-meter sampler height. Table 3-2 
shows a comparison of the maximum dosages for the two models at various times after 
the release, computed over the entire horizontal grid at the sampler height of 10 meters. 

Table 3-2. Maximum Dosages for MvM 15 at a Sampling Height of 10 Meters 


Time (min) 

Max HPAC Dosage 
(kg-s/m 3 ) 

Max NARAC Dosage 
(kg-s/m 3 ) 

30 

0 

0 

60 

1.2 x 10" 14 

0 

120 

4.6 x 10' 9 

1.2x1 O' 8 

180 

3.5 x 10' 8 

7.2 x 10 s 

240 

3.5 x 10' 8 

8.4 x 10 s 


At 30 minutes, neither model has a non-zero dosage at the 10-meter sampler 
height. At 60 minutes, the HP AC plume has descended to the sampling height, while the 
NARAC particles have not. At later times, the maximum NARAC dosage overtakes the 
maximum HP AC dosage and continues to increase, while the HP AC dosage levels off. 11 

Figures 3-31 and 3-32 show dosage contour plots 12 for MvM 15 at 120 and 240 
minutes, respectively (other figures for MvM 15 can be found in Appendix D). At 120 
minutes, all of the HP AC contours envelop the corresponding NARAC contours. At 240 
minutes, the higher-level contours compare more favorably; however, the low level 
NARAC contours are still enveloped by the HP AC contours. This trend is captured in 
the area MOE plot for MvM 15, shown in Figure 3-33. In the figure, symbol color 
denotes the time after the release (blue = 120, red = 180, and green = 240 minutes) and 
symbol size is proportional to contour level. Fewer contours were used at later times, due 
to the truncation of low-level contours at the horizontal domain boundary. It is clear 
from the figure that, at 120 minutes (blue), all of the HP AC contours surround the 


11 Recall, the dosage values shown in Table 3-2 reflect a 1 kg release. We expect that increasing the 
release mass would lead to proportionately increased dosages. 

12 Note that some low level contours are truncated at the boundary of the horizontal domain. This is an 
artifact of the chosen domain size and does not alter the conclusions drawn in this section. 
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Log(Dosage (kg s/m )) Crosswind (km) 


NARAC contours. At 180 minutes (red) and 240 minutes (green), the higher level 
NARAC and HP AC contours begin to intersect and the MOE values move away from the 
x = 1 axis. 


MvM15: Dosage (kg s/m 3 ), Height = 10 m, Time = 120 min. 




Figure 3-31. MvM 15 at 120 Minutes with Upper: Dosage Contours for HPAC (Blue - -) and 
NARAC (Red — ); Lower: Dosage on Downwind Symmetry Axis 
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Log (Dosage (kg s/m )) Crosswind (km) 


MvM15: Dosage (kg s/m 3 ), Height = 10 m, Time = 240 min. 




Downwind (km) 


Figure 3-32. MvM 15 at 240 Minutes with Upper: Dosage Contours for HPAC (Blue - -) and 
NARAC (Red — ); Lower: Dosage on Downwind Symmetry Axis 
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High Altitude Particle Run, MvM15: Area MOE 



Figure 3-33. Area-Based MOE for Particle MvM 15: Colors Indicate Time After Release and 
Circle Size is Proportional to Dosage Contour Level. Agreement is Better for Later Times 

and Higher Dosages. 


The trends observed for MvM 15 are similar to those observed in the other higher 
altitude releases near the boundary layer: MvM 2 with release height 50 m and MvM 3 
with release height 80 m, each with boundary layer height 100 m (recall Section A of this 
chapter). The cause of these differences is associated with the different vertical 
diffusivity parameterizations used by HPAC and NARAC. 

Figure 3-34 compares NARAC and HPAC instantaneous concentrations for MvM 
3 at 10-minute intervals after the release (t = 10, 20, 30, 40, 50 and 60 minutes). The 
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material in Figure 3-34 is moving from right to left. At each time there is a single HP AC 
contour (blue) and a single NARAC contour (red), and the concentration value for all 
contours is 10' kg/m . For times through 40 minutes, the horizontal positions and 
horizontal widths of the HP AC and NARAC concentration plumes agree well, however 
the vertical diffusion is much more pronounced in the HP AC plume. After 40 minutes 
the HP AC plume descends quickly (at 60 minutes the HP AC plume is on the ground near 
x = -27 km, and is barely visible), while the NARAC plume remains aloft. 


MvM03: Concentration (kg/m 3 ).Vertical Slice at y = 0 km 



Figure 3-34. NARAC and HPAC Instantaneous Concentrations (1.0x10 9 kg /m 3 ) for MvM 3 

at t = 10, 20, 30, 40, 50, and 60 Minutes. The Material is Moving from Right to Left. The 
HPAC Contours Show Much Greater Vertical Diffusion and Descend at a Greater Rate. 


The marked difference in the vertical dynamics for the two models points to a 
difference in the modeling of the vertical diffusivity. NARAC uses a significantly lower 
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free atmosphere diffusivity than HP AC [Ref. 3-3]. HP AC blends the free atmosphere 
and boundary layer diffusivities over roughly the top half of the boundary layer, while 
NARAC blends over a relative shallow layer at the boundary layer. Thus, a plausible 
hypothesis for the differences observed is that the 80 m release height (in the 100 m deep 
stable boundary layer) is associated with the vertical diffusivity blending layer (between 
the boundary layer and free atmosphere). 

It is difficult to characterize the turbulence in the stable boundary layer, especially 
above the surface and unfortunately there is no good data on elevated diffusion under 
stable conditions. There are well-known limitations to the similarity-theory turbulence 
scaling relationships that are used. Above the surface layer in the moderately stable 
(z;/L > 1) boundary layer, the local fluxes of heat and momentum, and not surface fluxes, 
can be the most important turbulence scaling parameters. However, local fluxes are often 
not available, and scaling relationships for the stable boundary layer based on surface 
fluxes are used. For very stable conditions (z/L > 6), turbulence may be intermittent in 
portions of the boundary layer, and no satisfactory theory has yet been developed. 
‘Upside down’ boundary layers, with higher turbulence aloft due to elevated wind 
maxima, are observed in very stable conditions, but not in less stable conditions. 

Comparison with the HP AC results also shows that, for releases occurring at a 
significant fraction of the boundary layer depth, the NARAC plumes show less horizontal 
(i.e., crosswind) diffusion than the HPAC results. This is most likely the result of 
different diffusivity parameterizations within the upper half of the boundary layer. (Note 
that this is in contrast to the results observed in the short-range, low-altitude release 
HPAC/NARAC Prairie Grass simulations, where the NARAC plumes are consistently 
wider than the HPAC plumes [Ref. 3-1].) 

C. EXCURSIONS 

In this final section of Chapter 3, we discuss the sensitivity of the model 
comparisons to input parameters that control spatial and temporal resolution. In the 
initial stages of this study, we arrived at the “baseline” parameters, outlined in Chapter 2, 
through a series of test comparisons, using a variety of parameter settings. This section 
serves to document some of the more important effects that were found in the course of 
these initial explorations. It also serves to underscore the need to understand and account 
for the intrinsic computational framework of the models that are being compared in a 
model-to-model study. 
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Three main effects were observed, two concerning HP AC spatial resolution and 
one concerning HP AC temporal resolution. 13 

1. HPAC Vertical Resolution 

Significant changes in HPAC downwind plume transport were observed when the 
HPAC vertical resolution was changed. This effect was most prominent in the low- 
altitude stable gas release, MvM 1. Figures 3-35 and 3-36 show the dosage contour 
comparisons with two values of the HPAC vertical resolution: a “coarse” value of 40 
meters (the default HPAC value for the specified grid top height) and a “fine” value of 10 
meters (the adopted “baseline” setting). The coarse resolution HPAC plume (blue), 
shown in Figure 3-35, is seen to propagate much farther downwind than the fine 
resolution HPAC plume, shown in Figure 3-36. For example, an examination of the 
HPAC dosage contour at 10~ 9 kg-s/m 3 shows that in the coarse resolution case, this 
contour extends approximately 37 km downwind, versus 29 km in the fine resolution 
case. 

The effect is caused by the inadequate resolution of the vertical wind profile for 
the coarse case, resulting in an artificially high transport speed at heights near the release 
height of 2 meters. 14 

Figure 3-37 shows the vertical wind profile for MvM 1 and a horizontal line at 40 
meters (the first sampling height for coarse sampling). Since the release occurs at 2 
meters, where the wind speed is approximately 3 m/s, and since the release is stable, 
implying a minimum of vertical mixing, it is reasonable to expect an effective transport 
speed near 3 m/s. It is clear from Figure 3-37 that sampling the vertical wind profile at a 
vertical resolution of 40 meters could artificially inflate the transport speed well beyond 3 
m/s, compared to sampling at 10-meter resolution. Although this effect was most 
dramatic in MvM 1, it was agreed to adopt the fine vertical resolution value of 10 meters 
as the baseline for all MvM runs (with the exception of MvM 15, which was run with 20- 
meter vertical resolution). This is a clear example of the need to move away from a 
default parameter setting in order to ensure a fair comparison of different model 
predictions. 


13 Related considerations apply to setting up NARAC meteorological and concentration grids but are not 
discussed here. (See Chapter 2 for additional discussion.) 

14 We note that the puffs use the wind speed at the puff centroid, but they assume a linear wind profile 
over the extent of the puff. The wind speed may vary significantly over the puff scale, hence requiring 
a finer resolution. 
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As noted above, the vertical resolution effect was most prominent for MvM 1. 
Recalling the vertical wind profiles of Figure 2-3, we can understand the origin of this 
relative prominence. First, the vertical wind speed profile for MvM 1 through 3 had the 
greatest speeds above the height of 10 meters, with a much more rapid increase in speed 
as a function of height than the profiles used for MvM 4 through 17. Therefore, the 
vertical profile associated with MvM 1 through 3 would be expected to be the most 
sensitive (i.e., artificially inflated) to the use of the lower (40-meter) vertical resolution. 


MvMOl: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Downwind (km) 


Figure 3-35. MvM 1 at 60 Minutes with Upper: Dosage Contours for HPAC (Blue - -) and 
NARAC (Red — ); Lower: Dosage on Downwind Symmetry Axis. HPAC Run Done at Coarse 
Vertical Resolution Value of 40 Meters (Default HPAC Value). Note Extended HPAC Plume 

in Downwind Direction. 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvMOl: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure 3-36. MvM 1 at 60 Minutes with Upper: Dosage Contours for HPAC (Blue - -) and 
NARAC (Red — ); Lower: Dosage on Downwind Symmetry Axis. HPAC Run Done at Fine 
Vertical Resolution Value of 10 Meters (Adopted Baseline HPAC Value). Effective 
Transport Speed of HPAC Plume is Lower With Finer Vertical Sampling. 
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Vertical Wind Profile for MvM 1 



Figure 3-37. Vertical Wind Profile for MvM 1 (Blue -) and 40-Meter Sampling Height 
(Red - -). Coarse Vertical Sampling Artificially Inflates the Effective Transport Speed. 


Also, recall that MvM 1, 2, and 3 had release heights of 2, 50, and 80 meters, 
respectively. Above 50 meters, all three vertical profiles that were considered were flat, 
that is, there was no change in wind speed with increasing altitude. Thus, only MvM l’s 
release occurred at a height where the associated wind speed was based on the first 
resolved 40-meter value (the significantly inflated value). For the conditions examined, 
the combination of low altitude release and steep vertical wind gradient led to a situation 
where substantially different predictions resulted from varying the vertical resolution. 

2. HPAC Horizontal Resolution 

A second resolution-induced effect was observed in HPAC dosage output, when 
the horizontal resolution (set in the HPAC Domain Editor) was changed from the default 
value of approximately 1 km to the finer value of 0.4 km. Figure 3-38 shows a screen 
capture of the HPAC surface dosage contours for MvM 1 at 60 minutes. The color-filled 
contours were computed at the finer 0.4-km horizontal resolution value, while the solid 
black contours span the same contour levels and were computed at the default horizontal 
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resolution value of 1 km. Both sets of contours were computed with a 40-meter vertical 
resolution (the HP AC default value, discussed above). 

Two differences between the contour sets are apparent. First, the downwind 
plume extent is reduced for the fine resolution contours, at least for the smaller dosages 
(e.g., 10' 7 and 10' 8 ). Next, there is the undesirable introduction of high-frequency spatial 
modulation in the fine resolution contours that distort the plume shape, particularly at 
larger dosages (see, for example, the boundary between the yellow and green regions). 


HP AC MvM 1, 0.4 km Horiz. Resolution 
Surface Dosage 

SF6 at 01-Apr-00 01:00:00Z (1.00 hrs) 



Figure 3-38. HPAC Surface Dosage Contours for MvM 1 at 60 Minutes. Color-Filled 
Contours Are For 0.4 km Horizontal Resolution and Solid Black Contours (at Same 
Contour Levels) Are For 1 km Horizontal Resolution. Note High-Frequency Distortions 
Introduced in Color Contours, Believed to Be an Artifact of “Puff-Splitting” Routines in 

SCIPUFF. 

The observed high-frequency spatial modulation is believed to be due to 
computational artifacts introduced by “puff-splitting” routines in SCIPUFF that control 
the creation of new Gaussian concentration puffs when plumes cross resolution cells. At 
the advice of the developer of SCIPUFF [Ref. 3-4], it was agreed to adopt the default 
HPAC horizontal resolution value of 1 km as the baseline value in all model-to-model 
comparisons. This choice, combined with the 10-meter vertical resolution choice, 
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produced HP AC contours that had acceptable transport speeds and acceptable (sensible) 
smoothness properties. 

It is interesting to note the asymmetry in the HP AC results with respect to 
changes in horizontal and vertical resolution. In the horizontal direction, coarsening the 
resolution improves the results by eliminating artificial high frequency modulation. In 
the vertical direction, making the resolution finer improves the comparative results with 
respect to transport speed. These effects were not anticipated and they underscore the 
complexities involved in making a fair model-to-model comparison. 

3. HPAC Time-Step Resolution 

The final parameter excursion investigated involves the HPAC maximum output 
time-step. This parameter is set within the Time Editor in the Project Editor, and, for all 
runs, the baseline value was chosen to be the default HPAC value of 60 seconds. The 
maximum time-step controls the time interval at which concentration values are listed in 
the HPAC sampler (.smp) file. The sampler file concentration data are then integrated, 
using Equation 2-2, to produce HPAC output dosages. 

It was noted in Section A.l that, for some releases, there are significant spatial 
oscillations 15 in the HPAC dosage near the release (see symmetry axis plots in Figures 
3-3 and 3-4 and the blue trace in Figure 3-39 below). These oscillations eventually decay 
with downwind distance (~ 1 to 3 km downwind). A survey of the downwind symmetry 
axis plots for all runs (see Appendices C and D) shows that the oscillations are present 
only for the 2-meter, instantaneous releases (MvM runs 1, 4, 7, 13, 16, and 17). 
Oscillations are present in both gas and particle runs, and at 30 and 60 minutes after the 
release. 

The magnitude of these oscillations is related to the output time-step chosen for 
this comparative study. Figure 3-39 shows the symmetry axis dosage for MvM 4 at 30 
minutes, computed with two time-steps: 60 seconds (the baseline value, in blue) and 30 
seconds (in red). Decreasing the time-step by a factor of two diminishes the magnitude 
of the oscillations by as much as a factor of 100. 


15 By viewing “raw” dosage values, it has been verified that the oscillations are not artifacts of the 
interpolation scheme used to compute dosages on the symmetry axis nor of the contouring routine used 
to generate dosage contours. 
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The SCIPUFF developer believes the dosage oscillations may be due to the 
following [Ref. 3-5]. At early times and large output time steps an instantaneous release 
plume can pass completely over an output sampler in one time step, with the result that 
the sampler records little or no concentration signal at that time step (this is especially 
true for the small release geometries considered in this study). Sampling the plume more 
frequently increases the chance that a sampler that previously recorded little or no 
concentration will now record a higher concentration. At later times (larger downwind 
distances), the plume has had time to grow dispersively and the undersampling effect is 
diminished. Future model-to-model studies that involve comparisons of short-range 
HP AC dosages would likely benefit from further reductions in the time-step. 

We note that the above oscillations are a consequence of our particular 
comparative methodology. For this study, we computed dosages by summing the output 
concentrations that were obtained from HPAC sampler files (i.e., Equation 2-2). We 
stress that the internally computed HPAC dosages, that are typically used to create a 
dosage plot, do not show these oscillations near the source 


MvM04, 30 Min: Dosage on Downwind Symmetry Axis 



Figure 3-39. HPAC Dosage Along Downwind Symmetry Axis for MvM 4 at 30 Minutes. 
Blue Trace Is for 60 Second Time-Step (Baseline Value) and Red Trace Is for 30 Second 
Time-Step. The Magnitude of the Dosage Oscillations Decreases Substantially by 

Decreasing the Time-Step. 
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APPENDIX A 
ACRONYMS 



2D 

ADAPT 

Ahpac 

Aol 

Aop 

Anarac 

ARAC 

ARAP 

Aup 

CB 

CBIAC 


deg C 
DOD 
DOE 
DTRA 

FOM 

GMU 

HP AC 
hr 

IDA 

kg 

km 

L 

LLNL 

LODI 

F 


APPENDIX A 
ACRONYMS 


T wo-dimensional 

Atmospheric Data Assimilation and Parameterization Techniques 

dosage region/area of an HP AC prediction 

region where model predictions agree 

region where model 1 overpredicts model 2 

dosage region/area of a NARAC prediction 

Atmospheric Release Advisory Center (same as NARAC) 

Aeronautical Research Associates of Princeton 

region where model 1 underpredicts model 2 

Chemical and Biological 

Chemical and Biological Defense Information Analysis Center 


degrees Celsius 
Department of Defense 
Department of Energy 
Defense Threat Reduction Agency 

Figure Of Merit 

George Mason University 

Hazard Prediction and Assessment Capability 
hour 

Institute for Defense Analyses 

kilogram 

kilometer 

Monin-Obukhov length 

Lawrence Livermore National Laboratory 

Lagrangian Operational Dispersion Integrator 

microns 


A-l 



m 

mg 

min 

MMD 

MOE 

m/s 

MvM 

NARAC 
NO A A 

y^m 

s 

SCIPUFF 

sec 

SF 6 or SF6 
t 

T&D 

u 

V&Y 

VV&A 

WMD 

z 

Zo 

zi or zi 
z r 

Zs 


meters 

milligram 

minutes 

Mass Median Diameter 
Measure of Effectiveness 
meters per second 
Model-Versus-Model 

National Atmospheric Release Advisory Center (same as ARAC) 
National Oceanic and Atmospheric Administration 

Stability Correction. 

seconds 

Second-Order Closure Integrated Puff 
seconds 

Sulfur Hexafluoride 
time 

Transport and Dispersion 

wind speed (m/s) 

Verification and Validation 
Verification, Validation, and Accreditation 

Weapons of Mass Destruction 

height (m) 

surface roughness factor (m) 
boundary layer height (m) 
release height (m) 
surface layer height (m) 
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APPENDIX B 

INITIAL CONDITIONS FOR HPAC AND NARAC PREDICTIONS 


This Appendix lists the input parameters and settings that were used for HPAC 
and NARAC for all seventeen model-to-model comparisons. For each of the seventeen 
runs, HPAC parameters are tabulated first, followed by NARAC “namelists” and other 
auxiliary NARAC input files. NARAC grid generation files are given at the end of this 
Appendix. 

The information provided in this appendix should allow for the reproduction of 
any and all of the model runs that were used for this study. Furthermore, starting with the 
input conditions described in this appendix, future analyses could easily extend these 
examinations by simply adjusting specific parameters. 
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HPAC MvMl 
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NARAC MvMl 


adapt_mvmO1.nml 

&adapt_control 
flag_map_adjust = .true. 
flag_debug = .true. 
opt_output_file = 'arac' 

/ 

&adapt_grid 

°pt_grid_file = 'gridgen' 

file_met_grid = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 

/ 

&adapt_metdata 
file_met_field 
opt_src_field 
opt_src_obs 
file_src_obs 
file_src_station 
flag_station_km 
nmethod 
/ 


&adapt_field2d 

zO = 0.008 

hgt_blend_layer = 10.0 

hgt_surface_layer = 10.0 

hgt_boundary_layer = 100.0 

inv_monin_obukhov_len = 0.04 

hgt_vert_coord = 'zAGL' 

/ 


= 'adapt_mtoml_2 0 0 0APR01_0 0 0 0 0 0.nc' 
= 'none ' 

= 'ascii2' 

= 'observ.met' 

== 'stnloc.met' 

= .false. 

= 2 


&adapt_method 
obs_date_time 
flag_mc_adjust 
flag_upr_in_sl 
opt_met_type 
opt_method 
opt_dist_wt 
opt_int_blend 
opt_int_sfc 
opt_int_upr 
opt_int_horz 
opt_vert_grid 
op t_wind_hor z 
opt_wind_vert 
max_veer_vert 
blend_max_veer 
/ 


'2000APR01_000000' 
.false. 

.true. 

' wind2d' 

' sparse_data' 

' inv_horz_dist_sq' 
'sfc_difference ' 

' log' 

'linear' 

'distance_wt' 

'sigmaZ' 

' uv' 

'spddir' 

240.0 
240.0 


&adapt_method 

opt_method 

opt_met_type 

obs_date_time 

/ 


' turb_nml' 

' turb_none' 

"200 0APR01 000000 


&adapt_turbulence 

turb_param_h = 'sigmav_simthry' 

turb_param_z = 'simthry' 


/ 


NARAC MvMl 

lodi files mvmO1.nml 
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&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/cone_1_40_- 

18000 0 001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvmOl/mgrid/adapt_mtoml_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 

NARAC MvMl 

lodi_mvm01.nml 

&prob_setup 
title 

tstart_str 
tstop_str 
dt_dump_str 
nbins 
nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 

&thist_param 

/ 

&src_param 

source_id = 'Source 1' 

geom_time_strs = "2000APR01_000000" 

max_num_part = 500000 

source_model = 'neutral' 

src_agl_flg = .true. 

geom_type = 2 

mean_x = 0.0 

mean_y = 0.0 

mean_z = 2.0 

std_x = 1.0 

std_y = 1.0 

std_z = 1.0 

cutoff_dx_min = 2.0 
cutoff_dy_min = 2.0 
cutoff_dz_min = 2.0 
cutoff_dx_max = 2.0 
cutoff_dy_max = 2.0 
cutoff dz max = 2.0 


= 'Model vs. Model: Source 1' 
= "200 0APR01_0 00000" 

= "200 0APR01_010000" 

= " 30 : 0 " 

= 2 
= 1 

= 'nongauss' 

= 1 

= "200 0APR01_0 00000" 

= 'arac' 
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start_time_str 

= 

"200 0APR01 000000" 



stop time str 

= 

"200 0APR01 000000" 



species 

= 

" SF6 " 



density 

= 

0.0 



er time strs 

= 

"200 0APR01 000000" 



emiss rates 

= 

1.0 



er_units_type 

= 

"mass" 



nset dep vel 

= 

0.0 



precip coeff 

= 

0.0 



m bin fract 

= 

1.0 



m bin diam min 

= 

o 

o 



m bin diam max 

= 

o 

o 



decay chain 

= 

.false. 



/ 





NARAC MvMl 





lodi mvmOl.nml (continued) 



&bin param 





bin id 

= 

'Bin 1: Integrated 

Air 

concentration xy' 

type 

= 

'air' 



source_list 

= 

'Source 1' 



species_name 

= 

" SF6" 



orientation 

= 

'xy' 



bin agl fig 

= 

.true. 



position 

= 

10.0 



width 

= 

20.0 



dt_samp_str 

= 

"30 : 00" 



dt_bin_out_str 

= 

"30 : 00" 



samp type 

= 

'integrated' 



/ 





&bin param 





bin id 

= 

'Bin 2: Integrated 

Air 

concentration xy' 

type 

= 

' air' 



source list 

= 

'Source 1' 



species_name 

= 

" SF6" 



orientation 

= 

'xy' 



bin agl fig 

= 

.true. 



position 

= 

10.0 



width 

= 

20.0 



dt samp str 

= 

o 

o 

o 

kO 



dt_bin_out_str 

= 

o 

o 

o 

KD 



samp type 

= 

'integrated' 



/ 





&turb_param 





turb param h 

= 

'sigmav simthry' 



turb param z 

= 

'simthry' 



read adapt turb 

/ 

= 

.true. 



&met_param 





/ 





NARAC MvMl 





stnloc mvmOl.met 





SFC 

'MtoMl' 0.0 0.0 : 

10 , 

. 0 
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UPR 


'MtoMl' 0.0 0.0 


NARAC MvMl 

observ mvmOl.met 


'MtoMl' 

METDATASET 

UPR 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 


0.008 

90 

2 0 0 0APR01 

0.5 

90 

1 

90 

1.5 

90 

2 

90 

2.5 

90 

3 

90 

3.5 

90 

4 

90 

4.5 

90 

5 

90 

5.5 

90 

6 

90 

6.5 

90 

7 

90 

7.5 

90 

8 

90 

8.5 

90 

9 

90 

9.5 

90 

10 

90 

10.5 

90 

11 

90 

11.5 

90 

12 

90 

12.5 

90 

13 

90 

13.5 

90 

14 

90 

14.5 

90 

15 

90 

15.5 

90 

16 

90 

16.5 

90 

17 

90 

17.5 

90 

18 

90 

18.5 

90 

19 

90 

19.5 

90 

20 

90 

20.5 

90 

21 

90 

21.5 

90 

22 

90 

22.5 

90 

23 

90 

23.5 

90 

24 

90 

24.5 

90 

25 

90 


0.0008 
000000' 

2.3467 
2.7835 
3.0606 
3.2724 
3.4484 
3.6017 
3.7394 
3.8657 
3.9832 
4.0938 
4.1988 
4.2993 
4.3959 
4.4891 
4.5796 
4.6675 
4.7533 
4.8372 
4.9194 
5 

5.0792 
5.1572 
5.234 
5.3098 
5.3846 
5.4577 
5.5288 
5.598 
5.6656 
5.7318 
5.7967 
5.8605 
5.9234 
5.9852 
6.0462 
6.1063 
6.1657 
6.2243 
6.2822 
6.3394 
6.396 
6.4519 
6.5071 
6.5618 
6.6159 
6.6694 
6.7223 
6.7747 
6.8266 
6.8779 
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observ mvmOl.met (continued) 


MtoMl' 

25.5 

90 

6.9287 

MtoMl' 

26 

90 

6.979 

MtoMl' 

26.5 

90 

7.0287 

MtoMl' 

27 

90 

7.078 

MtoMl' 

27.5 

90 

7.1268 

MtoMl' 

28 

90 

7.1751 

MtoMl' 

28.5 

90 

7.223 

MtoMl' 

29 

90 

7.2704 

MtoMl' 

29.5 

90 

7.3173 

MtoMl' 

30 

90 

7.3638 

MtoMl' 

30.5 

90 

7.4098 

MtoMl' 

31 

90 

7.4554 

MtoMl' 

31.5 

90 

7.5006 

MtoMl' 

32 

90 

7.5454 

MtoMl' 

32.5 

90 

7.5898 

MtoMl' 

33 

90 

7.6337 

MtoMl' 

33.5 

90 

7.6773 

MtoMl' 

34 

90 

7.7204 

MtoMl' 

34.5 

90 

7.7632 

MtoMl' 

35 

90 

7.8056 

MtoMl' 

35.5 

90 

7.8476 

MtoMl' 

36 

90 

7.8893 

MtoMl' 

36.5 

90 

7.9306 

MtoMl' 

37 

90 

7.9715 

MtoMl' 

37.5 

90 

8.0121 

MtoMl' 

38 

90 

8.0523 

MtoMl' 

38.5 

90 

8.0922 

MtoMl' 

39 

90 

8.1317 

MtoMl' 

39.5 

90 

8.171 

MtoMl' 

40 

90 

8.2099 

MtoMl' 

40.5 

90 

8.2484 

MtoMl' 

41 

90 

8.2867 

MtoMl' 

41.5 

90 

8.3246 

MtoMl' 

42 

90 

8.3623 

MtoMl' 

42.5 

90 

8.3996 

MtoMl' 

43 

90 

8.4366 

MtoMl' 

43.5 

90 

8.4734 

MtoMl' 

44 

90 

8.5098 

MtoMl' 

44.5 

90 

8.546 

MtoMl' 

45 

90 

8.5819 

MtoMl' 

45.5 

90 

8.6175 

MtoMl' 

46 

90 

8.6528 

MtoMl' 

46.5 

90 

8.6879 

MtoMl' 

47 

90 

8.7227 

MtoMl' 

47.5 

90 

8.7572 

MtoMl' 

48 

90 

8.7915 

MtoMl' 

48.5 

90 

8.8255 

MtoMl' 

49 

90 

8.8593 

MtoMl' 

49.5 

90 

8.8928 

MtoMl' 

50 

90 

8.9261 
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NARAC MvM2 


adapt_mvm02.nml (same as adapt_mvm01.nml) 

lodi_files_mvm02.nml 

&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/conc_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvmOl/mgrid/adapt_mtoml_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 

NARAC MvM2 

lodi_mvm02 .nml 

&prob_setup 
title 

tstart_str 
tstop_str 
dt_dump_str 
nbins 
nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 

&thist_param 

/ 

Sesrc_param 

source_id = 'Source 1' 

geom_time_strs = "200 0APR01_0 00000" 

max_num_part = 500000 

source_model = 'neutral' 

src_agl_flg = .true. 

geom_type = 2 

mean_x = 0.0 

mean_y = 0.0 

mean_z = 50.0 

std_x = 1.0 

std_y = 1.0 


= 'Model vs. Model: Source 2' 
= "200 0APR01_0 00000" 

= "200 0APR01_010000" 

= "30:0" 

= 2 
= 1 

= 'nongauss' 

= 1 

= "200 0APR01_0 00000" 

= 'arac' 
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std_z 

cutof f_dx_min 
cutof f_dy_min 
cutof f_dz_min 
cutof f_dx_max 
cutof f_dy_max 
cutof f_dz_max 

start_time_str 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min 

m_bin_diam_max 

decay_chain 


1.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 

"200 0APR01_0 00000 
"200 0APR01_0 00000 
" SF6 " 

6.035 

"200 0APR01_0 00000 

1.0 

"mass" 

0.0 
0.0 
1.0 
0.0 
0.0 

.false. 


NARAC MvM2 

lodi mvm02.nml (continued) 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 
dt_bin_out_str = 
samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 
dt_bin_out_str = 
samp_type 

/ 

&turb_param 
turb_param_h 
turb_param_z 
read_adapt_turb = 


"Bin 1: Integrated Ai 
' air" 

"Source 1' 

" SF6" 

"xy" 

.true. 

10.0 
20.0 
"30:00" 

"30:00" 

"integrated" 


"Bin 2: Integrated Ai 
"air" 

'Source 1' 

" SF6" 

"xy" 

.true. 

10.0 
20.0 
" 60 : 00 " 

"60:00" 

"integrated" 


'sigmav_simthry" 
'simthry" 

. true . 


concentration xy" 


concentration xy" 


B-10 



/ 

&met_param 

/ 

NARAC MvM2 

stnloc_mvm02.met (same as stnlocmvmOl.met) 
observ mvm02.met (same as observ mvmOl.met) 
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HPAC MvM3 


mvm3 

HPAC 3.2 parameter 

HPAC 3.2 value 

HPAC SCREEN 

agent 

SF6 

New Project Setup 

project file name 

mvm3.pri 


coordinates 

cartesian 


local oriqinx 

0 km 


local origin y 

0 km 


local origin latitude 

45 deqN 


local origin longitude 

Odeq 


reference times 

UTC 


local time of 00:00Z 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off 

New Roject Editor 

edit mode 

advanced 

Material Editor 

type 

qas 


name 

SF6 


material file 

Sf6.mtl 


units 

kg 


bin boundaries 

NA 


output flags 

surf, dosaqe, surf, depos. 


daytime decay rate 

0 


nighttime decay rate 

0 

Gas Parameter Editor 

density ratio (rho/riio air) 

5.02917 


deposition veloc ity 

0 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


height of release 

80 m 


source uncertainty 

no 


spec ification 

simple 


randomize location 

no 


release duration 

instantaneous 


release rate 

NA 


agent mass 

1 kq 


xsize 

1 m 


y size 

lm 


zsize 

1 m 


momentum 

0 


buoyancy 

0 

lime Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

100 


maximum time step 

60 sec 


output interval 

15 min 

Domain Editor 

xmin domain 

-50 km 


xmax domain 

1 km 


ymin domain 

-5 km 


ymax domain 

5 km 


horizontal resolution 

default 


vertical domain max height 

500 m 


vertical domain resolution 

10 m 

Weather Editor 

weatherdata type 

surface obs. only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability length scale 

NA 


large scale variability variance 

NA 


surface toughness 

0.008 m 


precipitation 

none 


obs. time bin size 

1 hr 

Meteorology Option Editor 

save meteorology fields 

no 


terrain/land cover file 

no 

Surface Observations File 

file name 

mvmSsfc 


elevation 

10 m 


zi (boundary layer height) 

100 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Length) 

25 m 

Options Editor 

puff split grid level 

2 


surface resolution 

default 


puff qrid resolution 

0 


boundary layer pts 

11 


stable atm. turbulence 

le-2 rrT2/s~2 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 m~2/s~2 


calm cond. scale 

1000 m 


surface dosage height 

0 


minimum puff mass 

le-20 


conditional averaging time 

1 hr 


lumped boundary layer 

no 
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NARAC MvM3 


adapt_mvm03.nml (same as adapt_mvm01.nml) 

lodi_files_mvm03.nml 

&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/conc_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvmOl/mgrid/adapt_mtoml_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 

NARAC MvM3 

lodi_mvm03.nml 

&prob_setup 
title 

tstart_str 
tstop_str 
dt_dump_str 
nbins 
nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 

&thist_param 

/ 

Sesrc_param 

source_id = 'Source 1' 

geom_time_strs = "2000APR01_000000" 

max_num_part = 500000 

source_model = 'neutral' 

src_agl_flg = .true. 

geom_type = 2 

mean_x = 0.0 

mean_y = 0.0 

mean_z = 80.0 

std_x = 1.0 

std_y = 1.0 

std z = 1.0 


= 'Model vs. Model: Source 3' 
= "200 0APR01_0 00000" 

= "200 0APR01_010000" 

= "30:0" 

= 2 
= 1 

= 'nongauss' 

= 1 

= "200 0APR01_0 00000" 

= 'arac' 
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cutof f_dx_min 
cutof f_dy_min 
cutof f_dz_min 
cutof f_dx_max 
cutof f_dy_max 
cutof f_dz_max 

start_time_str 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min 

m_bin_diam_max 

decay_chain 

/ 


2.0 
2.0 
2.0 
2.0 
2.0 
2.0 

"200 0APR01_0 00000 
"200 0APR01_0 00000 
" SF6 " 

6.035 

"200 0APR01_0 00000 
1.0 

"mass" 

0.0 
0.0 
1.0 
0.0 
0.0 

.false. 


NARAC MvM3 

lodi_mvm03.nml (continued) 

&bin_param 


bin id 

= 

'Bin 1: Integrated 

type 

= 

' air' 

source_list 

= 

'Source 1' 

species name 

= 

" SF6" 

orientation 

= 

'xy' 

bin agl fig 

= 

.true. 

position 

= 

10.0 

width 

= 

20.0 

dt_samp_str 

= 

"30 : 00" 

dt bin out str 

= 

"30 : 00" 

samp_type 

= 

'integrated' 

• in param 

bin id 

= 

'Bin 2: Integrated 

type 

= 

'air' 

source_list 

= 

'Source 1' 

species_name 

= 

" SF6" 

orientation 

= 

'xy' 

bin agl fig 

= 

.true. 

position 

= 

10.0 

width 

= 

20.0 

dt_samp_str 

= 

o 

o 

o 

KD 

dt_bin_out_str 

= 

O 

O 

O 

KD 

samp type 

= 

'integrated' 

urb param 

turb param h 

= 

'sigmav simthry' 

turb param z 

= 

'simthry' 

read adapt turb 

= 

.true. 


&met_param 

/ 
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NARAC MvM3 


stnloc_mvm03.met (same as stnloc_mvm01 .met) 
observ mvm03.met (same as observ mvmOl.met) 
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HPAC MvM4 


mvm4 

HPAC 3.2 parameter 

HPAC 3.2 value 

HPAC SCREEN 

agent 

SF6 

New Project Setup 

project file name 

mvm4.pri 


coordinates 

cartesian 


local oriqinx 

0 km 


local origin y 

0 km 


local origin latitude 

45 deqN 


local origin longitude 

Odeq 


reference times 

UTC 


local time of 00:00Z 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off 

New Roject Editor 

edit mode 

advanced 

Material Editor 

type 

qas 


name 

SF6 


material file 

Sf6.mtl 


units 

kg 


bin boundaries 

NA 


output flags 

surf, dosaqe, surf, depos. 


daytime decay rate 

0 


nighttime decay rate 

0 

Gas Parameter Editor 

density ratio (rho/riio air) 

5.02917 


deposition veloc ity 

0 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


height of release 

2 m 


source uncertainty 

no 


spec ification 

simple 


randomize location 

no 


release duration 

instantaneous 


release rate 

NA 


agent mass 

1 kq 


xsize 

1 m 


y size 

lm 


zsize 

1 m 


momentum 

0 


buoyancy 

0 

lime Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

100 


maximum time step 

60 sec 


output interval 

15 min 

Domain Editor 

xmin domain 

-50 km 


xmax domain 

1 km 


ymin domain 

-5 km 


ymax domain 

5 km 


horizontal resolution 

default 


vertical domain max height 

500 m 


vertical domain resolution 

10 m 

Weather Editor 

weatherdata type 

surface obs. only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability length scale 

NA 


large scale variability variance 

NA 


surface toughness 

0.008 m 


precipitation 

none 


obs. time bin size 

1 hr 

Meteorology Option Editor 

save meteorology fields 

no 


terrain/land cover file 

no 

Surface Observations File 

file name 

mvmN.sfc 


elevation 

10 m 


zi (boundary layer height) 

500 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Length) 

-500 m 

Options Editor 

puff split grid level 

2 


surface resolution 

default 


puff qrid resolution 

0 


boundary layer pts 

11 


stable atm. turbulence 

le-2 m~2/s~2 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 m~2/s~2 


calm cond. scale 

1000 m 


surface dosage height 

0 


minimum puff mass 

le-20 


conditional averaging time 

1 hr 


lumped boundary layer 

no 
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NARAC MvM4 


adapt_mvm04.nml 

&adapt_control 
flag_map_adjust = .true. 
flag_debug = .true. 
opt_output_file = 'arac' 

/ 

&adapt_grid 

°pt_grid_file = 'gridgen' 

file_met_grid = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 

/ 


&adapt_metdata 
file_met_field 
opt_src_field 
opt_src_obs 
file_src_obs 
file_src_station 
flag_station_km 
nmethod 


/ 


'adapt_mtom4_2 0 0 OAPRO1_0 0 0 0 0 0.nc' 
'none' 

'ascii2' 

'observ.met' 

' stnloc.met' 

.false. 

2 


&adapt_field2d 

zO = 0.008 

hgt_blend_layer = 10.0 

hgt_surface_layer = 10.0 

hgt_boundary_layer = 500.0 

inv_monin_obukhov_len = -0.002 

hgt_vert_coord = 'zAGL' 

/ 


&adapt_method 
obs_date_time 
flag_mc_adjust 
flag_upr_in_sl 
opt_met_type 
opt_method 
opt_dist_wt 
opt_int_blend 
opt_int_sfc 
opt_int_upr 
opt_int_horz 
opt_vert_grid 
op t_wind_hor z 
opt_wind_vert 
max_veer_vert 
blend_max_veer 
/ 


'2000APR01_000000' 
.false. 

.true. 

' wind2d' 

' sparse_data' 

' inv_horz_dist_sq' 
'sfc_difference ' 

' log' 

'linear' 

'distance_wt' 

'sigmaZ' 

' uv' 

'spddir' 

240.0 
240.0 


&adapt_method 

opt_method 

opt_met_type 

obs_date_time 

/ 


' turb_nml' 

'turb_none' 

"200 OAPRO1 000000" 


&adapt_turbulence 

turb_param_h = 'sigmav_simthry' 

turb_param_z = 'simthry' 

/ 


NARAC MvM4 


lodi_files_mvm04.nml 

&grid_name 

num_m_grids = 1 
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/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc 


m_grid_name = ' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/cone_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvm04/mgrid/adapt_mtom4_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 


NARAC \Iv\I4 


lodi mvm04 . nml 

&prob_setup 
title 

tstart_str 
tstop_str 
dt_dump_str 
nbins 
nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 

&thist_param 

/ 

&src_param 

source_id = 'Source 1' 

geom time strs = "2000APR01_000000" 
= 500000 
= 'neutral' 

= .true. 

= 2 
= 0.0 
= 0.0 
= 2.0 
= 1.0 
= 1.0 
= 1.0 
= 2.0 
= 2.0 
= 2.0 
= 2.0 
= 2.0 
= 2.0 


max_num_p art 
source_model 
src_agl_flg 

geom_type 

mean_x 

mean_y 

mean_z 

std_x 

std_y 

std_z 

cutof f_dx_min 
cutof f_dy_min 
cutof f_dz_min 
cutof f_dx_max 
cutof f_dy_max 
cutoff dz max 


= 'Model vs. Model: Source 4' 
= "200 0APR01_0 00000" 

= "200 0APR01_010000" 

=" 30 : 0 " 

= 2 
= 1 

= 'nongauss' 

= 1 

= "200 0APR01_0 00000" 

= 'arac' 
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"200 0APR01_0 00000" 
"200 OAPRO1_0 00000" 
" SF6 " 

6.035 

"200 OAPRO1_0 00000" 
1.0 

"mass" 

0.0 
0.0 
1.0 
0.0 
0.0 

.false. 


NARAC MvM4 

lodi_mvm04.nml (continued) 

&bin_param 
bin_id 
type 

source_list 
species_name 
orientation 
bin_agl_flg 
position 
width 

dt_samp_str 
dt_bin_out_str 
samp_type 

/ 

&bin_param 
bin_id 
type 

source_list 
species_name 
orientation 
bin_agl_flg 
position 
width 

dt_samp_str 
dt_bin_out_str 
samp_type 

/ 

&turb_param 

turb_param_h = 'sigmav_simthry' 
turb_param_z = 'simthry' 
read_adapt_turb = .true. 

/ 

&met_param 

/ 


NARAC MvM4 

stnloc_mvm04.met 

SFC 


= 'Bin 1: Integrated Air concentration xy 
= 'air' 

= 'Source 1' 

= "SF6" 

= 'xy' 

= .true. 

= 10.0 
= 20.0 
= "30:00" 

= "30:00" 

= 'integrated' 


= 'Bin 2: Integrated Air concentration xy 
= 'air' 

= 'Source 1' 

= "SF6" 

= 'xy' 

= .true. 

= 10.0 
= 20.0 
= "60:00" 

= "60:00" 

= 'integrated' 


start_time_str = 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min = 

m_bin_diam_max = 

decay_chain 
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'MtoMl' 0.0 0.0 10.0 

UPR 

'MtoMl' 0.0 0.0 


NARAC MvM4 


observ mvm04.met 


METDATASET '2000APR01_000000' 
UPR 


MtoMl' 

0.5 

90 

2.9251 

MtoMl' 

1 

90 

3.4132 

MtoMl' 

1.5 

90 

3.6977 

MtoMl' 

2 

90 

3.8988 

MtoMl' 

2.5 

90 

4.0543 

MtoMl' 

3 

90 

4.1808 

MtoMl' 

3.5 

90 

4.2875 

MtoMl' 

4 

90 

4.3795 

MtoMl' 

4.5 

90 

4.4605 

MtoMl' 

5 

90 

4.5326 

MtoMl' 

5.5 

90 

4.5977 

MtoMl' 

6 

90 

4.6569 

MtoMl' 

6.5 

90 

4.7112 

MtoMl' 

7 

90 

4.7613 

MtoMl' 

7.5 

90 

4.8078 

MtoMl' 

8 

90 

4.8511 

MtoMl' 

8.5 

90 

4.8918 

MtoMl' 

9 

90 

4.9299 

MtoMl' 

9.5 

90 

4.9659 

MtoMl' 

10 

90 

5 

MtoMl' 

10.5 

90 

5.0323 

MtoMl' 

11 

90 

5.063 

MtoMl' 

11.5 

90 

5.0923 

MtoMl' 

12 

90 

5.1202 

MtoMl' 

12.5 

90 

5.147 

MtoMl' 

13 

90 

5.1726 

MtoMl' 

13.5 

90 

5.1971 

MtoMl' 

14 

90 

5.2208 

MtoMl' 

14.5 

90 

5.2435 

MtoMl' 

15 

90 

5.2654 

MtoMl' 

15.5 

90 

5.2865 

MtoMl' 

16 

90 

5.3069 

MtoMl' 

16.5 

90 

5.3267 

MtoMl' 

17 

90 

5.3457 

MtoMl' 

17.5 

90 

5.3642 

MtoMl' 

18 

90 

5.3822 

MtoMl' 

18.5 

90 

5.3996 

MtoMl' 

19 

90 

5.4164 

MtoMl' 

19.5 

90 

5.4329 

MtoMl' 

20 

90 

5.4488 

MtoMl' 

20.5 

90 

5.4643 

MtoMl' 

21 

90 

5.4795 

MtoMl' 

21.5 

90 

5.4942 

MtoMl' 

22 

90 

5.5085 

MtoMl' 

22.5 

90 

5.5225 

MtoMl' 

23 

90 

5.5362 

MtoMl' 

23.5 

90 

5.5495 

MtoMl' 

24 

90 

5.5626 

MtoMl' 

24.5 

90 

5.5753 

MtoMl' 

25 

90 

5.5878 

MtoMl' 

25.5 

90 

5.5999 
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NARAC MvM4 

observ mvm04.met (continued) 


MtoMl' 

26 

90 

5.6118 

MtoMl' 

26.5 

90 

5.6235 

MtoMl' 

27 

90 

5.6349 

MtoMl' 

27.5 

90 

5.6461 

MtoMl' 

28 

90 

5.6571 

MtoMl' 

28.5 

90 

5.6678 

MtoMl' 

29 

90 

5.6784 

MtoMl' 

29.5 

90 

5.6887 

MtoMl' 

30 

90 

5.6989 

MtoMl' 

30.5 

90 

5.7088 

MtoMl' 

31 

90 

5.7186 

MtoMl' 

31.5 

90 

5.7282 

MtoMl' 

32 

90 

5.7376 

MtoMl' 

32.5 

90 

5.7469 

MtoMl' 

33 

90 

5.756 

MtoMl' 

33.5 

90 

5.765 

MtoMl' 

34 

90 

5.7738 

MtoMl' 

34.5 

90 

5.7824 

MtoMl' 

35 

90 

5.791 

MtoMl' 

35.5 

90 

5.7993 

MtoMl' 

36 

90 

5.8076 

MtoMl' 

36.5 

90 

5.8157 

MtoMl' 

37 

90 

5.8237 

MtoMl' 

37.5 

90 

5.8316 

MtoMl' 

38 

90 

5.8394 

MtoMl' 

38.5 

90 

5.847 

MtoMl' 

39 

90 

5.8545 

MtoMl' 

39.5 

90 

5.862 

MtoMl' 

40 

90 

5.8693 

MtoMl' 

40.5 

90 

5.8765 

MtoMl' 

41 

90 

5.8836 

MtoMl' 

41.5 

90 

5.8906 

MtoMl' 

42 

90 

5.8975 

MtoMl' 

42.5 

90 

5.9044 

MtoMl' 

43 

90 

5.9111 

MtoMl' 

43.5 

90 

5.9178 

MtoMl' 

44 

90 

5.9243 

MtoMl' 

44.5 

90 

5.9308 

MtoMl' 

45 

90 

5.9372 

MtoMl' 

45.5 

90 

5.9435 

MtoMl' 

46 

90 

5.9497 

MtoMl' 

46.5 

90 

5.9559 

MtoMl' 

47 

90 

5.962 

MtoMl' 

47.5 

90 

5.968 

MtoMl' 

48 

90 

5.9739 

MtoMl' 

48.5 

90 

5.9798 

MtoMl' 

49 

90 

5.9856 

MtoMl' 

49.5 

90 

5.9913 

MtoMl' 

50 

90 

5.997 
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HPAC MvM5 


mvm5 

HPAC 3.2 parameter 

HPAC 3.2 value 

HPAC SCREEN 

agent 

SF6 

New Project Setup 

project file name 

mvml l.prj 


coordinates 

cartesian 


local oriqinx 

0 km 


local origin y 

0 km 


local origin latitude 

45 deqN 


local origin longitude 

Odeq 


reference times 

UTC 


local time of 00:00Z 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off 

New Roject Editor 

edit mode 

advanced 

Material Editor 

type 

qas 


name 

SF6 


material file 

Sf6.mtl 


units 

kg 


bin boundaries 

NA 


output flags 

surf, dosaqe, surf, depos. 


daytime decay rate 

0 


nighttime decay rate 

0 

Gas Parameter Editor 

density ratio (rho/riio air) 

5.02917 


deposition veloc ity 

0 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


height of release 

250 m 


source uncertainty 

no 


spec ification 

simple 


randomize location 

no 


release duration 

instantaneous 


release rate 

NA 


agent mass 

1 kq 


xsize 

1 m 


y size 

lm 


zsize 

1 m 


momentum 

0 


buoyancy 

0 

lime Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

100 


maximum time step 

60 sec 


output interval 

15 min 

Domain Editor 

xmin domain 

-50 km 


xmax domain 

1 km 


ymin domain 

-5 km 


ymax domain 

5 km 


horizontal resolution 

default 


vertical domain max height 

500 m 


vertical domain resolution 

10 m 

Weather Editor 

weatherdata type 

surface obs. only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability length scale 

NA 


large scale variability variance 

NA 


surface toughness 

0.008 m 


precipitation 

none 


obs. time bin size 

1 hr 

Meteorology Option Editor 

save meteorology fields 

no 


terrain/land cover file 

no 

Surface Observations File 

file name 

mvmN.sfc 


elevation 

10 m 


zi (boundary layer height) 

500 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Length) 

-500 m 

Options Editor 

puff split grid level 

2 


surface resolution 

default 


puff qrid resolution 

0 


boundary layer pts 

11 


stable atm. turbulence 

le-2 m~2/s~2 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 m~2/s~2 


calm cond. scale 

1000 m 


surface dosage height 

0 


minimum puff mass 

le-20 


conditional averaging time 

1 hr 


lumped boundary layer 

no 
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NARAC MvM5 


adapt_mvm05 .nml (same as adapt_mvm04.nml) 

lodi_files_mvm05.nml 

&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/cone_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvm04/mgrid/adapt_mtom4_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 


NARAC \Iv\I5 

lodi mvmO5 . nml 


&prob_setup 

title 

tstart_str 

tstop_str 

dt_dump_str 

nbins 

nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met format 


/ 


'Model vs. Model: Source 5' 
"200 0APR01_0 00000" 

"200 0APR01_010000" 

"30:0" 

2 

1 

7 nongauss' 

1 

"200 0APR01_0 00000" 

'arac ' 


&thist_param 

/ 


Sc src_param 
source_id 
geom_time_strs 
max_num_p art 
source_model 
src_agl_flg 

geom_type 

mean_x 

mean_y 

mean_z 

std_x 

std_y 


'Source 1' 

"200 0APR01_0 00000" 
500000 
'neutral' 

.true. 

2 

0.0 
0.0 
250.0 
1.0 
1.0 
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std_z 

cutof f_dx_min 
cutof f_dy_min 
cutof f_dz_min 
cutof f_dx_max 
cutof f_dy_max 
cutof f_dz_max 

start_time_str 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min 

m_bin_diam_max 

decay_chain 


1 . 0 
2 . 0 
2 . 0 
2. 0 
2. 0 
2. 0 
2. 0 

"200 0APR01_0 00000 
"200 0APR01_0 00000 
" SF6 " 

6.035 

"200 0APR01_0 00000 

1.0 

"mass" 

0 . 0 
0 . 0 
1 . 0 
0 . 0 
0 . 0 

.false. 


NARAC MvM5 

lodi mvm05.nml (continued) 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 
dt_bin_out_str = 
samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 
dt_bin_out_str = 
samp_type 

/ 

&turb_param 
turb_param_h 
turb_param_z 
read_adapt_turb = 

/ 


'Bin 1: Integrated 
' air' 

'Source 1' 

" SF6" 

'xy' 

.true. 

10.0 
20.0 
"30:00" 

"30:00" 

'integrated' 


'Bin 2: Integrated 
' air' 

'Source 1' 

" SF6" 

'xy' 

.true. 

10.0 
20.0 
" 60 : 00 " 

"60:00" 

'integrated' 


'sigmav_simthry' 
'simthry' 

. true . 


Air concentration xy 


Air concentration xy 


&met_param 

/ 
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NARAC \Iv\I5 

stnloc_mvm05.met (same as stnloc_mvm04 .met) 
observ mvm05.met (same as observ mvm04.met) 
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HPAC MvM6 


mvm6 

HPAC 3.2 parameter 

HPAC 3.2 value 

HPAC SCREEN 

agent 

SF6 

New Project Setup 

project file name 

mvm6.pii 


coordinates 

cartesian 


local oriqinx 

0 km 


local origin y 

0 km 


local origin latitude 

45 deqN 


local origin longitude 

Odeq 


reference times 

UTC 


local time of 00:00Z 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off 

New Roject Editor 

edit mode 

advanced 

Material Editor 

type 

qas 


name 

SF6 


material file 

Sf6.mtl 


units 

kg 


bin boundaries 

NA 


output flags 

surf, dosaqe, surf, depos. 


daytime decay rate 

0 


nighttime decay rate 

0 

Gas Parameter Editor 

density ratio (rho/riio air) 

5.02917 


deposition veloc ity 

0 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


height of release 

400 m 


source uncertainty 

no 


spec ification 

simple 


randomize location 

no 


release duration 

instantaneous 


release rate 

NA 


agent mass 

1 kq 


xsize 

1 m 


y size 

lm 


zsize 

1 m 


momentum 

0 


buoyancy 

0 

lime Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

100 


maximum time step 

60 sec 


output interval 

15 min 

Domain Editor 

xmin domain 

-50 km 


xmax domain 

1 km 


ymin domain 

-5 km 


ymax domain 

5 km 


horizontal resolution 

default 


vertical domain max height 

500 m 


vertical domain resolution 

10 m 

Weather Editor 

weatherdata type 

surface obs. only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability length scale 

NA 


large scale variability variance 

NA 


surface toughness 

0.008 m 


precipitation 

none 


obs. time bin size 

1 hr 

Meteorology Option Editor 

save meteorology fields 

no 


terrain/land cover file 

no 

Surface Observations File 

file name 

mvmN.sfc 


elevation 

10 m 


zi (boundary layer height) 

500 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Length) 

-500 m 

Options Editor 

puff split grid level 

2 


surface resolution 

default 


puff qrid resolution 

0 


boundary layer pts 

11 


stable atm. turbulence 

le-2 m~2/s~2 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 m~2/s~2 


calm cond. scale 

1000 m 


surface dosage height 

0 


minimum puff mass 

le-20 


conditional averaging time 

1 hr 


lumped boundary layer 

no 
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NARAC MvM6 


adapt_mvm06.nml (same as adapt_mvm04.nml) 

lodi_filesmvmOS.nml 

&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/conc_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvm04/mgrid/adapt_mtom4_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 

NARAC MvM6 

lodi_mvm06.nml 

&prob_setup 
title 

tstart_str 
tstop_str 
dt_dump_str 
nbins 
nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 

&thist_param 

/ 

Sesrc_param 

source_id = 'Source 1' 

geom_time_strs = "200 0APR01_0 00000" 

max_num_part = 500000 

source_model = 'neutral' 

src_agl_flg = .true. 

geom_type = 2 

mean_x = 0.0 

mean_y = 0.0 

mean_z = 400.0 

std_x = 1.0 

std_y = 1.0 


= 'Model vs. Model: Source 6' 
= "200 0APR01_0 00000" 

= "200 0APR01_010000" 

= "30:0" 

= 2 
= 1 

= 'nongauss' 

= 1 

= "200 0APR01_0 00000" 

= 'arac' 
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std_z 

cutof f_dx_min 
cutof f_dy_min 
cutof f_dz_min 
cutof f_dx_max 
cutof f_dy_max 
cutof f_dz_max 

start_time_str 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min 

m_bin_diam_max 

decay_chain 


1 . 0 
2 . 0 
2 . 0 
2. 0 
2. 0 
2. 0 
2. 0 

"200 0APR01_0 00000 
"200 0APR01_0 00000 
" SF6 " 

6.035 

"200 0APR01_0 00000 

1.0 

"mass" 

0 . 0 
0 . 0 
1 . 0 
0 . 0 
0 . 0 

.false. 


NARAC MvM6 

lodi mvm06.nml (continued) 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 
dt_bin_out_str = 
samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 
dt_bin_out_str = 
samp_type 

/ 

&turb_param 
turb_param_h 
turb_param_z 
read_adapt_turb = 

/ 


'Bin 1: Integrated 
' air' 

'Source 1' 

" SF6" 

'xy' 

.true. 

10.0 
20.0 
"30:00" 

"30:00" 

'integrated' 


'Bin 2: Integrated 
' air' 

'Source 1' 

" SF6" 

'xy' 

.true. 

10.0 
20.0 
" 60 : 00 " 

"60:00" 

'integrated' 


'sigmav_simthry' 
'simthry' 

. true . 


Air concentration xy 


Air concentration xy 


&met_param 

/ 
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NARAC MvM6 

stnloc_mvm06.met (same as stnloc_mvm04.met) 
observ mvm06.met (same as observ mvm04.met) 
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HPAC MvM7 
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NARAC MvM7 


adapt_mvm07 .nml 

&adapt_control 
flag_map_adjust = .true. 
flag_debug = .true. 
opt_output_file = 'arac' 

/ 

&adapt_grid 

°pt_grid_file = 'gridgen' 

file_met_grid = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 

/ 


&adapt_metdata 
file_met_field 
opt_src_field 
opt_src_obs 
file_src_obs 
file_src_station 
flag_station_km 
nmethod 


/ 


'adapt_mtomO 7_2 0 0 0APRO1_0 0 0 0 0 0.nc' 
'none' 

'ascii2' 

'observ.met' 

' stnloc.met' 

.false. 

2 


&adapt_field2d 

zO = 0.008 

hgt_boundary_layer = 500.0 

inv_monin_obukhov_len = -0.02 

hgt_vert_coord = 'zAGL' 

/ 

hgt_blend_layer = 10.0 

hgt_surface_layer = 10.0 


&adapt_method 
obs_date_time 
flag_mc_adjust 
flag_upr_in_sl 
opt_met_type 
opt_method 
opt_dist_wt 
opt_int_blend 
opt_int_sfc 
opt_int_upr 
opt_int_horz 
opt_vert_grid 
op t_wind_hor z 
opt_wind_vert 
max_veer_vert 
blend_max_veer 
/ 


'2000APR01_000000' 
.false. 

.true. 

' wind2d' 

' sparse_data' 

' inv_horz_dist_sq' 
' sfc_difference ' 

' log' 

'linear' 

'distance_wt' 

'sigmaZ' 

' uv' 

'spddir' 

240.0 
240.0 


&adapt_method 

opt_method 

opt_met_type 

obs_date_time 

/ 


' turb_nml' 

'turb_none' 

"200 0APR01 000000" 


&adapt_turbulence 

turb_param_h = 'sigmav_simthry' 

turb_param_z = 'simthry' 

/ 


NARAC \Iv\I7 


lodi_files_mvm07. nml 

&grid_name 
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num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

'/ModelDevelopDiag/vandv/mtom/cgrid/conc 1 40 - 

18000 0 001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvmO 7/mgrid/adapt_mtom0 7_2 0 0 0APR01_0 0 0 0 0 0. 
nc' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 

NARAC MvM7 

lodi_mvm07 .nml 

&prob_setup 
title 

tstart_str 
tstop_str 
dt_dump_str 
nbins 
nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 

&thist_param 

/ 

&src_param 

source_id = 'Source 1' 

geom_time_strs = "2000APR01_000000" 

max_num_part = 500000 

source_model = 'neutral' 

src_agl_flg = .true. 

geom_type = 2 

mean_x = 0.0 

mean_y = 0.0 

mean_z = 2.0 

std_x = 1.0 

std_y = 1.0 

std_z = 1.0 

cutoff_dx_min = 2.0 
cutoff_dy_min = 2.0 
cutoff_dz_min = 2.0 
cutoff_dx_max = 2.0 
cutoff_dy_max = 2.0 
cutoff dz max = 2.0 


= 'Model vs. Model: Source 7' 
= "200 0APR01_0 00000" 

= "200 0APR01_010000" 

= " 30 : 0 " 

= 2 
= 1 

= 'nongauss' 

= 1 

= "200 0APR01_0 00000" 

= 'arac' 
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/ 


start_time_str 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min 

m_bin_diam_max 

decay_chain 


"200 0APR01_0 00000" 
"200 OAPRO1_0 00000" 
" SF6 " 

6.035 

"200 OAPRO1_0 00000" 
1.0 

"mass" 

0.0 
0.0 
1.0 
0.0 
0.0 

.false. 


NARAC MvM7 

lodi_mvm07.nml (continued) 

&bin_param 


bin id 

= 

'Bin 1: Integrated 

type 

= 

'air' 

source_list 

= 

'Source 1' 

species_name 

= 

" SF6" 

orientation 

= 

'xy' 

bin agl fig 

= 

.true. 

position 

= 

10.0 

width 

= 

20.0 

dt_samp_str 

= 

"30 : 00" 

dt_bin_out_str 

= 

"30 : 00" 

samp type 

= 

'integrated' 

• in param 

bin id 

= 

'Bin 2: Integrated 

type 

= 

' air' 

source list 

= 

'Source 1' 

species_name 

= 

" SF6" 

orientation 

= 

'xy' 

bin agl fig 

= 

.true. 

position 

= 

10.0 

width 

= 

20.0 

dt samp str 

= 

o 

o 

o 

kO 

dt_bin_out_str 

= 

o 

o 

o 

KD 

samp_type 

= 

'integrated' 

urb_param 

turb param h 

= 

'sigmav simthry' 

turb param z 

= 

'simthry' 

read_adapt_turb 

= 

.true. 


&met_param 

/ 


NARAC MvM7 


stnloc_mvm07.met 

SFC 

'MtoMl' 0.0 0.0 10.0 

UPR 
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MtoMl 


0.0 0.0 


NARAC MvM7 

observ mvm07.met 


METDATASET '2000APR01_000000' 
UPR 


MtoMl' 

0.5 

90 

3.0643 

MtoMl' 

1 

90 

3.5578 

MtoMl' 

1.5 

90 

3.838 

MtoMl' 

2 

90 

4.0316 

MtoMl' 

2.5 

90 

4.1783 

MtoMl' 

3 

90 

4.2956 

MtoMl' 

3.5 

90 

4.3929 

MtoMl' 

4 

90 

4.4756 

MtoMl' 

4.5 

90 

4.5473 

MtoMl' 

5 

90 

4.6104 

MtoMl' 

5.5 

90 

4.6666 

MtoMl' 

6 

90 

4.7173 

MtoMl' 

6.5 

90 

4.7632 

MtoMl' 

7 

90 

4.8052 

MtoMl' 

7.5 

90 

4.8439 

MtoMl' 

8 

90 

4.8796 

MtoMl' 

8.5 

90 

4.9128 

MtoMl' 

9 

90 

4.9438 

MtoMl' 

9.5 

90 

4.9728 

MtoMl' 

10 

90 

5 

MtoMl' 

10.5 

90 

5.0257 

MtoMl' 

11 

90 

5.0499 

MtoMl' 

11.5 

90 

5.0729 

MtoMl' 

12 

90 

5.0947 

MtoMl' 

12.5 

90 

5.1155 

MtoMl' 

13 

90 

5.1353 

MtoMl' 

13.5 

90 

5.1541 

MtoMl' 

14 

90 

5.1722 

MtoMl' 

14.5 

90 

5.1895 

MtoMl' 

15 

90 

5.2061 

MtoMl' 

15.5 

90 

5.2221 

MtoMl' 

16 

90 

5.2374 

MtoMl' 

16.5 

90 

5.2522 

MtoMl' 

17 

90 

5.2664 

MtoMl' 

17.5 

90 

5.2802 

MtoMl' 

18 

90 

5.2935 

MtoMl' 

18.5 

90 

5.3063 

MtoMl' 

19 

90 

5.3187 

MtoMl' 

19.5 

90 

5.3308 

MtoMl' 

20 

90 

5.3424 

MtoMl' 

20.5 

90 

5.3538 

MtoMl' 

21 

90 

5.3647 

MtoMl' 

21.5 

90 

5.3754 

MtoMl' 

22 

90 

5.3858 

MtoMl' 

22.5 

90 

5.3959 

MtoMl' 

23 

90 

5.4057 

MtoMl' 

23.5 

90 

5.4153 


NARAC MvM7 

observ_mvm07.met (continued) 

'MtoMl' 24 90 5.4246 

'MtoMl' 24.5 90 5.4337 
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' MtoMl' 
' MtoMl' 
' MtoMl' 
' MtoMl' 
' MtoMl' 
' MtoMl' 
7 MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 
'MtoMl' 


25 90 

25.5 90 

26 90 

26.5 90 

27 90 

27.5 90 

28 90 

28.5 90 

29 90 

29.5 90 

30 90 

30.5 90 

31 90 

31.5 90 

32 90 

32.5 90 

33 90 

33.5 90 

34 90 

34.5 90 

35 90 

35.5 90 

36 90 

36.5 90 

37 90 

37.5 90 

38 90 

38.5 90 

39 90 

39.5 90 

40 90 

40.5 90 

41 90 

41.5 90 

42 90 

42.5 90 

43 90 

43.5 90 

44 90 

44.5 90 

45 90 

45.5 90 

46 90 

46.5 90 

47 90 

47.5 90 

48 90 

48.5 90 

49 90 

49.5 90 

50 90 


5.4425 
5.4512 
5.4596 
5.4679 
5.476 
5.4838 
5.4915 
5.4991 
5.5065 
5.5137 
5.5208 
5.5277 
5.5345 
5.5412 
5.5477 
5.5541 
5.5604 
5.5666 
5.5726 
5.5786 
5.5845 
5.5902 
5.5959 
5.6014 
5.6069 
5.6122 
5.6175 
5.6227 
5.6278 
5.6329 
5.6378 
5.6427 
5.6475 
5.6523 
5.6569 
5.6615 
5.6661 
5.6705 
5.6749 
5.6793 
5.6836 
5.6878 
5.692 
5.6961 
5.7002 
5.7042 
5.7081 
5.712 
5.7159 
5.7197 
5.7234 
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HPAC MvM8 


HPAC SCREEN 

agent 

SF6 

New Project Setup 

project file name 

mvm8.pii 


coordinates 

cartesian 


local oriqinx 

0 km 


local oriqiny 

0 km 


local origin latitude 

45 deqN 


local oriqin longitude 

Odeq 


reference times 

UTC 


local time ofOO:OOZ 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off 

New Project Editor 

edit mode 

advanced 

Material Editor 

type 

qas 


name 

SF6 


material file 

Sf6.mtl 


units 

kg 


bin boundaries 

NA 


output fiaqs 

surf, dosaqe, surf, depos 


daytime decay rate 

0 


niqhttime decay rate 

0 

Gas Parameter Editor 

density ratio (rho/rho air) 

5.02917 


deposition veloc ity 

0 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


heiqht of release 

250 m 


source uncertainly 

no 


specification 

simple 


randomize location 

no 


release duration 

instantaneous 


release rate 

NA 


aqent mass 

lkq 


xsize 

1 m 


ysize 

1 m 


zsize 

lm 


momentum 

0 


buoyancy 

0 

lime Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

100 


maximum time step 

60 sec 


output interval 

15 min 

Domain Editor 

xmin domain 

-50 km 


xmax domain 

1 km 


vmin domain 

-5 km 


ymax domain 

5 km 


horizontal resolution 

default 


vertical domain max heiqht 

500 m 


vertical domain resolution 

10 m 

Weather Editor 

weather data type 

surface obs only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability lenqth scale 

NA 


large scale variability variance 

NA 


surface touqhness 

0.008 m 


prec ipitation 

none 


obs. time bin size 

1 hr 

Meteoroloqv Option Editor 

save meteoroloqy fields 

no 


terrain/land cover file 

no 

Surface Observations File 

file name 

mvmU.sfc 


elevation 

10 m 


zi (boundary layer heiqht) 

500 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Lenqth) 

-50 m 

Options Editor 

puff split q rid level 

2 


surface resolution 

default 


puff q rid resolution 

0 


boundary layer pts 

11 


stable atm. turbulence 

le-2 rrT2/s~2 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 rrT2/s~2 


calm cond. scale 

1000 m 


surface dosaqe heiqht 

0 


minimum puff mass 

le-20 


conditional averaqinq time 

1 hr 1 


lumped boundary layer 

no 
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NARAC MvM8 


adapt_mvm08.nml (same as adapt_mvm07.nml) 

lodi_files_mvm08.nml 

&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/conc_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvm07/mgrid/adapt_mtom07_2000APR01_000000. 
nc' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 

NARAC MvM8 

lodi_mvm08.nml 

&prob_setup 
title 

tstart_str 
tstop_str 
dt_dump_str 
nbins 
nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 

&thist_param 

/ 

Sesrc_param 

source_id = 'Source 1' 

geom_time_strs = "2000APR01_000000" 

max_num_part = 500000 

source_model = 'neutral' 

src_agl_flg = .true. 

geom_type = 2 

mean_x = 0.0 

mean_y = 0.0 

mean_z = 250.0 

std_x = 1.0 

std_y = 1.0 

std z = 1.0 


= 'Model vs. Model: Source 8' 
= "200 0APR01_0 00000" 

= "200 0APR01_010000" 

= "30:0" 

= 2 
= 1 

= 'nongauss' 

= 1 

= "200 0APR01_0 00000" 

= 'arac' 
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cutof f_dx_min 
cutof f_dy_min 
cutof f_dz_min 
cutof f_dx_max 
cutof f_dy_max 
cutof f_dz_max 

start_time_str 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min 

m_bin_diam_max 

decay_chain 

/ 


2.0 
2.0 
2.0 
2.0 
2.0 
2.0 

"200 0APR01_0 00000 
"200 0APR01_0 00000 
" SF6 " 

6.035 

"200 0APR01_0 00000 
1.0 

"mass" 

0.0 
0.0 
1.0 
0.0 
0.0 

.false. 


NARAC MvM8 

lodi_mvm 08 .nml (continued) 

&bin_param 


bin id 

= 

'Bin 1: Integrated 

type 

= 

' air' 

source_list 

= 

'Source 1' 

species name 

= 

" SF6" 

orientation 

= 

'xy' 

bin agl fig 

= 

.true. 

position 

= 

10.0 

width 

= 

20.0 

dt_samp_str 

= 

"30 : 00" 

dt bin out str 

= 

"30 : 00" 

samp_type 

= 

'integrated' 

• in param 

bin id 

= 

'Bin 2: Integrated 

type 

= 

'air' 

source_list 

= 

'Source 1' 

species_name 

= 

"SF6" 

orientation 

= 

'xy' 

bin agl fig 

= 

.true. 

position 

= 

10.0 

width 

= 

20.0 

dt_samp_str 

= 

o 

o 

o 

KD 

dt_bin_out_str 

= 

O 

O 

O 

KD 

samp type 

= 

'integrated' 

urb param 

turb param h 

= 

'sigmav simthry' 

turb param z 

= 

'simthry' 

read adapt turb 

= 

.true. 


&met_param 

/ 
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NARAC MvM8 

stnloc_mvm08.met (same as stnloc_mvm07.met) 
observ mvm08.met (same as observ mvm07.met) 
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HPAC MvM9 


mvm9 

HPAC 3.2 parameter 

HPAC 3.2 value 

HPAC SCREEN 

agent 

SF6 

New Project Setup 

project file name 

mvm9.pii 


coordinates 

cartesian 


local oriqinx 

0 km 


local origin y 

0 km 


local origin latitude 

45 deqN 


local origin longitude 

Odeq 


reference times 

UTC 


local time of 00:00Z 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off 

New Roject Editor 

edit mode 

advanced 

Material Editor 

type 

qas 


name 

SF6 


material file 

Sf6.mtl 


units 

kg 


bin boundaries 

NA 


output flags 

surf, dosaqe, surf, depos. 


daytime decay rate 

0 


nighttime decay rate 

0 

Gas Parameter Editor 

density ratio (rho/riio air) 

5.02917 


deposition veloc ity 

0 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


height of release 

400 m 


source uncertainty 

no 


spec ification 

simple 


randomize location 

no 


release duration 

instantaneous 


release rate 

NA 


agent mass 

1 kq 


xsize 

1 m 


y size 

lm 


zsize 

1 m 


momentum 

0 


buoyancy 

0 

lime Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

100 


maximum time step 

60 sec 


output interval 

15 min 

Domain Editor 

xmin domain 

-50 km 


xmax domain 

1 km 


ymin domain 

-5 km 


ymax domain 

5 km 


horizontal resolution 

default 


vertical domain max height 

500 m 


vertical domain resolution 

10 m 

Weather Editor 

weatherdata type 

surface obs. only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability length scale 

NA 


large scale variability variance 

NA 


surface toughness 

0.008 m 


precipitation 

none 


obs. time bin size 

1 hr 

Meteorology Option Editor 

save meteorology fields 

no 


terrain/land cover file 

no 

Surface Observations File 

file name 

mvmU.sfc 


elevation 

10 m 


zi (boundary layer height) 

500 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Length) 

-50 m 

Options Editor 

puff split grid level 

2 


surface resolution 

default 


puff qrid resolution 

0 


boundary layer pts 

11 


stable atm. turbulence 

le-2 m~2/s~2 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 m~2/s~2 


calm cond. scale 

1000 m 


surface dosage height 

0 


minimum puff mass 

le-20 


conditional averaging time 

1 hr 


lumped boundary layer 

no 
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NARAC MvM9 


adapt_mvm09.nml (same as adapt_mvm07.nml) 

lodi_files_mvm09.nml 

&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/conc_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvm07/mgrid/adapt_mtom07_2000APR01_000000. 
nc' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 

NARAC MvM9 

lodi_mvm09.nml 

&prob_setup 
title 

tstart_str 
tstop_str 
dt_dump_str 
nbins 
nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 

&thist_param 

/ 

Sesrc_param 

source_id = 'Source 1' 

geom_time_strs = "2000APR01_000000" 

max_num_part = 500000 

source_model = 'neutral' 

src_agl_flg = .true. 

geom_type = 2 

mean_x = 0.0 

mean_y = 0.0 

mean_z = 400.0 

std_x = 1.0 

std_y = 1.0 

std z = 1.0 


= 'Model vs. Model: Source 9' 
= "200 0APR01_0 00000" 

= "200 0APR01_010000" 

= "30:0" 

= 2 
= 1 

= 'nongauss' 

= 1 

= "200 0APR01_0 00000" 

= 'arac' 
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cutof f_dx_min 
cutof f_dy_min 
cutof f_dz_min 
cutof f_dx_max 
cutof f_dy_max 
cutof f_dz_max 

start_time_str 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min 

m_bin_diam_max 

decay_chain 

/ 


2.0 
2.0 
2.0 
2.0 
2.0 
2.0 

"200 0APR01_0 00000 
"200 0APR01_0 00000 
" SF6 " 

6.035 

"200 0APR01_0 00000 
1.0 

"mass" 

0.0 
0.0 
1.0 
0.0 
0.0 

.false. 


NARAC MvM9 

lodi_mvm09.nml (continued) 

&bin_param 


bin id 

= 

'Bin 1: Integrated 

type 

= 

' air' 

source_list 

= 

'Source 1' 

species name 

= 

" SF6" 

orientation 

= 

'xy' 

bin agl fig 

= 

.true. 

position 

= 

10.0 

width 

= 

20.0 

dt_samp_str 

= 

"30 : 00" 

dt bin out str 

= 

"30 : 00" 

samp_type 

= 

'integrated' 

• in param 

bin id 

= 

'Bin 2: Integrated 

type 

= 

'air' 

source_list 

= 

'Source 1' 

species_name 

= 

"SF6" 

orientation 

= 

'xy' 

bin agl fig 

= 

.true. 

position 

= 

10.0 

width 

= 

20.0 

dt_samp_str 

= 

o 

o 

o 

KD 

dt_bin_out_str 

= 

O 

O 

O 

KD 

samp type 

= 

'integrated' 

urb param 

turb param h 

= 

'sigmav simthry' 

turb param z 

= 

'simthry' 

read adapt turb 

= 

.true. 


&met_param 

/ 
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NARAC MvM9 

stnloc_mvm09.met (same as stnloc_mvm07 .met) 
observ mvm09.met (same as observ mvm07.met) 
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HPAC MvMlO 


mvmlO 

HPAC 3.2 parameter 

HPAC 3.2 value 

HPAC SCREEN 

agent 

5 micron particles 

New Project Setup 

project file name 

mvmlO.pri 


coordinates 

cartesian 


local oriqinx 

0 km 


local origin y 

0 km 


local origin latitude 

45 deqN 


local origin longitude 

Odeq 


reference times 

UTC 


local time of 00:00Z 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off 

New Roject Editor 

edit mode 

advanced 

Material Editor 

type 

partic le 


name 

STUFF5 


material file name 

Stuff5.mtl 


units 

kg 


bin boundaries 

4.89 n , 5.12 n 


output flags 

surf, dosaqe, surf, depos. 


daytime decay rate 

0 


nighttime decay rate 

0 

Particle Parameter Editor 

density 

500 kq/m~3 


deposition veloc ity 

NA 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


height of release 

2 m 


source uncertainty 

no 


spec ification 

simple 


randomize location 

NA 


release duration 

15 min 


release rate 

0.067 kq/min 


agent mass 

1 kq 


xsize 

NA 


y size 

lm 


zsize 

1 m 


momentum 

0 


buoyancy 

0 

Time Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

100 


maximum time step 

60 sec 


output interval 

15 min 

Domain Editor 

xmin domain 

-50 km 


xmax domain 

1 km 


ymin domain 

-5 km 


ymax domain 

5 km 


horizontal resolution 

default 


vertical domain max height 

500 m 


vertical domain resolution 

10 m 

Weather Editor 

weatherdata type 

surface obs. only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability length scale 

NA 


large scale variability variance 

NA 


surface toughness 

0.008 m 


precipitation 

none 


obs. time bin size 

1 hr 

Meteorology Option Editor 

save meteorology fields 

no 


terrain/land cover file 

no 

Surface Observations File 

file name 

mvmN.sfc 


elevation 

10 m 


zi (boundary layer height) 

500 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90 deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Length) 

-500 m 

Options Editor 

puff split grid level 

2 


surface resolution 

default 


puff qrid resolution 

0 


boundary layer pts 

11 


stable atm. turbulence 

le-2 m~2/s~2 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 m~2/s~2 


calm cond. scale 

1000 m 


surface dosage height 

0 


minimum puff mass 

le-20 


conditional averaging time 

1 hr 


lumped boundary layer 

no 
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NARAC MvMlO 


adapt_mvmlO .nml (same as adapt_mvm07 .nml) 

lodi_filesmvmlO.nml 

&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/conc_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvm04/mgrid/adapt_mtom4_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 

NARAC MvMlO 

lodi_mvml0.nml 

&prob_setup 
title 

tstart_str 
tstop_str 
dt_dump_str 
nbins 
nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 

&thist_param 

/ 

Sesrc_param 

source_id = 'Source 1' 

geom_time_strs = "200 0APR01_0 00000" 

max_num_part = 500000 

source_model = 'neutral' 

src_agl_flg = .true. 

geom_type = 2 

mean_x = 0.0 

mean_y = 0.0 

mean_z = 2.0 

std_x = 1.0 

std_y = 1.0 


= 'Model vs. Model: Source 10' 
= "200 0APR01_0 00000" 

= "200 0APR01_010000" 

= "30:0" 

= 3 
= 1 

= 'nongauss' 

= 1 

= "200 0APR01_0 00000" 

= 'arac' 
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std_z 

cutoff_dx_min 
cutoff_dy_min 
cutoff_dz_min 
cutoff_dx_max 
cutoff_dy_max 
cutoff_dz_max 

start_time_str = 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min = 

m_bin_diam_max = 

decay_chain 


1.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 

"200 0APR01_0 00000" 
"200 0APR01_0 01500" 
"Stuff5" 

500.0 

"200 0APR01_0 00000" 
1.lllle-3 
"mass" 

0.003 
0.0 
1.0 

5 . Oe-6 
5 . Oe-6 
.false. 


NARAC MvMlO 

lodi mvmlO.nml (continued) 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

position 

dt_bin_out_str 

samp_type 


'Bin 1: Integrated Air concentration xy 
' air' 

'Source 1' 

"Stuff5" 

'xy' 

.true. 

10.0 
20.0 
"30:00" 

"30:00" 

'integrated' 


'Bin 2: Integrated Air concentration xy 
' air' 

'Source 1' 

"Stuff5" 

'xy' 

.true. 

10.0 
20.0 
"60:00" 

"60:00" 

'integrated' 


'Bin 3: Instantaneous deposition' 
'dep_dry' 

'Source 1' 

"Stuff5" 

'xy' 

0.0 

"30 : 00" 

' instantaneous' 
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' sigmav_simthry 
' simthry , 

. true . 


/ 

&turb_param 
turb_param_h 
turb_param_z 
read_adapt_turb = 


&met_param 

/ 

NARAC MvMlO 

stnlocmvmlO.met (same as stnloc_mvm07 .met) 
observ mvmlO.met (same as observ mvm07.met) 
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HPAC MvMll 


mvmll 

HPAC 3.2 parameter 

HPAC 3.2 value 

HPAC SCREEN 

agent 

5 micron particles 

New Project Setup 

project file name 

mvmll.pii 


coordinates 

cartesian 


local oriqinx 

0 km 


local origin y 

0 km 


local origin latitude 

45 deqN 


local origin longitude 

Odeq 


reference times 

UTC 


local time of 00:00Z 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off 

New Roject Editor 

edit mode 

advanced 

Material Editor 

type 

partic le 


name 

STUFF5 


material file name 

Stuff5.mtl 


units 

kg 


bin boundaries 

4.89 n , 5.12 n 


output flags 

surf, dosaqe, surf, depos. 


daytime decay rate 

0 


nighttime decay rate 

0 

Particle Parameter Editor 

density 

500 kq/m~3 


deposition veloc ity 

NA 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


height of release 

250 m 


source uncertainty 

no 


spec ification 

simple 


randomize location 

NA 


release duration 

15 min 


release rate 

0.067 kq/min 


agent mass 

1 kq 


xsize 

NA 


y size 

lm 


zsize 

1 m 


momentum 

0 


buoyancy 

0 

Time Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

100 


maximum time step 

60 sec 


output interval 

15 min 

Domain Editor 

xmin domain 

-50 km 


xmax domain 

1 km 


ymin domain 

-5 km 


ymax domain 

5 km 


horizontal resolution 

default 


vertical domain max height 

500 m 


vertical domain resolution 

10 m 

Weather Editor 

weatherdata type 

surface obs. only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability length scale 

NA 


large scale variability variance 

NA 


surface toughness 

0.008 m 


precipitation 

none 


obs. time bin size 

1 hr 

Meteorology Option Editor 

save meteorology fields 

no 


terrain/land cover file 

no 

Surface Observations File 

file name 

mvmN.sfc 


elevation 

10 m 


zi (boundary layer height) 

500 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90 deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Length) 

-500 m 

Options Editor 

puff split grid level 

2 


surface resolution 

default 


puff qrid resolution 

0 


boundary layer pts 

11 


stable atm. turbulence 

le-2 m~2/s~2 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 m~2/s~2 


calm cond. scale 

1000 m 


surface dosage height 

0 


minimum puff mass 

le-20 


conditional averaging time 

1 hr 


lumped boundary layer 

no 
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NARAC MvMll 


adapt_mvmll.nml (same as adapt_mvm07.nml) 

lodi_filesmvmll.nml 

&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/conc_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvm04/mgrid/adapt_mtom4_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 

NARAC MvMll 

lodi_mvmll.nml 

&prob_setup 
title 

tstart_str 
tstop_str 
dt_dump_str 
nbins 
nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 

&thist_param 

/ 

Sesrc_param 

source_id = 'Source 1' 

geom_time_strs = "200 0APR01_0 00000" 

max_num_part = 500000 

source_model = 'neutral' 

src_agl_flg = .true. 

geom_type = 2 

mean_x = 0.0 

mean_y = 0.0 

mean_z = 250.0 

std_x = 1.0 

std_y = 1.0 


= 'Model vs. Model: Source 11' 
= "200 0APR01_0 00000" 

= "200 0APR01_010000" 

= "30:0" 

= 3 
= 1 

= 'nongauss' 

= 1 

= "200 0APR01_0 00000" 

= 'arac' 
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std_z 

cutoff_dx_min 
cutoff_dy_min 
cutoff_dz_min 
cutoff_dx_max 
cutoff_dy_max 
cutoff_dz_max 

start_time_str = 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min = 

m_bin_diam_max = 

decay_chain 


1.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 

"200 0APR01_0 00000" 
"200 0APR01_0 01500" 
"Stuff5" 

500.0 

"200 0APR01_0 00000" 
1.lllle-3 
"mass" 

0.003 
0.0 
1.0 

5 . Oe-6 
5 . Oe-6 
.false. 


NARAC MvMll 

lodi mvmll.nml (continued) 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

position 

dt_bin_out_str 

samp_type 


'Bin 1: Integrated Air concentration xy 
' air' 

'Source 1' 

"Stuff5" 

'xy' 

.true. 

10.0 
20.0 
"30:00" 

"30:00" 

'integrated' 


'Bin 2: Integrated Air concentration xy 
' air' 

'Source 1' 

"Stuff5" 

'xy' 

.true. 

10.0 
20.0 
"60:00" 

"60:00" 

'integrated' 


'Bin 3: Instantaneous deposition' 
'dep_dry' 

'Source 1' 

"Stuff5" 

'xy' 

0.0 

"30 : 00" 

' instantaneous' 
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/ 

&turb_param 

turb_param_h = ' sigmav_simthry 7 
turb_param_z = 'simthry' 
read_adapt_turb = .true. 


&met_param 

/ 


NARAC MvMll 

stnlocmvmll.met (same as stnloc_mvm07.met) 
observ mvmll.met (same as observ mvm07.met) 
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HPAC MvM12 


mvml2 

HPAC 3.2 parameter 

HPAC 3.2 value 

HPAC SCREEN 

agent 

5 micron particles 

New Project Setup 

project file name 

mvml2.pri 


coordinates 

cartesian 


local oriqinx 

0 km 


local origin y 

0 km 


local origin latitude 

45 deqN 


local origin longitude 

Odeq 


reference times 

UTC 


local time of 00:00Z 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off 

New Roject Editor 

edit mode 

advanced 

Material Editor 

type 

partic le 


name 

STUFF5 


material file name 

Stuff5.mtl 


units 

kg 


bin boundaries 

4.89 n , 5.12 n 


output flags 

surf, dosaqe, surf, depos. 


daytime decay rate 

0 


nighttime decay rate 

0 

Particle Parameter Editor 

density 

500 kq/m~3 


deposition veloc ity 

NA 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


height of release 

400 m 


source uncertainty 

no 


spec ification 

simple 


randomize location 

NA 


release duration 

15 min 


release rate 

0.067 kq/min 


agent mass 

1 kq 


xsize 

NA 


y size 

lm 


zsize 

1 m 


momentum 

0 


buoyancy 

0 

Time Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

100 


maximum time step 

60 sec 


output interval 

15 min 

Domain Editor 

xmin domain 

-50 km 


xmax domain 

1 km 


ymin domain 

-5 km 


ymax domain 

5 km 


horizontal resolution 

default 


vertical domain max height 

500 m 


vertical domain resolution 

10 m 

Weather Editor 

weatherdata type 

surface obs. only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability length scale 

NA 


large scale variability variance 

NA 


surface toughness 

0.008 m 


precipitation 

none 


obs. time bin size 

1 hr 

Meteorology Option Editor 

save meteorology fields 

no 


terrain/land cover file 

no 

Surface Observations File 

file name 

mvmN.sfc 


elevation 

10 m 


zi (boundary layer height) 

500 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90 deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Length) 

-500 m 

Options Editor 

puff split grid level 

2 


surface resolution 

default 


puff qrid resolution 

0 


boundary layer pts 

11 


stable atm. turbulence 

le-2 m~2/s~2 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 m~2/s~2 


calm cond. scale 

1000 m 


surface dosage height 

0 


minimum puff mass 

le-20 


conditional averaging time 

1 hr 


lumped boundary layer 

no 
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NARAC MvM12 


adapt_mvml2 .nml (same as adapt_mvm07.nml) 
lodi files mvml2.nml 


&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/cone_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvm04/mgrid/adapt_mtom4_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 


NARAC \Iv\I12 

lodi mvml2 . nml 


&prob_setup 

title 

tstart_str 

tstop_str 

dt_dump_str 

nbins 

nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met format 


/ 


'Model vs. Model: Source 12' 
"200 0APR01_0 00000" 

"200 0APR01_010000" 

"30:0" 

3 

1 

'nongauss' 

1 

"200 0APR01_0 00000" 

'arac' 


&thist_param 

/ 


Sesrc_param 
source_id 
geom_time_strs 
max_num_p art 
source_model 
src_agl_flg 

geom_type 
mean_x 
mean_y 
mean_z 
std x 


= 'Source 1' 

= "200 0APR01_0 00000" 
= 500000 
= 'neutral' 

= .true. 

= 2 
= 0.0 
= 0.0 
= 400.0 
= 1.0 
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std_y 

std_z 

cutoff_dx_min 
cutoff_dy_min 
cutoff_dz_min 
cutoff_dx_max 
cutoff_dy_max 
cutoff_dz_max 

start_time_str = 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min = 

m_bin_diam_max = 

decay_chain 


1.0 
1.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 

"200 0APR01_0 00000" 
"2000APR01_001500" 
"Stuff5" 

500.0 

"200 0APR01_0 00000" 
1.lllle-3 
"mass" 

0.003 
0.0 
1.0 

5 . Oe-6 
5 . Oe-6 
.false. 


NARAC MvM12 

lodi mvml2.nml (continued) 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 
dt_bin_out_str = 
samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 
dt_bin_out_str = 
samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

position 

dt bin out str = 


'Bin 1: Integrated Air 
'air' 

'Source 1' 

"Stuff5" 

'xy' 

.true. 

10.0 
20.0 
"30:00" 

"30:00" 

'integrated' 


'Bin 2: Integrated Ai 
' air' 

'Source 1' 

"Stuff5" 

'xy' 

.true. 

10.0 
20.0 
"60:00" 

"60:00" 

'integrated' 


'Bin 3: Instantaneous 
'dep_dry' 

'Source 1' 

"Stuff5" 

'xy' 

0.0 

"30 : 00" 


concentration xy 


concentration xy 


deposition' 
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samp_type = 'instantaneous' 


&turb_param 

turb_param_h = 'sigmav_simthry' 
turb_param_z = 'simthry' 
read_adapt_turb = .true. 


&met_param 

/ 


NARAC MvM12 

stnloc_mvml2.met (same as stnloc_mvm07.met) 
observ mvm!2.met (same as observ mvm07.met) 
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HPAC MvM13 


mvml3 

HPAC 3.2 parameter 

HPAC 3.2 value 

HPAC SCREEN 

agent 

50 micron particles 

New Project Setup 

project file name 

mvml3.pri 


coordinates 

cartesian 


local oriqinx 

0 km 


local origin y 

0 km 


local origin latitude 

45 deqN 


local origin longitude 

Odeq 


reference times 

UTC 


local time of 00:00Z 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off 

New Roject Editor 

edit mode 

advanced 

Material Editor 

type 

partic le 


name 

STUFF50 


material file name 

Stuff50.mtl 


units 

kg 


bin boundaries 

48.9 n , 51.1 n 


output flags 

surf, dosaqe, surf, depos. 


daytime decay rate 

0 


nighttime decay rate 

0 

Particle Parameter Editor 

density 

500 kq/m~3 


deposition veloc ity 

NA 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


height of release 

2 m 


source uncertainty 

no 


spec ification 

simple 


randomize location 

no 


release duration 

instantaneous 


release rate 

NA 


agent mass 

1 kq 


xsize 

1 m 


y size 

lm 


zsize 

1 m 


momentum 

0 


buoyancy 

0 

lime Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

100 


maximum time step 

60 sec 


output interval 

15 min 

Domain Editor 

xmin domain 

-50 km 


xmax domain 

1 km 


ymin domain 

-5 km 


ymax domain 

5 km 


horizontal resolution 

default 


vertical domain max height 

500 m 


vertical domain resolution 

10 m 

Weather Editor 

weatherdata type 

surface obs. only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability length scale 

NA 


large scale variability variance 

NA 


surface toughness 

0.008 m 


precipitation 

none 


obs. time bin size 

1 hr 

Meteorology Option Editor 

save meteorology fields 

no 


terrain/land cover file 

no 

Surface Observations File 

file name 

mvmN.sfc 


elevation 

10 m 


zi (boundary layer height) 

500 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Length) 

-500 m 

Options Editor 

puff split grid level 

2 


surface resolution 

default 


puff qrid resolution 

0 


boundary layer pts 

11 


stable atm. turbulence 

le-2 m~2/s~2 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 m~2/s~2 


calm cond. scale 

1000 m 


surface dosage height 

0 


minimum puff mass 

le-20 


conditional averaging time 

1 hr 


lumped boundary layer 

no 
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NARAC MvM13 


adapt_mvml3.nml (same as adapt_mvm07.nml) 
lodi files mvml3.nml 


&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/cone_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvm04/mgrid/adapt_mtom4_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 


NARAC \Iv\I13 

lodi mvml3 . nml 


&prob_setup 

title 

tstart_str 

tstop_str 

dt_dump_str 

nbins 

nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met format 


/ 


'Model vs. Model: Source 13' 
"200 0APR01_0 00000" 

"200 0APR01_010000" 

"30:0" 

3 

1 

'nongauss' 

1 

"200 0APR01_0 00000" 

'arac' 


&thist_param 

/ 


Sesrc_param 
source_id 
geom_time_strs 
max_num_p art 
source_model 
src_agl_flg 

geom_type 
mean_x 
mean_y 
mean_z 
std x 


= 'Source 1' 

= "200 0APR01_0 00000" 
= 500000 
= 'neutral' 

= .true. 

= 2 
= 0.0 
= 0.0 
= 2.0 
= 1.0 
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std_y 

std_z 

cutoff_dx_min 
cutoff_dy_min 
cutoff_dz_min 
cutoff_dx_max 
cutoff_dy_max 
cutoff_dz_max 

start_time_str = 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min = 

m_bin_diam_max = 

decay_chain 


NARAC MvM13 


1 . o 
1 . o 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 

"200 0APR01_0 00000" 
"200 0APR01_0 00000" 
"Stuff50" 

500.0 

"200 0APR01_0 00000" 
1.0 

"mass" 

0.003 
0.0 
1.0 

50 . Oe-6 
50 . Oe-6 
.false. 


lodi mvml3.nml (continued) 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 
dt_bin_out_str = 
samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 
dt_bin_out_str = 
samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

position 

dt bin out str = 


'Bin 1: Integrated Air 
'air' 

'Source 1' 

"Stuff50" 

'xy' 

.true. 

10.0 
20.0 
"30:00" 

"30:00" 

'integrated' 


'Bin 2: Integrated Ai 
' air' 

'Source 1' 

"Stuff50" 

'xy' 

.true. 

10.0 
20.0 
"60:00" 

"60:00" 

'integrated' 


'Bin 3: Instantaneous 
'dep_dry' 

'Source 1' 

"Stuff50" 

'xy' 

0.0 

"30 : 00" 


concentration xy 


concentration xy 


deposition' 
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samp_type = 'instantaneous' 


&turb_param 

turb_param_h = 'sigmav_simthry' 
turb_param_z = 'simthry' 
read_adapt_turb = .true. 


&met_param 

/ 


NARAC MvM13 

stnloc_mvml3.met (same as stnloc_mvm07.met) 
observ mvml3.met (same as observ mvm07.met) 
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HPAC MvM14 


mvml4 

HPAC 3.2 parameter 

HPAC 3.2 value 

HPAC SCREEN 

agent 

50 micron particles 

New Project Setup 

project file name 

mvml4.pri 


coordinates 

cartesian 


local oriqinx 

0 km 


local oriqiny 

0 km 


local origin latitude 

45 deqN 


local oriqin lonqitude 

Odeq 


reference times 

UTC 


local time of 00:00Z 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off i 

New Roject Editor 

edit mode 

advanced 

Material Editor 

type 

partic le 


name 

STUFF50 


material file name 

Stuff50.mtl 


units 

kg 


bin boundaries 

48.9 p, , 51.1 p. 


output flaqs 

surf, dosaqe, surf, depos 


daytime decay rate 

0 


nighttime decay rate 

0 

Particle Parameter Editor 

density 

500kq/m~3 


deposition veloc ity 

NA 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


height of release 

2 m 


source uncertainty 

no 


spec ification 

simple 


randomize location 

NA 


release duration 

15 min 


release rate 

0.067 kq/min 


agent mass 

1 kq 


xsize 

NA 


y size 

1 m 


zsize 

1 m 


momentum 

0 


buoyancy 

0 

lime Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

100 


maximum time step 

60 sec 


output interval 

15 min 

Domain Editor 

xmin domain 

-50 km 


xmax domain 

1 km 


ymin domain 

-5 km 


ymax domain 

5 km 


horizontal resolution 

default 


vertical domain max heiqht 

500 m 


vertical domain resolution 

10 m 

Weather Editor 

weather data type 

surface obs only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability length scale 

NA 


large scale variability variance 

NA 


surface roughness 

0.008 m 


precipitation 

none 


obs. time bin size 

1 hr 

Meteorology Option Editor 

save meteoroloqy fields 

no 


terrain/land cover file 

no 

Surface Observations Rle 

file name 

mvmN.sfc 


elevation 

10 m 


zi (boundary layer height) 

500 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90 deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Length) 

-500 m 

Options Editor 

puff split grid level 

2 


surface resolution 

default 


puff qrid resolution 

0 i 


boundary layer pts 

11 


stable atm. turbulence 

le-2 m'V/sTZ 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 m~2/s~2 


calm cond. scale 

1000 m 


surface dosaqe height 

0 


minimum puff mass 

le-20 


conditional averaging time 

1 hr 


lumped boundary layer 

no 
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NARAC MvM14 


adapt_mvml4 .nml (same as adapt_mvm07.nml) 
lodi files mvml4.nml 


&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/cone_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvm04/mgrid/adapt_mtom4_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 

NARAC \I v \I14 

lodi mvml4 . nml 


&prob_setup 

title 

tstart_str 

tstop_str 

dt_dump_str 

nbins 

nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 


'Model vs. Model: Source 14' 
"200 0APR01_0 00000" 

"200 0APR01_010000" 

"30:0" 

3 

1 

'nongauss' 

1 

"200 0APR01_0 00000" 

'arac' 


&thist_param 

/ 


Sesrc_param 
source_id 
geom_time_strs 
max_num_p art 
source_model 
src_agl_flg 

geom_type 
mean_x 
mean_y 
mean_z 
std x 


'Source 1' 

"200 0APR01_0 00000 
500000 
'neutral' 

.true. 

2 

0.0 
0.0 
2.0 
1.0 
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std_y 

std_z 

cutoff_dx_min 
cutoff_dy_min 
cutoff_dz_min 
cutoff_dx_max 
cutoff_dy_max 
cutoff_dz_max 

start_time_str = 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

m_bin_fract 

m_bin_diam_min = 

m_bin_diam_max = 

decay_chain 


NARAC MvM14 


1 . o 
1 . o 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 

"200 0APR01_0 00000" 
"2000APR01_001500" 
"Stuff50" 

500.0 

"200 0APR01_0 00000" 
1.lllle-3 
"mass" 

0.003 
0.0 
1.0 

50 . Oe-6 
50 . Oe-6 
.false. 


lodi mvml4.nml (continued) 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

position 

dt bin out str 


'Bin 1: Integrated Air concentration xy 
'air' 

'Source 1' 

"Stuff50" 

'xy' 

.true. 

10.0 
20.0 
"30:00" 

"30:00" 

'integrated' 


'Bin 2: Integrated Air concentration xy 
' air' 

'Source 1' 

"Stuff50" 

'xy' 

.true. 

10.0 
20.0 
"60:00" 

"60:00" 

'integrated' 


'Bin 3: Instantaneous deposition' 
'dep_dry' 

'Source 1' 

"Stuff50" 

'xy' 

0.0 

"30 : 00" 
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samp_type = 'instantaneous' 


&turb_param 

turb_param_h = 'sigmav_simthry' 
turb_param_z = 'simthry' 
read_adapt_turb = .true. 


&met_param 

/ 

NARAC \Iv\I 14 

stnloc_mvml4.met (same as stnloc_mvm07 .met) 
observ mvm!4.met (same as observ mvm07.met) 


B-63 



HPAC MvM15 


mvml5 

HPAC 3.2 parameter 

HPAC 3.2 value 

HPAC SCREEN 

agent 

50 micron particles 

New Project Setup 

project file name 

mvml5.pri 


coordinates 

cartesian 


local oriqinx 

0 km 


local origin y 

0 km 


local origin latitude 

45 deqN 


local origin longitude 

Odeq 


reference times 

UTC 


local time of 00:00Z 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off 

New Roject Editor 

edit mode 

advanced 

Material Editor 

type 

partic le 


name 

STUFF50 


material file name 

Stuff50.mtl 


units 

kg 


bin boundaries 

48.9 n , 51.1 n 


output flags 

surf, dosaqe, surf, depos. 


daytime decay rate 

0 


nighttime decay rate 

0 

Gas Parameter Editor 

density ratio (rho/riio air) 

500 kq/m~3 


deposition veloc ity 

NA 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


height of release 

750 m 


source uncertainty 

no 


spec ification 

simple 


randomize location 

NA 


release duration 

instantaneous 


release rate 

NA 


agent mass 

1 kq 


xsize 

1 m 


y size 

lm 


zsize 

1 m 


momentum 

0 


buoyancy 

0 

lime Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

400 


maximum time step 

60 sec 


output interval 

15 min 

Domain Editor 

xmin domain 

-100 


xmax domain 

1 km 


ymin domain 

-20 km 


ymax domain 

20 km 


horizontal resolution 

default 


vertical domain max height 

1000 m 


vertical domain resolution 

20 m 

Weather Editor 

weatherdata type 

surface obs only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability length scale 

NA 


large scale variability variance 

NA 


surface toughness 

0.008 m 


precipitation 

none 


obs. time bin size 

1 hr 

Meteorology Option Editor 

save meteorology fields 

no 


terrain/land cover file 

no 

Surface Observations File 

file name 

mvmN.sfc 


elevation 

10 m 


zi (boundary layer height) 

500 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90 deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Length) 

-500 m 

Options Editor 

puff split grid level 

2 


surface resolution 

default 


puff qrid resolution 

0 


boundary layer pts 

11 


stable atm. turbulence 

le-2 m~2/s~2 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 m~2/s~2 


calm cond. scale 

1000 m 


surface dosage height 

0 


minimum puff mass 

le-20 


conditional averaging time 

1 hr 


lumped boundary layer 

no 
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NARAC MvM15 


adapt_mvml5. nml 

&adapt_control 
flag_map_adjust = .true. 
flag_debug = .true. 
opt_output_file = 'arac' 

/ 

&adapt_grid 

°pt_grid_file = 'gridgen' 

file_met_grid = '/ModelDevelopDiag/vandv/mtom/mvml5/grid/mvml5_mgrid.nc' 

/ 


&adapt_metdata 
file_met_field 
opt_src_field 
opt_src_obs 
file_src_obs 
file_src_station 
flag_station_km 
nmethod 


/ 


'adapt_mtoml5_2 0 0 0APRO1_0 0 0 0 0 0.nc' 
'none' 

'ascii2' 

'observ.met' 

' stnloc.met' 

.false. 

2 


&adapt_field2d 

zO = 0.008 

hgt_boundary_layer = 500.0 

inv_monin_obukhov_len = -0.002 
hgt_vert_coord = 'zAGL' 

/ 

hgt_blend_layer = 10.0 

hgt_surface_layer = 10.0 


&adapt_method 
obs_date_time 
flag_mc_adjust 
flag_upr_in_sl 
opt_met_type 
opt_method 
opt_dist_wt 
opt_int_blend 
opt_int_sfc 
opt_int_upr 
opt_int_horz 
opt_vert_grid 
op t_wind_hor z 
opt_wind_vert 
max_veer_vert 
blend_max_veer 
/ 


'2000APR01_000000' 
.false. 

.true. 

' wind2d' 

' sparse_data' 

' inv_horz_dist_sq' 
' sfc_difference ' 

' log' 

'linear' 

'distance_wt' 

'sigmaZ' 

' uv' 

'spddir' 

240.0 
240.0 


&adapt_method 

opt_method 

opt_met_type 

obs_date_time 

/ 


' turb_nml' 

'turb_none' 

"200 0APR01 000000" 


&adapt_turbulence 

turb_param_h = 'sigmav_simthry' 

turb_param_z = 'simthry' 

/ 


NARAC \Iv\I15 


lodi_files_mvml5. nml 

&grid_name 
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1 


num_m_grids = 
m_grid_name = 

'/ModelDevelopDiag/vandv/mtom/mvml5/grid/mvml5_mgrid.nc' 
c_grid_name = 

'/ModelDevelopDiag/vandv/mtom/mvml5/cgrid/cone_1_100_- 

48000 0 001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvml5/mgrid/adapt_mtoml5_2 0 0 0APR01_0 0 0 0 0 0. 

nc' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 

NARAC \Iv\I15 


lodi mvml5. nml 

&prob_setup 

title 

tstart_str 

tstop_str 

dt_dump_str 

nbins 

nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 

&thist_param 

/ 

&src_param 


source id 

= 

'Source 1' 


geom_time_strs 

= 

"200 0APR01 

000000 

max num part 

= 

500000 


source_model 

= 

'neutral' 


sre agl fig 

= 

.true. 


geom type 

= 

2 


mean x 

= 

o 

o 


mean y 

= 

o 

o 


mean z 

= 

750.0 


std x 

= 

1.0 


std_y 

= 

1 . 0 


std z 

= 

1 . 0 


cutoff dx min 

= 

ISJ 

o 


cutoff dy min 

= 

2.0 


cutoff dz min 

= 

ISJ 

o 


cutoff dx max 

= 

2.0 


cutof f_dy_max 

= 

ISJ 

o 


cutoff dz max 

= 

2.0 


start time str 

- 

"200 0APR01 

000000 

stop_time_str 

= 

"200 0APR01 

^000000 

species 

= 

"Stuff50" 


density 

= 

500 . 0 



= 'Model vs. Model: Source 15' 
= "200 0APR01_0 00000" 

= "200 0APR01_04 0 000" 

= " 30 : 0 " 

= 6 
= 1 

= 'nongauss' 

= 1 

= "200 0APR01_0 00000" 

= 'arac' 
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er_time_strs 
emiss_rates 
er_units_type 
nset_dep_vel 
precip_coeff 
m_bin_fract 
m_bin_diam_min 
m_bin_diam_max 
decay_chain 

/ 

NARAC MvM15 

lodi_mvml5.nml 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 


= "200 0APR01_0 00000" 
= 1.0 
= "mass" 

= 0.003 
= 0.0 
= 1.0 
= 50.0e-6 
= 50.Oe-6 
= .false. 


(continued) 


= 'Bin 1: Integrated Air concentration xy' 
= 'air' 

= 'Source 1' 

= "Stuff50" 

= 'xy' 

= .true. 

= 10.0 
= 20.0 
= "30:00" 

= "30:00" 

= 'integrated' 


= 'Bin 2: Integrated Air concentration xy' 
= 'air' 

= 'Source 1' 

= "Stuff50" 

= ' xy' 

= .true. 

= 10.0 
= 20.0 
= "60:00" 
a "60:00" 

= 'integrated' 


= 'Bin 3: Integrated Air concentration xy' 
= 'air' 

= 'Source 1' 

= "Stuff50" 

= ' xy' 

= .true. 

= 10.0 
= 20.0 
= " 2 : 00 : 00 " 

= " 2 : 00 : 00 " 

= 'integrated' 


= 'Bin 4: Integrated Air concentration xy' 
= 'air' 

= 'Source 1' 

= "Stuff50" 

= ' xy' 

= .true. 

= 10.0 
= 20.0 
= "3:00:00" 

= "3:00:00" 

= 'integrated' 
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NARAC MvM15 

lodi mvml5.nml (continued) 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

position 

dt_bin_out_str 

samp_type 

/ 

&turb_param 

turb_param_h 

turb_param_z 

read_adapt_turb 

/ 

&met_param 

/ 


'Bin 5: Integrated Air concentration xy' 
' air' 

'Source 1' 

"Stuff50" 

' xy' 

.true. 

10.0 
20.0 

" 4:00 : 00 " 

" 4:00 : 00 " 

'integrated' 


'Bin 6: Instantaneous deposition' 
'dep_dry' 

'Source 1' 

"Stuff50" 

' xy' 

0.0 

"30 : 00" 

'instantaneous' 


' sigmav_simthry' 
'simthry' 

. true. 


NARAC MvM15 


stnloc_mvml5.met 

SFC 

'MtoMl' 0.0 0.0 10.0 

UPR 

'MtoMl' 0.0 0.0 


NARAC MvM15 


observ mvml5.met 


METDATASET 

UPR 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 

'MtoMl' 


' 2 000APR01 

0.5 90 

1 90 

1.5 90 

2 90 

2.5 90 

3 90 

3.5 90 

4 90 

4.5 90 

5 90 

5.5 90 

6 90 

6.5 90 


000000 ' 

2.9251 
3.4132 
3.6977 
3.8988 
4.0543 
4.1808 
4.2875 
4.3795 
4.4605 
4.5326 
4.5977 
4.6569 
4.7112 
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MtoMl' 

7 

90 

4.7613 

MtoMl' 

7.5 

90 

4.8078 

MtoMl' 

8 

90 

4.8511 

MtoMl' 

8.5 

90 

4.8918 

MtoMl' 

9 

90 

4.9299 

MtoMl' 

9.5 

90 

4.9659 

MtoMl' 

10 

90 

5 

MtoMl' 

10.5 

90 

5.0323 

MtoMl' 

11 

90 

5.063 

MtoMl' 

11.5 

90 

5.0923 

MtoMl' 

12 

90 

5.1202 

MtoMl' 

12.5 

90 

5.147 

MtoMl' 

13 

90 

5.1726 

MtoMl' 

13.5 

90 

5.1971 

MtoMl' 

14 

90 

5.2208 

MtoMl' 

14.5 

90 

5.2435 

MtoMl' 

15 

90 

5.2654 

MtoMl' 

15.5 

90 

5.2865 

MtoMl' 

16 

90 

5.3069 

MtoMl' 

16.5 

90 

5.3267 

MtoMl' 

17 

90 

5.3457 

MtoMl' 

17.5 

90 

5.3642 

MtoMl' 

18 

90 

5.3822 

MtoMl' 

18.5 

90 

5.3996 

MtoMl' 

19 

90 

5.4164 

MtoMl' 

19.5 

90 

5.4329 

MtoMl' 

20 

90 

5.4488 

MtoMl' 

20.5 

90 

5.4643 

MtoMl' 

21 

90 

5.4795 

MtoMl' 

21.5 

90 

5.4942 

MtoMl' 

22 

90 

5.5085 

MtoMl' 

22.5 

90 

5.5225 

MtoMl' 

23 

90 

5.5362 

MtoMl' 

23.5 

90 

5.5495 

MtoMl' 

24 

90 

5.5626 

MtoMl' 

24.5 

90 

5.5753 

MtoMl' 

25 

90 

5.5878 

MtoMl' 

25.5 

90 

5.5999 


NARAC YIvM 15 

observ mvm015.met (continued) 


MtoMl' 

26 

90 

5.6118 

MtoMl' 

26.5 

90 

5.6235 

MtoMl' 

27 

90 

5.6349 

MtoMl' 

27.5 

90 

5.6461 

MtoMl' 

28 

90 

5.6571 

MtoMl' 

28.5 

90 

5.6678 

MtoMl' 

29 

90 

5.6784 

MtoMl' 

29.5 

90 

5.6887 

MtoMl' 

30 

90 

5.6989 

MtoMl' 

30.5 

90 

5.7088 

MtoMl' 

31 

90 

5.7186 

MtoMl' 

31.5 

90 

5.7282 

MtoMl' 

32 

90 

5.7376 

MtoMl' 

32.5 

90 

5.7469 

MtoMl' 

33 

90 

5.756 

MtoMl' 

33.5 

90 

5.765 

MtoMl' 

34 

90 

5.7738 

MtoMl' 

34.5 

90 

5.7824 

MtoMl' 

35 

90 

5.791 

MtoMl' 

35.5 

90 

5.7993 
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MtoMl' 

36 

90 

5.8076 

MtoMl' 

36.5 

90 

5.8157 

MtoMl' 

37 

90 

5.8237 

MtoMl' 

37.5 

90 

5.8316 

MtoMl' 

38 

90 

5.8394 

MtoMl' 

38.5 

90 

5.847 

MtoMl' 

39 

90 

5.8545 

MtoMl' 

39.5 

90 

5.862 

MtoMl' 

40 

90 

5.8693 

MtoMl' 

40.5 

90 

5.8765 

MtoMl' 

41 

90 

5.8836 

MtoMl' 

41.5 

90 

5.8906 

MtoMl' 

42 

90 

5.8975 

MtoMl' 

42.5 

90 

5.9044 

MtoMl' 

43 

90 

5.9111 

MtoMl' 

43.5 

90 

5.9178 

MtoMl' 

44 

90 

5.9243 

MtoMl' 

44.5 

90 

5.9308 

MtoMl' 

45 

90 

5.9372 

MtoMl' 

45.5 

90 

5.9435 

MtoMl' 

46 

90 

5.9497 

MtoMl' 

46.5 

90 

5.9559 

MtoMl' 

47 

90 

5.962 

MtoMl' 

47.5 

90 

5.968 

MtoMl' 

48 

90 

5.9739 

MtoMl' 

48.5 

90 

5.9798 

MtoMl' 

49 

90 

5.9856 

MtoMl' 

49.5 

90 

5.9913 

MtoMl' 

50 

90 

5.997 
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HPAC MvM16 


mvml6 

HPAC 3.2 parameter 

HPAC 3.2 value 

HPAC SCREEN 

agent 

500 micron particles 

New Project Setup 

project file name 

mvml6.pii 


coordinates 

cartesian 


local oriqinx 

0 km 


local origin y 

0 km 


local origin latitude 

45 deqN 


local origin longitude 

Odeq 


reference times 

UTC 


local time of 00:00Z 

0 


mode 

standard 


dynamics 

dense qas 


static puffs 

enabled 


hazard area 

off 

New Roject Editor 

edit mode 

advanced 

Material Editor 

type 

partic le 


name 

STUFF500 


material file name 

Stuff500.mtl 


units 

kg 


bin boundaries 

see comment 


output flags 

surf, dosaqe, surf, depos. 


daytime decay rate 

0 


nighttime decay rate 

0 

Gas Parameter Editor 

density ratio (rho/riio air) 

500 kq/m~3 


deposition veloc ity 

NA 


minimum concentration 

0 

Release Editor 

time 

0 


X 

0 


y 

0 


height of release 

2 m 


source uncertainty 

no 


spec ification 

simple 


randomize location 

NA 


release duration 

instantaneous 


release rate 

NA 


agent mass 

1 kq 


xsize 

1 m 


V size 

lm 


zsize 

1 m 


MMD 

500 \i 


sigma 

2 

lime Editor 

start day 

4/1/00 


start time 

0 


stop day 

4/1/00 


stop time 

100 


maximum time step 

60 sec 


output interval 

5 min 

Domain Editor 

xmin domain 

-3 km 


xmax domain 

0.5 km 


ymin domain 

-0.1 km 


ymax domain 

0.1 km 


horizontal resolution 

default 


vertical domain max height 

500 m 


vertical domain resolution 

10 m 

Weather Editor 

weatherdata type 

surface obs. only 


boundary layer type 

observations 


large scale variability 

none 


large scale variability length scale 

NA 


large scale variability variance 

NA 


surface toughness 

0.008 m 


precipitation 

none 


obs. time bin size 

1 hr 

Meteorology Option Editor 

save meteorology fields 

no 


terrain/land cover file 

no 

Surface Observations File 

file name 

mvmN.sfc 


elevation 

10 m 


zi (boundary layer height) 

500 m 


z (altitude for wind measurement) 

10 m 


wind direction (from) 

90 deq (east) 


wind speed 

5 m/s 


MOL (Monin-Obukhov Length) 

-500 m 

Options Editor 

puff split grid level 

2 


surface resolution 

default 


puff qrid resolution 

0 


boundary layer pts 

11 


stable atm. turbulence 

le-2 m~2/s~2 


stable atm. scale 

10 m 


stable dissipation 

4e-4 m~2/s~3 


calm cond. turbulence 

0.25 m~2/s~2 


calm cond. scale 

1000 m 


surface dosage height 

0 


minimum puff mass 

le-20 


conditional averaging time 

1 hr 


lumped boundary layer 

no 
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NARAC MvM16 


adapt_mvml6.nml (same as adapt_mvm07.nml) 
lodi files mvmlS.nml 


&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/cone_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvm04/mgrid/adapt_mtom4_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 


NARAC MvM16 

lodi mvml6 . nml 


&prob_setup 

title 

tstart_str 

tstop_str 

dt_dump_str 

nbins 

nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 


'Model vs. Model: Source 16' 
"200 0APR01_0 00000" 

"200 0APR01_010000" 

"30:0" 

3 

1 

'nongauss' 

1 

"200 0APR01_0 00000" 

'arac' 


&thist_param 

/ 


Sesrc_param 
source_id 
geom_time_strs 
max_num_p art 
source_model 
src_agl_flg 

geom_type 
mean_x 
mean_y 
mean_z 
std x 


'Source 1' 

"200 0APR01_0 00000 
500000 
'neutral' 

.true. 

2 

0.0 
0.0 
2.0 
1.0 
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std_y 

std_z 

cutoff_dx_min 
cutoff_dy_min 
cutoff_dz_min 
cutoff_dx_max 
cutoff_dy_max 
cutoff_dz_max 

start_time_str = 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

mass_distrib 

mmd 

gsd 

decay_chain 


1.0 
1.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 

"200 0APR01_0 00000" 
"200 0APR01_0 00000" 
"Stuff500" 

500.0 

"200 0APR01_0 00000" 
1.0 

"mass" 

0.003 
0.0 

"lognormal" 

500.Oe-6 
2 

.false. 


NARAC MvMlo 

lodi mvml6.nml (continued) 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 
dt_bin_out_str = 
samp_type 

/ 


'Bin 1: Integrated 
'air' 

'Source 1' 

"Stuff500" 

'xy' 

.true. 

10.0 
20.0 
"30:00" 

"30:00" 

'integrated' 


ir concentration xy' 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 

/ 


'Bin 2: Integrated Air concentration xy' 
' air' 

'Source 1' 

"Stuff500" 

'xy' 

.true. 

10.0 
20.0 
" 60 : 00 " 

"60:00" 

'integrated' 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

position 

dt bin out str 


'Bin 3: Instantaneous deposition' 
' dep_dry' 

'Source 1' 

" Stuff500" 

'xy' 

0.0 

"30 : 00" 
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samp_type = 'instantaneous' 

/ 

&turb_param 

turb_param_h = 'sigmav_simthry' 

turb_param_z = 'simthry' 

read_adapt_turb = .true. 

/ 

&met_param 

/ 

NARAC MvM16 

stnloc_mvml6.met (same as stnloc_mvm07.met) 
observ mvmlS.met (same as observ mvm07.met) 
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HPAC MvM17 



B-75 

























































































NARAC MvM17 


adapt_mvml7.nml (same as adapt_mvm07.nml) 
lodi_files_mvml7.nml 

&grid_name 

num_m_grids = 1 

m_grid_name = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
c_grid_name = 

' /ModelDevelopDiag/vandv/mtom/cgrid/cone_1_40_- 

18000 0001 grd.nc' 

/ 

&metfiles 

grid_num = 1 

met_file_name = 

'/ModelDevelopDiag/vandv/mtom/mvm04/mgrid/adapt_mtom4_2000APR01_000000.n 
c' 

/ 

&particle_file 

particle_file_name = 'part.nc' 

/ 

&decay_chains_file 

decay_chains_file_name = 'none' 

/ 

NARAC \I v \I17 

lodi_mvml7.nml 

&prob_setup 
title 

tstart_str 
tstop_str 
dt_dump_str 
nbins 
nsrc 

rdm_dist 
num_met_times 
met_time_strs 
met_format 

/ 

&thist_param 

/ 

&src_param 
source_id 
geom_time_strs 
max_num_p art 
source_model 
src_agl_flg 

geom_type = 1 

xl =0.0 

yl = -250.0 

zl = 2.0 

x2 =0.0 

y2 = 250.0 


= 'Source 1' 

= "200 0APR01_0 00000" 
= 500000 
= 'neutral' 

= .true. 


= 'Model vs. Model: Source 17' 
= "200 0APR01_0 00000" 

= "200 0APR01_010000" 

= " 30 : 0 " 

= 3 
= 1 

= 'nongauss' 

= 1 

= "200 0APR01_0 00000" 

= 'arac' 
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z2 


2.0 


start_time_str 

stop_time_str 

species 

density 

er_time_strs 

emiss_rates 

er_units_type 

nset_dep_vel 

precip_coeff 

mass_distrib 

mmd 

gsd 

decay_chain 


NARAC MvM17 


"200 0APR01_0 00000" 
"200 OAPRO1_0 00000" 
"Stuff500" 

500.0 

"200 OAPRO1_0 00000" 
1.0 

"mass" 

0.003 
0.0 

"lognormal" 

500.Oe-6 
2.0 

.false. 


lodi mvml7.nml (continued) 


&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

bin_agl_flg 

position 

width 

dt_samp_str 

dt_bin_out_str 

samp_type 

/ 

&bin_param 

bin_id 

type 

source_list 

species_name 

orientation 

position 

dt_bin_out_str 

samp_type 

/ 

&turb_param 
t u rb_p a r am_h 
turb_param_z 
read_adapt_turb 


'Bin 1: Integrated Air concentration xy 
' air' 

'Source 1' 

"Stuff500" 

'xy' 

.true. 

10.0 
20.0 
"30:00" 

"30:00" 

'integrated' 


'Bin 2: Integrated Air concentration xy 
' air' 

'Source 1' 

"Stuff500" 

'xy' 

.true. 

10.0 
20.0 
" 60 : 00 " 

"60:00" 

'integrated' 


'Bin 3: instantaneous deposition' 
'dep_dry' 

'Source 1' 

"Stuff500" 

'xy' 

0.0 

"30 : 00" 

' instantaneous' 


' sigmav_simthry' 
' simthry' 

. true . 
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/ 


&met_param 

/ 

NARAC \Iv\I17 

stnloc_mvml7.met (same as stnloc_mvm07 .met) 
observ_mvml7.met (same as observ_mvm07.met) 

NARAC Gridgen Files, MvM 1-14,16,17 

gridgencgrid.nml 

&gridgen_control 

/ 

&gridgen_definition 

parent_grid_file = '/ModelDevelopDiag/vandv/mtom/grid/mvml_mgrid.nc' 
gradient_option = 'centered' 
grid_type = 'cone' 


number x points 

= 131 

number y points 

= 97 


rel x reference 

= 0.95 

rel y reference 

= 0.5 


x reference 

= 0.000 

y reference 

= 0.0 


x gridrange 

= 40.000 

y gridrange 

= 8.0 


number x basepoints = 41 

number y basepoints = 

41 

x basegrid 

= 0.0 

y basegrid 

= 

0 . 0 

x_basestep 

= 0.05 

y_basestep 

= 

0.05 


number_vert_points = 2 6 

vert_basegrid = 0.0 

vert_basestep = 0.004 

z_gridtop = 2500.0 

/ 


gridgenmgrid.nml 

&gridgen_control 

/ 

&gridgen_definition 
flat_topo = .true. 
grid_file = 'mvml_mgrid.nc' 
gradient_option = 'centered' 
grid_type = 'main' 

number_x_points = 81 number_y_points = 17 

x_reference_index = 77 y_reference_index = 9 

x_reference = 0.0 y_reference = 0.0 

x_gridrange = 40.00 y_gridrange = 8.00 

number_vert_points = 26 

vert_basegrid = 0.0 

vert_basestep = 0.004 

z_gridtop = 2500.0 

/ 
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NARAC Gridgen Files, MvM15 


gridgen_cgrid_mvml5.nml 

&gridgen_control 

/ 


&gridgen_definition 
parent_grid_file = 

'/ModelDevelopDiag/vandv/mtom/mvml5/grid/mvml5_mgrid.nc' 
gradient_option = 'centered' 
grid_type = 'cone' 


number_x_points = 117 
rel_x_reference = 0.98 
x_reference = 0.000 

x_gridrange = 100.0 


number_y_points = 137 
rel_y_reference = 0.5 
y_reference = 0.0 

y_gridrange = 20.0 


number_x_basepoints 

x_basegrid 

x_basestep 


21 number_y_basepoints = 41 

-0.5 y_basegrid = 0.0 

0.05 y_basestep = 0.05 


number_vert_points = 26 

vert_basegrid = 0.0 

vert_basestep = 0.004 

z_gridtop = 2500.0 

/ 


gridgen_mgrid_mvml5.nml 

&gridgen_control 

/ 


&gridgen_definition 
flat_topo = .true. 
grid_file = 'mvml5_mgrid.nc' 
gradient_option = 'centered' 
grid_type = 'main' 


number_x_points = 201 

x_reference_index = 197 

x_reference = 0.0 

x_gridrange = 100.0 


number_y_points = 41 

y_reference_index = 21 

y_reference = 0.0 

y_gridrange = 20.00 


number_vert_points 

vert_basegrid 

vert_basestep 

z_gridtop 

/ 


26 
0.0 
0.004 
2500.0 
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APPENDIX C 

GAS RELEASE COMPARISON PLOTS 



Log (Dosage (kg s/m )) Crosswind (km) 


APPENDIX C 

GAS RELEASE COMPARISON PLOTS 


This appendix contains dosage contours, crosswind dosage profiles, and dosage 
histograms for all of the gas releases. The 54 figures of this appendix are grouped by run 
number and time (30 and 60 minutes after the release). 


MvMOl: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure C-1. MvM 1 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 


C-1 


















































x10 s MvMOl: Crosswind Dosage Profiles, Time = 30 min. 



Crosswind Distance (km) 


Figure C-2. MvM 1 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED — ) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvMOl: Dosage Histogram Comparison, Time = 30 min. 



Figure C-3. MvM 1 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvMOl: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure C-4. MvM 1 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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x 10 -s MvMOl: Crosswind Dosage Profiles, Time = 60 min. 



Crosswind Distance (km) 


Figure C-5. MvM 1 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvMOl: Dosage Histogram Comparison, Time = 60 min. 



Figure C-6. MvM 1 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM02: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure C-7. MvM 2 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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Crosswind Distance (km) 


Figure C-8. MvM 2 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED — ) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM02: Dosage Histogram Comparison, Time = 30 min. 



Figure C-9. MvM 2 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM02: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure C-10. MvM 2 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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Frequency 


MvM02: Crosswind Dosage Profiles, Time = 60 min. 



Figure C-11. MvM 2 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM02: Dosage Histogram Comparison, Time = 60 min. 



Figure C-12. MvM 2 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM03: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure C-13. MvM 3 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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MvM03: Crosswind Dosage Profiles, Time = 30 min. 
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Figure C-14. MvM 3 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM03: Dosage Histogram Comparison, Time = 30 min. 



Figure C-15. MvM 3 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM03: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure C-16. MvM 3 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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MvM03: Crosswind Dosage Profiles, Time = 60 min. 



05 



Figure C-17. MvM 3 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM03: Dosage Histogram Comparison, Time = 60 min. 



Figure C-18. MvM 3 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM04: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure C-19. MvM 4 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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Crosswind Distance (km) 

Figure C-20. MvM 4 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED — ) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM04: Dosage Histogram Comparison, Time = 30 min. 



Figure C-21. MvM 4 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED-) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM04: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure C-22. MvM 4 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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MvM04: Crosswind Dosage Profiles, Time = 60 min. 
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Figure C-23. MvM 4 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM04: Dosage Histogram Comparison, Time = 60 min. 



Figure C-24. MvM 4 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM05: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure C-25. MvM 5 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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MvM05: Crosswind Dosage Profiles, Time = 30 min. 
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Figure C-26. MvM 5 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM05: Dosage Histogram Comparison, Time = 30 min. 



Figure C-27. MvM 5 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM05: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure C-28. MvM 5 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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MvM05: Crosswind Dosage Profiles, Time = 60 min. 
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Figure C-29. MvM 5 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM05: Dosage Histogram Comparison, Time = 60 min. 



Figure C-30. MvM 5 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM06: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure C-31. MvM 6 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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x 10 » MvM06: Crosswind Dosage Profiles, Time = 30 min. 
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Figure C-32. MvM 6 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED — ) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM06: Dosage Histogram Comparison, Time = 30 min. 



Figure C-33. MvM 6 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log (Dosage (kg s/m )) Crosswind (km) 


MvM06: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure C-34. MvM 6 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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x10 8 MvM06: Crosswind Dosage Profiles, Time = 60 min. 
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Figure C-35. MvM 6 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM06: Dosage Histogram Comparison, Time = 60 min. 



Figure C-36. MvM 6 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM07: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure C-37. MvM 7 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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MvM07: Crosswind Dosage Profiles, Time = 30 min. 
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Figure C-38. MvM 7 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM07: Dosage Histogram Comparison, Time = 30 min. 



Figure C-39. MvM 7 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM07: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure C-40. MvM 7 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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MvM07: Crosswind Dosage Profiles, Time = 60 min. 
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Figure C-41. MvM 7 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM07: Dosage Histogram Comparison, Time = 60 min. 



Figure C-42. MvM 7 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM08: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure C-43. MvM 8 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 


C-29 


















































x10 ? MvM08: Crosswind Dosage Profiles, Time = 30 min. 
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Figure C-44. MvM 8 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM08: Dosage Histogram Comparison, Time = 30 min. 



Figure C-45. MvM 8 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM08: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure C-46. MvM 8 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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x 10 ? MvM08: Crosswind Dosage Profiles, Time = 60 min. 
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Figure C-47. MvM 8 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM08: Dosage Histogram Comparison, Time = 60 min. 



Figure C-48. MvM 8 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 


C-32 











































































Log(Dosage (kg s/m )) Crosswind (km) 


MvM09: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure C-49. MvM 9 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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Figure C-50. MvM 9 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM09: Dosage Histogram Comparison, Time = 30 min. 



Figure C-51. MvM 9 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log (Dosage (kg s/m )) Crosswind (km) 


MvM09: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure C-52. MvM 9 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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MvM09: Crosswind Dosage Profiles, Time = 60 min. 
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Figure C-53. MvM 9 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM09: Dosage Histogram Comparison, Time = 60 min. 



Figure C-54. MvM 9 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 


C-36 

















































































APPENDIX D 

PARTICLE RELEASE COMPARISON PLOTS 



APPENDIX D 

PARTICLE RELEASE COMPARISON PLOTS 


This appendix contains dosage contours, crosswind dosage profiles, and dosage 
histograms for all of the particle releases. The 51 figures shown in this appendix are 
grouped by run number and time: 30 and 60 minutes after the release (120, 180, and 240 
minutes after the release for the high altitude release, MvM 15). 


MvMI 0: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure D-1. MvM 10 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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1Q s MvMI 0: Crosswind Dosage Profiles, Time = 30 min. 
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Figure D-2. MvM 10 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvMIO: Dosage Histogram Comparison, Time = 30 min. 



Figure D-3. MvM 10 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log (Dosage (kg s/m )) Crosswind (km) 


MvMIO: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure D-4. MvM 10 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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x 10 


MvMIO: Crosswind Dosage Profiles, Time = 60 min. 



Crosswind Distance (km) 

Figure D-5. MvM 10 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvMIO: Dosage Histogram Comparison, Time = 60 min. 



Figure D-6. MvM 10 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM11: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure D-7. MvM 11 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 













































x 10 -' MvM11: Crosswind Dosage Profiles, Time = 30 min. 
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Figure D-8. MvM 11 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM11: Dosage Histogram Comparison, Time = 30 min. 



Figure D-9. MvM 11 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log (Dosage (kg s/m )) Crosswind (km) 


MvM11: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure D-10. MvM 11 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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x 10 -' MvMI1: Crosswind Dosage Profiles, Time = 60 min. 



Crosswind Distance (km) 


Figure D-11. MvM 11 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED — ) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvMII: Dosage Histogram Comparison, Time = 60 min. 



Figure D-12. MvM 11 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM12: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure D-13. MvM 12 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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MvM03: Crosswind Dosage Profiles, Time = 30 min. 
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Figure D-14. MvM 12 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED — ) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM12: Dosage Histogram Comparison, Time = 30 min. 



Figure C-15. MvM 3 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log (Dosage (kg s/m )) Crosswind (km) 


MvM12: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure D-16. MvM 12 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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Figure D-17. MvM 12 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM12: Dosage Histogram Comparison, Time = 60 min. 



Figure D-18. MvM 12 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM13: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure D-19. MvM 13 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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MvM13: Crosswind Dosage Profiles, Time = 30 min. 
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Figure D-20. MvM 13 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED — ) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM13: Dosage Histogram Comparison, Time = 30 min. 



Figure D-21. MvM 13 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log (Dosage (kg s/m )) Crosswind (km) 


MvM13: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure D-22. MvM 13 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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Figure D-23. MvM 13 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM13: Dosage Histogram Comparison, Time = 60 min. 



Figure D-24. MvM 13 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM14: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure D-25. MvM 14 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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x 10 -5 MvMI 4: Crosswind Dosage Profiles, Time = 30 min. 
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Figure D-26. MvM 14 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED — ) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM14: Dosage Histogram Comparison, Time = 30 min. 



Figure D-27. MvM 14 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM14: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure D-28. MvM 14 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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MvM14: Crosswind Dosage Profiles, Time = 60 min. 



Crosswind Distance (km) 


Figure D-29. MvM 14 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED — ) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM14: Dosage Histogram Comparison, Time = 60 min. 



Figure D-30. MvM 14 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


10 


MvM15: Dosage (kg s/m 3 ), Height = 10 m, Time = 120 min. 




Figure D-31. MvM 15 at 120 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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X 10 


MvM15: Dosage Plume Widths, Time = 120 min. 
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Figure D-32. MvM 15 at 120 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM15: Dosage Histogram Comparison, Time = 120 min. 



Figure D-33. MvM 15 at 120 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM15: Dosage (kg s/m 3 ), Height = 10 m, Time = 180 min. 




Figure D-34. MvM 15 at 180 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 


D-23 



















































x10 -io MvMI5: Dosage Plume Widths, Time = 180 min. 
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Figure D-35. MvM 15 at 180 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM15: Dosage Histogram Comparison, Time = 180 min. 



Figure D-36. MvM 15 at 180 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM15: Dosage (kg s/m 3 ), Height = 10 m, Time = 240 min. 
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Figure D-37. MvM 15 at 240 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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X 10 


MvM15: Dosage Plume Widths, Time = 240 min. 



Figure D-38. MvM 15 at 240 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM15: Dosage Histogram Comparison, Time = 240 min. 



Figure D-39. MvM 15 at 240 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM16: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure D-40. MvM 16 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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MvM16: Crosswind Dosage Profiles, Time = 30 min. 
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Figure D-41. MvM 16 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED — ) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM16: Dosage Histogram Comparison, Time = 30 min. 



Figure D-42. MvM 16 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM16: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure D-43. MvM 16 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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Figure D-44. MvM 16 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED — ) 
Crosswind Dosage Plumes for Various Downwind Distances 
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Figure D-45. MvM 16 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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Log(Dosage (kg s/m )) Crosswind (km) 


MvM17: Dosage (kg s/m 3 ), Height = 10 m, Time = 30 min. 




Figure D-46. MvM 17 at 30 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 


D-31 

















































Frequency 


i 


MvM17: Crosswind Dosage Profiles, Time = 30 min. 
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Figure D-47. MvM 17 at 30 Minutes: HPAC (BLUE - -) and NARAC (RED -) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM17: Dosage Histogram Comparison, Time = 30 min. 



Figure D-48. MvM 17 at 30 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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MvM17: Dosage (kg s/m 3 ), Height = 10 m, Time = 60 min. 




Figure D-49. MvM 17 at 60 Minutes with Upper: Dosage Contours for HPAC (BLUE - -) and 
NARAC (RED — ); Lower: Dosage on Downwind Symmetry Axis 
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x 10 « MvMI 7: Crosswind Dosage Profiles, Time = 60 min. 
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Figure D-50. MvM 17 at 60 Minutes: HPAC (BLUE - -) and NARAC (RED — ) 
Crosswind Dosage Plumes for Various Downwind Distances 


MvM17: Dosage Histogram Comparison, Time = 60 min. 



Figure D-51. MvM 17 at 60 Minutes: Dosage Histograms for HPAC (BLUE - -) and 

NARAC (RED -) 
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APPENDIX E 
TASK ORDER EXTRACT 



APPENDIX E 

TASK ORDER EXTRACT 


DC-9-1797 


TITLE: Support for DTRA and LLNL in the Validation Analysis of Hazardous 
Material Transport and Dispersion Prediction Models 

This task order is for work to be performed by the Institute for Defense Analyses 
(IDA) under Contract DASW01-98-C-0067, for the Defense Threat Reduction Agency 
(DTRA) in the Department of Defense. This task order authorizes funding for FY 2000. 


1. BACKGROUND: 

The Hazard Prediction and Assessment Capability (HPAC) is a suite of codes that 
predicts the effects of hazardous material releases into the atmosphere and their impact 
on civilian and military populations. The software can use integrated source terms, high- 
resolution weather forecasts, and particulate transport models to predict hazard areas 
produced by battlefield or terrorist use of weapons of mass destruction (WMD), by 
conventional counterforce attacks against WMD facilities, or by military and industrial 
accidents. HPAC is a forward deployable, counterproliferation and counterforce 
capability software tool available for government, government-related, or academic use. 
This tool assists warfighters in selecting weapon mixes for targets containing WMD and 
in emergency response to hazardous agent release. HPAC’s relatively fast-running, 
physics-based algorithms enable users to model and predict hazard areas and human 
collateral effects in minutes. 

The DTRA Verification and Validation (V&V) Program represents ongoing 
activities performed in parallel with development of all predictive codes in support 
of HPAC. One element of V&V is to perform code-on-code comparisons. In this 
strategy, each code receives the same input. In this manner, differences in the 
output predictions can lead to the identification of software bugs, or help to assess 
technical strengths and weaknesses of component algorithms within each code. In 
addition, a certain amount of credibility for both models is achieved when their 
predictions agree. When the inputs are simple, such as for fixed winds and simple 
terrain, the predictions tend to be dominated by the dispersion algorithms. 
Comparisons at this level of complexity are important to establish fundamental 
dispersion algorithm veracity, and to help discover software bugs. As more 
complex terrain and weather is included as input, the number of physical processes 
responsible for transport and dispersion increases and the predictions become the 
result of many interdependent algorithm calculations. 
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Code-on-code comparisons will be performed using the DTRA code HP AC, 
the Lawrence Livermore National Laboratory (LLNL) code Atmospheric Release 
Advisory Capability (ARAC), and, possibly, Sandia-developed codes. These codes 
represent major national investments in transport and dispersion modeling within 
their respective applications. The comparisons will provide information from 
which to validate the HP AC and ARAC models (and perhaps others), as well as 
provide an opportunity to advance both technologies. The code comparisons will 
include short, medium, and long-range transport distances. Complex terrain and 
weather will also be included. 

It is very difficult to separate meteorological uncertainty from the transport 
and dispersion model accuracy when comparing predictions to field-trial validation 
quality or real-world data. The validation challenge is to assess whether a model 
performs well over different field trials, and ultimately reflects real-world 
phenomena. Some codes perform better under certain conditions and specific 
scenarios. Hazard prediction models are generally developed for a range of user 
communities and applications. Each user community has a different set of 
requirements. Thus, the corresponding hazard models tend to be optimized for 
specific applications. The process of accrediting a model is always couched in 
terms of the end-user requirements. 

Various figures-of-merit (FOM) are used to express model performance 
relative to observed data. Most FOMs tend to use manifestations of a ratio 
(geometric or arithmetic) between the predicted and measured quantities. The 
compared quantities are usually peak, plume-centerline, and off-axis concentration 
or dosage, as well as crosswind and along-wind spread and area coverage. Other 
FOMs may include the second-moment of the dosage and concentration values at a 
sampler location. All these FOMs are reasonable validation performance measures, 
but none of them explicitly expresses an application-oriented performance measure. 
A “yardstick” is needed that measures application-oriented model performance. 
The scale on this yardstick would clearly and directly relate to specific user’s 
concerns and needs. The pursuit of this “accreditation” performance measure is a 
new initiative at DTRA. 


2. OBJECTIVE: 

IDA will conduct independent analysis and special studies associated with 
verification and validation of the suite of models associated with the Hazard 
Assessment and Prediction Capability. IDA will support development of user- 
oriented performance measures of effectiveness (MOE) using validation quality 
field trial data sets; coordinate scenario definition and arbitration for code-on-code 
V&V activities; and assist DTRA and the Department of Energy in identifying the 
V&V parameter space associated with various hazard assessment and collateral 
effects communities. 

The objectives of verification and validation analysis and coordination are: 
(1) to ensure that a consistent analysis approach is used when comparing model 
predictions, and assist DTRA in the implementation of code-on-code analysis. 
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comparisons, and interpretation; and (2) to define measures of effectiveness in 
terms of user-specific objectives and applications. 

The scope of this effort may be expanded to other programs as directed by 

DTRA. 

3. STATEMENT OF WORK: 

As required by DTRA technical representatives, IDA will perform the following 

tasks: 


a. Support the planning, implementation, arbitration, and evaluation of code-on- 
code comparison activities. The purpose of these activities is to compare the transport 
and dispersion algorithms and corresponding output predictions between DTRA’s suite of 
models, the DOE ARAC’s suite of transport and dispersion models, and other models as 
called for. IDA will support code-on-code scenario definition, coordinate the 
identification and implementation of common performance measures, and support 
development of a common analysis approach. IDA will conduct independent analysis, as 
needed, to support the code-on-code analysis and interpretation of the results. 

b. Explore validation and accreditation MOEs given a framework that includes 
quantification of false positive and false negative predictions. This exploration would 
include the computation of MOE values for various formulations based on short-range 
comparisons of HP AC and ARAC predictions to field trial data. A key to interpreting the 
results of this effort will be obtaining a sense for what are the acceptable user 
requirements. These requirements will differ among potential user groups (military 
targeting, passive CB defense, civilian first responders, military versus civilian 
population human effects, etc.). 


4. CORE STATEMENT : 

This research is consistent with IDA’s mission in that it will support specific 
analytical requirements of the sponsor and will assist the sponsor with planning efforts. 
Accomplishment of this task order requires an organization with experience in 
operationally oriented issues from a joint and combined perspective, which IDA, a 
Federally Funded Research and Development Center, is able to provide. It draws upon 
IDA’s core competencies in Systems Evaluations and Operational Test and Evaluation. 
Performance of this task order will benefit from and contribute to the long-term 
continuity of IDA’s research program. 


E-3 




REPORT DOCUMENTATION PAGE 

Form Approved 

OMB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data 
sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any 
other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations 
and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22203-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704- 
0188), Washington, DC 20503. 

1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. 1 

May 2001 

REPORT TYPE AND DATES COVERED 

Pinal: Jan. 2000 - May 2001 

4. TITLE AND SUBTITLE 

Evaluation of Transport and Dispersion Models: A Controlled Comparison of 

HP AC and NARAC Predictions 

5. FUNDING NUMBERS 

DASW01-98-C-0067, OUSD(A) 

DC-9-1797 

6. AUTHOR(s) 

Dr. Steve Warner, Dr. Nathan Platt, Dr. James L. Heagy, et al. 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Institute for Defense Analyses 

1801 N. Beauregard St. 

Alexandria, VA 22311 

8. PERFORMING ORGANIZATION REPORT NUMBER 

IDA Paper P-3555 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

Dr. Allan Reiter 

Defense Threat Reduction Agency (DTRA) 

6801 Telegraph Road 

Alexandria, VA 22310-3398 

10. SPONSORING/MONITORING AGENCY REPORT 
NUMBER 

11. SUPPLEMENTARY NOTES 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Approved for public release; distribution unlimited. Directorate for Lreedom of 
Information and Security Review, 23 May 2001. 

12b. DISTRIBUTION CODE 

13. ABSTRACT ( Maximum 200 words) 


The goal of this task is to improve the verification, validation, and accreditation (VV&A) of hazard prediction and assessment 
models. These studies are part of a larger joint VV&A effort that DTRA and the Department of Energy, via the Lawrence 
Livermore National Laboratory (LLNL), are conducting. This joint effort includes comparisons of the LLNL and DTRA transport 
and dispersion modeling systems - National Atmospheric Release Advisory Center (NARAC) and Hazard Prediction and 
Assessment Capability (HPAC), respectively. The Institute for Defense Analyses’ role includes conducting comparisons between 
the models and providing analysis and discussions associated with these examinations. This paper provides detailed comparisons 
of HP AC and NARAC predictions for 17 releases. The 17 trials were chosen with the desire for simplicity, with the expectation 
that simple releases would test only the core transport and dispersion algorithms within each model. Consequently, all releases 
were chosen to have idealized weather - a single unidirectional wind observation of 5 meters per second at 10 meters altitude - and 
only flat terrain was considered. The 17 trials did allow for the comparison of model predictions as a function of atmospheric 
stability condition, source type (gas or particle), release duration, release altitude, and particle size. 

14. SUBJECT TERMS 15. NUMBER OF PAGES 

Model validation; model intercomparison, hazardous material transport and dispersion; HP AC; 276 

NARAC; measure of effectiveness 16. price code 


17. SECURITY CLASSIFICATION OF 

18. SECURITY CLASSIFICATION OF 

19. SECURITY CLASSIFICATION OF 

20. LIMITATION OF 

REPORT 

THIS PAGE 

ABSTRACT 

ABSTRACT 

Unclassified 

Unclassified 

Unclassified 

Unlimited 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 




